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OPDOT: A COMPUTER PROGRAM FOR THE OPTIMUM PRELIMINARY 

DESIGN OF A TRANSPORT AIRPLANE 


Steven M. Sliva and P, Douglas Arbuckle 


ABSTRACT 


A description of a computer program, OPDOT, for the optimal preliminary 
design of transport aircraft is given. OPDOT utilizes constrained parameter 
optimization to minimize a performance index (e.g. direct operating cost per 
block hour) while satisfying operating constraints. The approach in OPDOT 
uses geometric descriptors as independent design variables. The independent 
design variables are systematically iterated to find the optimum design. The 
technical development of the program is provided and a program listing with 
sample input and output are utilized to illustrate its use in preliminary 
design. This is not meant to be a user’s guide, but rather a description of a 
useful design tool developed for studying the application of new technologies 
to transport airplanes. 


INTRODUCTION 


When new technologies in aircraft design, fabrication and operation are 
evaluated, current practice requires engineering judgment in making compromises. 
An engineer utilizes a combination of limited analyses , experience and intuition 
to combine new technologies (e.g., aerodynamics, controls or structures) to 
maximize the benefits. This approach is imprecise because it involves extra- 
polating experience from previous designs and because the improvements are 
usually made to augment multiple, ill-defined criteria (e.g., weight, cost, or 
performance). To properly evaluate any changes in design concepts, the 
airplane configuration should he allowed to evolve to optimize a single, well- 
defined performance index. 

This report describes OPDOT (Optimum Preliminary Design of Transports), a 
computer program written to perform preliminary design and evaluation of trans- 
port aircraft using nonlinear programing techniques . A set of independent 
design variables is iterated upon until a minimum of a performance index which 
satisfies a series of constraint functions has been calculated. The design 
variables usually consist of geometry characteristics and mission parameters, 
while the constraint functions include, for example, regulatory performance 
requirements and handling quality design criteria. A slightly modified, 
constrained sequential optimizer is utilized in the program. 

This program, therefore, allows the evaluation of new technologies 
incorporated into an aircraft design in an optimal fashion. The degree of 
detail in the analyses when the performance function and the constraint func- 
tions are evaluated is at the preliminary design or classical aeronautics level . 
That is, the precision in some phases of the calculations is expected to be as 



poor as 5-10 percent. Hence, whereas the predictive capabilities are expected 
to be marginal, the accuracy of the relative comparisons of designs is expected 
to be good. 

This report, which includes a program listing, sample input and sample 
output, is a description of a useful analytical tool for analyzing the effects 
of new technologies on the preliminary design and sizing of transport airplanes. 
It describes the methods of calculation, program organization and some of the 
various options available, but it is not meant to be a comprehensive user's 
manual. The program code was written to expeditiously obtain answers for a 
study of the impact of active controls upon transport design. This limited the 
amount of effort that could be spent on developing user flexibility and on 
integrating into the program a high degree of self-annotation. 


SYMBOLS 


A wing aspect ratio 

A horizontal tail aspect ratio 

b span, m 

b ratio of tail 

B Brequet range 

c specific fuel 

c mean aerodynamic chord of ving 

c ith constraint value 
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C chordvise force coefficient 

drag coefficient, drag/qS 

C lift coefficient, lift/qS 
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C approach lift coefficient 

C tail lift coefficient 
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design lift coefficient for ving airfoil section 
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lift coefficient at takeoff 

lift coefficient during second segment climb 

pitching moment coefficient, pitching moment/qSc 

normal force coefficient 

thrust coefficient 

unaugmented performance index 
augmented performance index 
ith constraint function 

vertical tail height, m 

angle of incidence , deg 

penalizing weight 

pitch feedback gain, 36/30 

pitch rate feedback gain, 36/3- 

0 

distance from wing aerodynamic center to tail aerodynamic center, m 

glide ratio 
mass (slugs) 

Mach number 

pitching moment derivative due to pitching velocity, sec*^ 
pitching moment derivative due to unit vertical velocity, sec”^ 
pitching moment derivative due to vertical acceleration, sec"^ 
pitching moment derivative due to control deflection, sec"^ 
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number of engines 
airplane vertical gain, g’ s/rad 

ith constraint normalizing factor, 

design parameter (constant for optimization) 

dynamic pressure 

airplane range, kilometers 

lower and upper boundaries of ith constraint 

tail area, 

2 

wing area, m 

thrust at altitude, N 
installed thrust at sea level, N 

Take-off parameter 
reference velocity, m/sec 

cruise velocity, m/sec 

lower and upper boundaries of ith independent design variable 

aircraft weight at altitude, N 
maximum aircraft weight at take-off, U 

initial aircraft weight during cruise segment, N 

final aircraft weight after cruise segment, N 

ith independent design variable 

% MAC from datum to c.g. in x-direction 




independent design variable transformed 


ac 


dimensional vertical force derivative due to elevator deflection. 

-1 

sec 

% MAC from datum to c.g, in z— direction 

height of thrust vector from c.g. MAC) 

Oswald's efficiency factor for 0 ^ C < C 

L L 


Oswald's efficiency factor for C > C^ 

L L 


0 ) 


sp 


flight path angle, rad 
short period frequency, sec 

Munk's interference factor 

\ 

tail efficiency, — 


-1 


Superscript 
* optimum 


Subscripts 

ac aerodynamic center 

cr cruise 

q. pitch rate 

t tail 

u velocity 

w wing 

a angle of attack 

a rate of change in angle of attack 
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elevator deflection 


PROGRAM DESCRIPTION 
General 


The overall flow of the program is depicted in figure 1, A set of 
starting values for the selected independent design variables and design 
constants is input and used to initialize the optimizer and the data base. 

Initially, the program was written with seven independent variables (wing area, 
wing aspect ratio, fuselage length, horizontal tail area, horizontal tail aspect 
ratio, aft-most center-of-gravity position, and installed thrust), but it has 
the inherent capability to handle more and has successfully converged with 
thirteen. Typical design constants include nonvarying geometries, mission 
parameters, economic constants, nonlinear aerodynamics data and some levels of 
technology. An extensive list of design constants that were used in one study 
is shown in Appendix VII. 

Design constants are prime candidates for being changed to independent 
design variables. Both design constants and independent design variables are 
held constant for each call to the performance function evaluation routines. 

Independent design variables are typically altered each function call by the 
optimizer, while design constants are not allowed to vary for the entire opti- 
mization. A method for augmenting the set of independent design variables with 
design constants will be described in a later section. 

The inputs (the current value of independent design variables and the 
design constants) are utilized by a sequence of subroutines that calculate a 
performance index which is selected by the user. Typically, minimum direct 
operating cost per block hour is chosen, but minimum direct operating cost per 
flight, maximum return-on-investment per year, minimtim income required for a 
15 percent return on investment, maximum L/D and minimum take-off gross weight 
are also available as criteria to be optimized. During the series of subroutine 
calls, data is exchanged with and stored in the data base for future use. The 
program has been constructed in a modular fashion to allow users to replace 
routines with preferred versions to allow significant configuration changes or 
to improve the level of accuracy. 

Next a series of subroutines is called to calculate the constraint 
functions. Those that are calculated for cruising flight utilize data stored in 
the data base during the performance function evaluations. Many subroutines a 

were written in such a fashion as to provide data in slow flight configurations 
as well as in cruising flight. These are called to yield take-off and landing 
performance data. As a byproduct, the longitudinal stability derivatives are • 

generated. These nondimensional derivatives, for both approach and cruise, are 
converted to dimensional derivatives and are then used to determine the roots of 
a fourth order linear model of the longitudinal dynamics. These roots are used 
to calculate the damping and frequency in the short period and phugoid modes. 
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The program determines which constraint functions are violated and adds a 
penalty term for each violation to the performance index to create an augmented 
performance index. The optimizer then iterates upon the design variables to 
minimize the augmented function. If the weights on the penalty terms are 
sufficiently large, the violations will be driven to zero. A convergence of the 
optimizer resxilts in the minimum unaugmented performance index that satisfies 
the constraint functions. 


Optimization Code 


The optimization is performed by a sequential simplex method (Ref. 1 and 2) 
which utilizes a continuous penalty function. This direct search algorithm has 
the advantage of not using gradient evaluations, and hence does not perform 
poorly near "ridges" in the performance index. Additionally, the penalty scheme 
is independent of the number of active constraints. Its chief disadvantage is 
slow convergence in large regions of small gradients of the augmented 
performance function with respect to the independent design variables. 

The general problem is formulated as follows: 

Let the unaugmented performance index, f, be a function of the independent 
design variables, x, and design parameters, p. 


and 


then 


f = fcn(x,p) 
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( 1 ) 


(2) 


The goal is to' find the minimum of the augmented function, F, with the gains, K, 
large . 


A variable transformation (Ref. 3) is used to automatically scale the 
variables and apply "side" constraints, which are inequality constraints applied 
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directly on the design variables. This resvilted in a reduction in the number 
of iterations required for convergence. 

The form of the transformation is as follows: 


X, 

1 



2 h 



(3) 


where V and V are the ith upper and lower independent design variable 

XJ. , A» , 

1 1 

boundaries. So the simplex optimizer iterates on the transformed variable Z, 
which spans the set of allowable values of the independent design variables with 
the range in Z of 1 to -1. This allows consistency in step size selection and 
limits the allowed values of the independent design variables. 

A version of the program is listed in Appendix I. The main program 
SIMPACT, the subroutines NELMIN and SETUP with the function FN are used to 
perform the optimization. Some key variables and a description of the pertinent 
labelled common blocks are shown in Appendices V and VI, respectively. Prior to 
the optimization, a series of inputs to initialize the optimization blocks is 
read in and XINPUT is used to initialize the aircraft data. NELMIN, the 
subroutine which returns the constrained minimxim, is called several times 
(usually two) with increasing weights and diminishing convergence criteria and 
initial step sizes. This is to help in obtaining a satisfactory local minimum 
with no constraints violated and, ideally, with the active constraints resting 
against their boundaries. 

NELMIN calls FN which returns the augmented performance index. FN calls 
SETUP which performs the variable transformations, obtains the unaugmented 
performance index, calls the constraint evaluation routines, determines the 
penalty terms and then assembles the augmented performance index. The unaug- 
mented performance index is determined by calling DOCOST and the constraint 
functions are calculated from CNSTRN. 


Evaluation of Unaugmented Performance Index 


A flow diagram showing the general procedure for evaluating the unaugmented 
performance index is shown in figure 2. DOCOST (included in Appendix I) 
assembles various cost components by first calling GEO to calculate and store 
some geometry constants and then calling CGCAL to assign the center-of-gravity 
positions for the various phases of flight and the landing gear position (if 
variable). Then WEIGHT is called which is used to estimate the airplane's 
operating weights, the amount of fuel burned during the mission and a variety of 
other parameters required from the cruise portion of the flight. 

In WEIGHT, an initial estimate of take-off weight and fuel fraction is made. 

The individual weight components are determined using statistical relationships 
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from references through 6. The primary source was reference 4, hut the 
critical components for the intended uses of the program (i.e,, wing, horizontal 
tail and fuselage) were limited to geometric ranges to maintain validity. To 
improve the capability of predicting the weights of these components (e.g., at 
high aspect ratios) an average of values calculated from references U through 6 
was made. After the component weights are summed, FUELCAL is used to determine 
the weight of the fuel required to fly the passenger mission and the reserve 
mission. This fuel weight is used to estimate the weight of the fuel systems. 

The sum of the individual estimated weight components is compared with the 
initial estimate of take-off weight; and, if the difference is greater than some 
convergence criterion (usually about .2 Newtons), a new estimate is made and the 
components are summed again. This continues until the weight loop convergence 
criterion is satisfied. The new estimate for the gross take-off weight is made 
through a weighting scheme based on the number of current iterations. The 
total and average number of iterations is displayed to the user to provide 
guidance in possible programing changes in the event of slow weight loop conver- 
gence. Usually WEIGHT averages between 3 and 5 iterations per function call 
during an optimization run. 

FUELCAL assumes a flight profile schematically illustrated in Figure 3. A 

fixed percentage of the total fuel burnoff is attributed to the following tasks: 

taxi, take-off, initial climb, climb to cruise, descent and landing. The 

remainder of the flight (the cruise portion) is divided into ten equal segments. 

During the first segment the transport is flown at a C_ for maximum range 

L 

factor, B*. The initial cruise altitude is 11000 m. (36000 ft), and CRUALT is 
called to find the desired altitude at the end of the first segment to maintain 
the same for a new weight, while insuring the aircraft is also cruising at 

the desired Mdch Number. The required excess thrust to generate the calculated 
climb gradient is then saved for future use in the constraint functions. 

Segments 2 through 5 are flown in a cruise/climb mode at and B*, which 

can be calculated from classical relationships. At segment 6, however, the 
climb is increased so that segments 7 through 10 can be flown at cruise Mach 
number, M , and 9^% L/D . The cruise is backed off L/D slightly to help 

provide some speed stability. 

Thus, as modelled above, the independent design variables only impact the 
cruise portion of the flight. To simulate the complex reserve mission require- 
ment, the transport is flown for an additional lUOO kilometers (lOOO nautical 
miles) at 9100 meters (30,000 feet) at the speed for maximum range. 

CRUFUEL calculates the amount of fuel burned during each segment as well as 
the time required to fly it and the altitude change to satisfy the cruise/climb 
assumptions. As previously described, the aircraft is flown at the speed for 
maximum Brequet range factor during the first five segments provided the 
resulting Mach number is less than or equal to the desired cruise Mach number. 
The solution comes from classical aeronautics, for example, reference U. 
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L/Dg^ = .9h3 L/D 


ik) 


from aeronautics and assuming parabolic drag polars 



.79 C, 


-‘L/D 


max 


(5) 


CRUALT returns the required altitude to fly at the specified weight, lift 
coefficient and Mach nurnher at the end of each segment. CRUFUEL then estimates 
a rate-of-climb slightly greater than that which would maintain the maximum 
range factor cruise for the given altitudes. The eventual goal is to achieve 
a cruise at 98 ^ of maximum L/D for the last four segments of the cruise distance 
at the cruise Mach number. Holding the Mach number fixed results in increasing 
lift coefficients as altitude increases. This is continued until the airplane 
attains maximum L/D. 

The assumed mission profile, although patently suboptimal, varies less 
than 3 percent in fuel consumption from some optimal profiles (Ref. T)* Given 
the level of accuracy of the program and the desire to compare designs rather 
than predict the performance of one design, this level of precision was deemed 
acceptable. 

XLOD is used to estimate the aerodynamic performance of the airplane. The 
parasite drag is obtained from CDZL. CDZL performs a drag buildup by estimating 
the Reynolds number, friction coefficients, and various nonlinear constants as 
illustrated in references k and 8. Increments in drag are included for ”crud" 
drag and flap deflections. XLOD then calls STABCOD to estimate the stability 
and control derivatives while in the indicated flight configuration. These 
nondimensional derivatives are obtained from a combination of empirical and 
analytical relations developed from references 8 through 10 for transport 
airplanes. Some aeroelastic correction factors are applied to the derivatives 
based on observations of data in references 10 through 12. 

XLOD then utilizes the stability and control data as it calls TRIM. The 
desired airplane lift coefficient with the specified Mach number, parasite drag, 
flap configuration, center-of- gravity position and phase of flight are input to 
TRIM. 


The following classical non-linear trim equations (Ref. 13) were used in 
TRIM to represent the normal and chordwise forces and to solve for the required 
tail or wing lift coefficients: 
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C = C cos (a - 1 ) - C_ sin (a - i ) 
CD V w L V V 

w V w 


(7) 


C. - C cos (a - 1 ) + C- sin (a - i ) 
N L w w D V V 

w w w 


( 8 ) 


An iterative scheme utilizing the above equations is used whereby a new 
tail lift coefficient is estimated until a convergence criterion is satisfied. 
Direct substitution into the vertical force, horizontal force and pitching 
moment equation is not possible since it has been deemed inappropriate to 
linearize the transcendental functions. It would also require neglecting the 
vertical offset of the center-of-gravity from the aerodynamic center and thrust- 
line and neglecting the contributions due to tail drag. Typically three or four 
iterations are required to satisfy the trim convergence criteria (AC < .003). 

^t “ 

TRIM is used in one of two fashions. First, if a desired airplane lift 
coefficient is input, the routine iterates to find the required lift coeffi- 
cients for the tail and wing. Alternatively, if a wing lift coefficient is 
input, the required tail lift coefficient is output along with the resulting 
airplane lift coefficient. The latter mode is used to determine the maximum 
trimmed lift coefficient for approach or take-off configurations where stalling 
of the wing is a concern. 

The wing compressibility drag contribution is calculated in XLOD by using 

the empirical relationships found in reference li+, which were derived from 

supercritical aerodynamics wind tunnel data. The fuselage compressibility drag 

term is modelled from the graphs in reference 13. It should be noted that it is 

assumed that the fuselage is not area ruled and hence calculated drag will be 

pessimistic for transonic configurations (l.O > M > 0.9). 

cr 

The induced drag contribution is obtained as follows: 
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The first two terms are the wing contribution including an offset for the 
design lift coefficient of the highly cambered wing. The third term represents 
the interference drag between the lift vectors of the tail and the wing. Notice 
how the interference term could be negative if the tail lift were downward. The 
fourth term is the drag contribution of the tail lift (positive for a tail load 
in any direction). The interference factor, a, is a function of the gap ratio, 
h^/b^, and the span ratio, b. This term is calculated from a least squares 

polynomial fit (Ref. 15) of the curves in reference l6. 

The total drag, calculated in XLOD, is the sum of the induced drag, the 
drag due to elevator deflection (C estimated from Ref. 13), the compress- 

ibility drag and the parasite drag. The L/D is obviously calculated as C /C . 

li D 

Additionally, the lift coefficient for L/D is estimated and stored for future 

max 

use. XLOD, CDZL, STABCOD and TRIM are generalized to function for both cruise 
and approach conditions. 

CRUFUEL then calls ENGINE to determine the thrust and specific fuel 
consumption as a function of altitude and Mach number. The engine performance 
comes from a normalized model of the baseline engine from reference 7* The 
engine weight and size are scaled according to reference 4 based upon the 
installed thrust. The specific fuel consumption obtained from ENGINE and the 
L/D from XLOD are substituted into the classical Brequet range relationship for 
each cruise segment to determine the fuel consumption. 



exp 


c R ' 

V L/D 
cr 


( 10 ) 


After WEIGHT has converged upon the aircraft operating weights for the 
desired mission, DOCOST continues with the cost estimates. AIRCOST used the 
weight, some production assumptions (number of prototypes, number of production, 
time for development, etc.) and the statistical relationships of reference 4 to 
predict the purchase cost of the airplane. Some cost increases based on 
references 17 and l8 are arbitrarily applied to account for the inclusion of 
active controls. 

MAINCST uses statistical relationships found in references 19 and 20 to 
determine the cost of airplane maintenance. A mnnber of configuration assump- 
tions have to be made (e.g., number of APU*s, windows and IMU’s) to utilize 
these equations (see Appendix VII ). .The equations for estimating the other 
direct operating cost terms come from references 17 and 20 . Indirect operating 
cost is predicted using the statistical, relationships from references 17 and 21. 
An annual rate of return on investment (ROI) is calculated and the remaining 
performance indices are saved in the data base for future use by the optimizer. 



EVALUATION OF CONSTRAINT FUNCTIONS 


The program version included herein has 52 constraint functions that can be 
applied to the transport design. The designer chooses an upper and lower bound- 
ary for each function as an input. The program does a test on all constraint 
lower boundaries; and, if -999 is input for the lower boundary of a constraint 
function, the constraint is not included in the penalty function even if it is a 
violation. The constraint functions are of two general types, design or opera- 
tional constraints and handling quality constraints . The first set restricts 
the design to avoid infeasible geometries or to insure satisfying performance 
regulations and mission requirements. The second set is used in the study of 
tail sizing and the impact of flying qualities design criteria upon transports 
with relaxed static stability augmentation systems. 

CNSTRN returns the values of the constraint functions to SETUP, where they 
are identified as violated or not violated, normalized and assembled into a 
penalty function. The ratio of cruise thrust available to cruise thrust 
required is obtained from the data base as are the cruise altitudes and the 
cruise wing lift coefficient. The geometry constraints include insuring that 
the aft center-of-gravity is far enough forward of the main landing gear to 
provide sufficient nose wheel steering and that there is enough floor space to 
seat the passengers. 

The missed approach climb gradient and the second segment climb gradient 
are engine-out performance requirements specified by the Federal Aviation 
Regulations, FAR's, (Ref. 22). The required thrust to weight ratio is calcu- 
lated as follows : 




1 

L/D 


+ sin a 



( 11 ) 


The flight path angle is specified by the FAR's and the L/D is obtained by 
calling XLOD with the proper speed and configuration specified. The second 
segment climb is performed at maximum gross weight and at a lift coefficient 
defined by 




/1.1+U 


(12) 


The missed approach climb is performed at maximum landing weight and at a lift 
coefficient defined by 




/1. 69 


(13) 
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C and C are determined by specifying the maximtam lift coefficient that the 
^TO ^so 

wing can support in each flap configuration and then calling XLOD, which for 
this case trims the airplane maintaining the wing lift coefficient. Since the 
tail of conventional configurations is generally carrying a download at this 
point, the aircraft will usually trim at an overall lift coefficient less than 
the one specified for the wing alone. 

The landing and take-off field length are determined using empirical 
relationships from reference 23. The lauding field length utilizes approach 
speed as the independent parameter. TOP, which is defined as 


TOP 


»TO''^W 

T /W 


(lit) 


is used as the independent parameter for the take-off analysis. 

Several of the flying quality constraints are control power requirements. 
One is to maintain a lift coefficient on the tail greater than -.8 during 
approach (Ref. 24). This is to provide adequate margin from the maximum down- 
load capable of being supported by the tail (generally C = —1.2) to insure 

t 

max 

a capability to rotate and trim the aircraft for landing. 

The tail is also required to be able to rotate the airplane for take-off. 
The maximum available download, the tail can produce during take-off roll is 
calculated using the relationships in reference l4 modified for ground effect 
using the geometric angle-of-attack method of reference 8. The required down- 
load at the tail is determined from statics , such as the development in 
reference 25. The constraint specifically requires the ratio of the available 
tail download to the required tail download to be greater than 1. 

The flying quality analysis is initiated by trimming the airplane in 
approach configuration with an altitude of 150 meters by calling XLOD. The 
nondimensional stability derivatives for cruise and approach, which are stored 
in the data base, are converted to dimensional stability derivatives by DIMDER. 
The characteristic equation for the fourth order longitudinal set of equations 
(Ref. 26) is assembled by LONGRT. The four roots are determined by using RPOLY, 
a system routine for finding roots of polynomials on Langley Research Center's 
FORTRAN Math Library. 

The preceding analysis is used to assign the following constraint functions 
for both cruise and approach: static stability, maneuver stability, dynsmic 

stability, phugoid mode frequency and damping and the short period mode 
frequency and damping. The dimensional stability derivatives are used to 
estimate the following parameters which have been suggested as usefiil for flying 



qualities analysis: time-to-double , time-to-half , 

2 

approach, vertical gain and 


flight path stability in 


The tail is configured with a trimmable stabilizer, maintaining the 
elevator for maneuvering. If the stabilizer "hits" a control stop in either 
cruise or approach, the elevator is deflected to satisfy the remaining trim 
requirements. The amount of this trim deflection is stored as a constraint 
function and is usually required to be zero. Otherwise, a control deflection 
would indicate a loss of control authority, and in some cases, an increase in 
trim drag. 


Since one intended use of the program is to study unaugmented flying 
qualities design criteria, it is desirable to insure that the airplane is 
capable of being practically augmented to excellent flying qualities. A pitch- 
attitude-hold with pitch-rate-command autopilot was chosen as a conservative 
estimate of an augmentation system. The airplane is arbitrarily augmented to 
to 2 

have: = 1 and C = .7. An extension of reference 27 is used to calculate 

n^/a sp 

the feedback gains and K*. In reference 27 it is assumed that and 

are negligible and hence zero. If these assumptions are removed, the following 
relations are derived utilizing the short period approximation to the 
longitudinal dynamics : 
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These gains are then substituted in equations B-31 and B -38 of reference 27 
for estimating the variance of the elevator position and elevator position rate 
in cruise and approach. The turbulence is assumed to have a characteristic 
length of 760 meters (2500 feet) with an RMS gust level of .9 and 2.13 m/sec 
(3 and 7 ft/sec) in cruise and approach, respectively. These autopilot calcul- 
ations are used to assign the following quantities in cruise and approach to 

available constraint functions: K_, K*, a. and a*. The constraint functions 

u o o 0 

are used to insure that enough aerodynamic control exists to stabilize the 
airplane to excellent flying qualities and that enough hydraulic capability is 
available to prevent control surface rate saturation in heavy turbulence. 
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PROGRAM USE 


A listing of the computer program set up to optimize seven design 
variables is included as Appendix I. Appendices II and III show sample input 
and output, respectively, for the program. Appendix IV contains a listing of a 
procedure file that will execute the program on the Langley Research Center 
computer system. As an aid in understanding the coding, a list of key program 
variables by routine and descriptions of their values are presented in 
Appendix V. Appendix VI is a compendium describing the variables in the 
common blocks. 

The procedure file listed in Appendix IV contains a call to PPB, a program 
for executing a geometry preprocessor upon the output data placed on TAPEU by 
subroutine XOUTPUT. This preprocessor puts on TAPE? a data set suitable for 
executing ABS2290, an airplane graphics package described in reference 28. It 
is useful during conceptual design trade studies to see pictures of the 
configurations being generated. An example of this feature is shown in 
Figure 4. 

Typically, with a case similar to the one contained in the appendices, 
approximately 500 function calls, or iterations, are required for a convergence 
of NELMIN. A function call averages about 1 second in execution time on the 
Langley Research Center Cyber 1T5* 

If the user desired to add more design variables for the optimizer to 
iterate upon, these can be added as assignment statements beneath the transform- 
ation in SETUP (see Appendix l). Sample statements are left for adding cruise 
Mach number, wing sweep angle, wing thickness ratio and fuselage diameter as 
design variables. Usually all that is necessary to add an independent design 
variable is to equate it to a variable in the system of common blocks, which 
should contain degrees of freedom adequate for studies at the preliminary design 
level. 

An array in the common block GEOM named PX has been included to aid in the 
study of certain changes representative of technological improvements. The 
specifics of its use are described in Appendix VI. For example, the following 
parameters could be studied during a design series: engine fuel efficiency, 

wing drag reduction, pitching moment reduction and structural efficiency. 
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APPENDIX I - PROGRAM LISTING 


PROGRAM SIMPACT ( INPUT , OUTPUT ,TAPE5=INPUT , TAPE6=OUTPUT ,TAPE4 ) 
************************************** 

* OPTIMIZATION COMMON BLOCKS * 
************************************** 

COMMON /AVOID/FACT (59) ,GNORM(59) 

COMMON /CONSTR/SU(59) ,SL(59) ,XINEQ(59) 

COMMON /DEBUG/IDEBUG,IDEBUG2 

COMMCW /FCOUNT/ICNT, JCNT,KCNT,NDAV,LCNT, IWT,NOIT 
COMMON /LABELP/ARUN(8) 

COMMON /PERFAINAUG,SCF,NVAR,MINEQ 
COMMON /STRAIN/ CON (59) 

COMMON /VARIAB/AMP(15),AVE(15) 

COMMCW AlOL/CAYY,MC,ILINE,IOUT 
******************************** 

* DESIGN COMMON BLOCKS * 

******************************** 

COMMON /DEVARA>ESIGN(15) ,ITERM(10) ,CST(10) 

COMMON /ORA3/CDS (6 ) , CDSAP (6 ) 

COMMOJ /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /GRAVIIY/CG(6) 

COMMOI /STAB/DERIVCR(15) ,DERIVAP(15) ,STOR(20) 

COMMON /WrSVE/WTS(20) 

********************** 

* INPUT DATA * 

********************** 

REAL STEP(15),XMIN(15),XSEC(15) 

REAL XBAR(15) ,XBARO(15) 

REAL XL (15) ,XU(15) 

READ(5,18) ARUN 

18 FORMAT (8 AlO) 

READ (5,*) NVAR,MINEQ 
READ(5,*) (XBARO(I) ,I=1,NVAR) 

IDEBUG=0 

IWT=0 

IDEBUG2=0 

DO 19 J=1,NVAR 

READ(5,*) JI,XL(J),XU(J) 

AVE ( J ) = (XL ( J ) +XU ( J ) ) /2 . 0 
AMP (J) = (XU (J) -XL (J) ) /2 .0 

19 CONTINUE 

DO 20 I=1,MINEQ 
20READ(5,*) IJ,SL(I),SU(I) 

WRITE (6, 927) ARUN 
927 FORMAT ("1"/5X*RUN NO=*8A10/) 

WRITE (6, 718) NVAR,MINEQ 
718 FORMAT (/10X*NO. OF YARIABLES=* I5/10X 
$*NO. OF CONSTRAINTS=*I5) 

WRITE (6, 720) (XBARO(L) ,L=1,NVAR) 
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720 FORMAT ( 10X*VARIABLES=*4 (T30 , 5F10 . 4/) ) 

DO 24 J=1,NVAR 

WRITE(6,22) J, XL(J),XU(J),AVE(J),AMP(J) 

22 P0RMAT(5X*J,XL,XU,AVE,AMP=*I5,4F15.4) 

24 CONTINUE 

DO 26 I=1,MINBQ 
WRITE(6,928) I,SL(I) ,SU(I) 

928 FORMAT(5X*I,SL(I) ,SU(I)=*I5,2F10.2) 

FACT(I)=1. 

XINEQ(I)=1. 

IF(SL(I) .EQ.-999.) FACT(I)=0. 

GNORM (I) = (ABS (SU (I) )+ABS (SL (I) ) )/2.0 

26 CONTINUE 

DO 27 I=1,NVAR 

27 XBAR(I)=XBARO(I) 

READ(5,*) NOJIEL 
WRITE (6, 39) NONEL 

39 FORMAT (10X*NUMBER OF REQUESTED NELMINS=*I5) 

READ(5,*) SCF,RBQMIN,CAYY,STEP1,ILINE 

CALL XINPUT 

IOUT=0 

ILINER=ILINE 

ILINE=0 

JCNT=0 

****************************** 

* START OPTIMIZATION * 

****************************** 

DO 1020 1=1, NONEL 
DO 1000 K=1,NVAR 
1000 STEP(K)=STEP1/I 
CAYY=CAYY*10. 

REQMIN=REQMIN/10 

KCNT=0 

WRITE(6,904) (XBAR(L) ,L=1,NVAR) 

904 FORMAT (/10X*INITAIL XBAR*5(T35,5F12.4/) ) 

XSIGN=1. 

DO 915 IJ=1,NVAR 

IF(ABS(XBAR(IJ)).LT.0.8) GO TO 915 
XO=XBAR(IJ) 

XSIGN=SIGN (XSIGN,XBAR(IJ) ) 

DXB=ABS (XBAR ( I J ) ) -1 . 

IF(DXB.GE.l. .OR.DXB.EQ.O.)XBAR(IJ)=XBAR(IJ)/(ABS(XBAR(IJ) )+.25) 
IF(DXB.GT.O.O.AND.DXB.LT.l.O) XBAR(IJ)=XBAR(IJ)-XSIGN*DXB*2. 
STEP(IJ)=SIGN(STEP(IJ) ,-XBAR(IJ) ) 

WRITE(6,911) IJ,XO,XBAR(IJ) ,STEP(IJ) ,XSIGN 

911 FORMAT(3X" RESET (VAR,XO,XBAR,STEP,XSIGN)="I5,4F12.4) 

915 CONTINUE 

WRITE(6,906) (STEP(L) ,L=1,NVAR) 

906 FORMAT (10X*INITAIL STEPS*5(T35,5F12.4/) ) 

WRITE (6,908) REQMIN,SCF,CAYY 
908 FORMAT(10X*REQMIN,SCF,CAYY*T35,3E12.4) 
************************************^* 11 * 

* START NELDER-MEAD SUBROUTINE * 
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ic’k'kiticic'k'k’kiticit’kici^icit'kic'k'kicic’kicitisicirftic'k'kiiiticicirkit 


ICOUNT=1500 

CALL NELMIN (NVAR , XBAR, XMIN , XSEC , YNEWLO , YSEC , REQ^lIN , STEP , 

$ ICOUNT) 

WRITE (6, 806) 

806 FORMAT (/5X*/ / / NELMIN COMPLETE/ / / */) . 

WRITE(6,810) (XMIN(L),L=1,NVAR) 

810 FORMAT (10X*XMIN=*3(T35,5F15. 8/)) 

WRITE (6, 812) YNEWLO 
812 FORMAT(10X*YNEWLO=*T35,E15.6) 

WRITE (6, 814) ICOUNT, JCNT 

814 FORMAT (10X*ICOUNT=*T35,I5/10X*TOT. FUNCTIOJ CALLS=*T35,I5/) 

xiwi^iv/r 

XJCNT=JCNT 
RAIWT=XIWr/XJCNT 
WRITE (6, 830) XIWT^RAIWT 

830 FORMAT (/10X,22HT0T. WEIGHT ITERATIONS, T35,F12. 2/ 

$ 10X,28HAVE. WT. ITERATIONS PER CALL,T35,F12.3) 

DO 1010 K=1,NVAR 
1010 XBAR(K)=XMIN(K)*1.0 
IOUT=l 
ILINE=0 
IOUT=0 

ILINE=ILINER 

*1rk1rk1rk1t1t1i1fk1fk1rk-k-k1rk1i1t1rk1t1rk*-k1cirk1i1t1tirkirk'kitic1iirlcicirk1rk-k1i1e1t1t1titicicit**1ftc1rk‘k1rk1c1fk 


* THE FOLLOWING STATEMENTS CAN BE USED TO FIND THE GRADIENTS * 

* AT THE OPTIMAL SOLUTION POINT; * 

* * 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 

* 


DX=1.0E-5 

DO 610 JI=1,NVAR 

IDX=0 

603 DO 605 J=1,NVAR 
XBARO(J)=XMIN(J) 

IF(J.EQ.JI) XBARO(J)=XMIN(J)+DX 


605 CONTINUE * 

VAL=EN(XBARO) * 

FX ( JI ) = (VAL-ORIG) /DX * 

IF(FX(JI).LT.1.0E4) GO TO 610 * 

IF(IDX.GT.O) GO TO 610 * 

IDX=1 * 

DX=-DX * 

GO TO 603 * 

610 CONTINUE * 

DO 620 JJI=1,NVAR * 

WRITE (6, 607) JJI,FX.(JJI) * 

607 FORMAT(10X*DERIVATIVE WITH RESPECT TO VARIABLE NO.* * 
$ I5,10X,E15.5) * 

620 CONTINUE * 


•k**1t1t*-k*1rk-k1rk1t1c*1t*1:'k'k1fk**1t1rlrk1rk'k1e1e-k*1e1e1rk1t*'^*irkii1e*1c1t1rk1tif1e*1rk1e1c1rk'k1f1rk'k'k*1t1t 


1020 CONTINUE 

CALL XOUTPUT(4) 

WRITE(6,810) (XMIN(L) ,L=1,NVAR) 
STOP 
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END 

SUBROUTINE NELMIN (N , START , XMIN , XSEC , YNEWLO , 
lYSEC , REQMIN , STEP , ICOUNT) 

REAL START (N) ,STEP(N) ,XMIN(N) , 

IXSEC(N) ,YNEWLO,YSEC, REQMIN, P(20, 21) ,PSTAR(20) , 

2P2STAR(20) ,PBAR(20) ,Y(20) ,DN,Z,YLO,RCOEFF, 
3YSTAR,ECOEFF,Y2STAR,CCOEFF,FN,DABIT,DCHK, 

4COORDl,COORD2 

DATA RCOEFF/1 . 0/ , ECOEFF/2 . 0/, CCOEFF/0 . 5/ 

DATA PSTAR,P2STAR,PBAR /60*0./ 

KCOUNT=ICOUNT 

ICOUNT^O 

IF( REQMIN .LE. 0.) ICOUNT=ICOUNT-l 

IF(N .LE. 0) ICOUNT=ICOUNT-10 

IF(N .GT. 20 ) ICOUNT=ICOUNT-10 

IF(ICOUNT .LT. 0) RETURN 

DABIT=2.04607E-35 

BIGNUM=1.0E38 

KONVGE=5 

XN=FLOAT(N) 

DN=DBLE (XN) 

NN=N+1 

icicit'k'k'kic'k'k’kicicicic'kic'kicicic'kicicif'kiiciic'kicieicicicic'kiciciciticitiicicif'k'k 

* CONSTRUCTION OF INITIAL SIMPLEX * 

1c********************************************* 

DO 1 1=1 ,N 

1 P(I,NN)=START(I) 

Y(NN)=FW (START) 

ICOUNT=ICOUNT+l 
DO 2 J=1,N 
DCHK=START(J) 

START ( J) =DCHK+STEP ( J ) 

DO 3 1=1 ,N 
3 P(I,J)=START(I) 

Y(J)=FN (START) 

ICOUNT=ICOUNT+l 

2 START (J) =DCHK 

************************************************************************ 


* SIMPLEX CONSTRUCTION COMPLETE * 

* * 

* FIND HIGHEST AND LOWEST Y VALUES * 

* * 

* YNEWLO (Y(IHI)) INDICATES THE VERTEX OF THE SIMPLEX TO BE * 

* REPLACED. * 


************************************************************************ 
1000 YLO=Y(l) 

YNEWLO=YLO 

ILO=l 

IHI=1 

DO 5 1=2 ,NN 

IF(Y(I) .GE. YLO) GO TO 4 
YLO=Y(I) 

ILO=I 
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4 IF(Y(I) .LE. YNEWLO) GO TO 5 
YNEWLO=Y(I) 

IHI=I 

5 CONTINUE 

*************************************************** 

* PERFORM CONVERGENCE CHECKS Ol FUNCTION * 

*************************************************** 

DCHK= ( YNEWLOfDABIT) / (YLO+DABIT) -1 . 

IF(ABS(DCHK) .LT. REQMIN) GO TO 900 

KONVGE=KONVGE-l 

IF(KONVGE .NE. 0) GO TO 2020 

KONVGE=5 

************************************1(*1rk1t1t**********************it1c* 

* CHECK CONVERGENCE OF COORDINATES OMLY EVERY 5 SIMPLEXES * 
******************************************************************* 

DO 2015 1=1 ,N 
COORDl=P(I,l) 

COORD2=COORDl 
DO 2010 J=2,NN 

IF(P(I,J) .GE. COORDl) GO TO 2005 
COORDl=P(I,J) 

2005 IF(P(I,J) .LE. COORD2) GO TO 2010 
COORD2=P(I,J) 

2010 CONTINUE 

DCHK= (COORD2+DABIT) / (COORDl+DABIT) -1 . 

IF(ABS(DCHK) .GT. REQMIN) GO TO 2020 
2015 CONTINUE 
GO TO 900 

2020 IF(ICOUNT .GE. KCOUNT) GO TO 900 
************************************************************************ 

* CALCULATE PEAR, THE CENTROID OF THE SIMPLEX VERTICES EXCEPTING * 

* THAT WITH Y VALUE YNEWLO. * 

**************************************************./fit*.),.)c1c**************** 

DO 7 1=1 ,N 
Z=0.0 

DO 6 J=1,NN 

6 Z=Z+P(I,J) 

Z=Z-P(I,IHI) 

7 PBAR(I)=Z/DN 

************************************************** 

* REFLECTIOI THROUGH THE CENTROID * 

************************************************** 

DO 8 1=1 ,N 

8 PSTAR ( I ) = ( 1 . O+RCOEFF) *PBAR ( I ) -RCOEFF*P ( I , IHI ) 

YSTAR=FN (PSTAR) 

ICOUNT=ICOUNT+l 

IF(YSTAR .GE. YLO) GO TO 12 

IFdCOUNT .GE. KCOUNT) GO TO 19 

****************************** ****************ifi,i,ic** 

* SUCCESSFUL REFLECTION, SO EXTENSION * 

**************************************************** 

DO 9 I=1,N 

9 P2STAR ( I ) =ECOEFF*PSTAR ( I ) + ( 1 . O-ECOEFF) *PBAR ( I ) 
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Y2STAR=EN(P2STAR) 

IC0UNT=IC0UNT+1 

************************************************ 

* RETAIN EXTENSION OR C(»ITRACTI(»J * 

************************************************ 

IF(Y2STAR .GE. YSTAR) GO TO 19 

10 DO 11 1=1 ,N 

11 P(I,IHI)=P2STAR(I) 

Y(IHI)=Y2STAR 

GO TO 1000 

***************************** 

* NO EXTENSION * 

* * ifr * 4r 4r 4r 4r ir ★ 4r lAr * * * * * * 

12 L=0 

DO 13 1=1 ,NN 

IF(Y(I) .GT. YSTAR) L=L+1 

13 CONTINUE 

IF(L .GT. 1) GO TO 19 
IF(L .EQ. 0) GO TO 15 

*********************************************************Kit***i( 

* CONTRACTIOM ON THE REFLECTION SIDE OF THE CENTROID * 

********************************************************i,1,ic***ic 

DO 14 1=1 ,N 

14 P(I,IHI)=PSTAR(I) 

Y(IHI)=YSTAR 

*************************************************************** 

* CONTRACTIC^I ON THE Y(IHI) SIDE OF THE CENTROID * 

*************************************************************** 

15 IF(ICOUNT .GE. KCOUNT) GO TO 900 
DO 16 1=1 ,N 

16 P2STAR ( I ) =CCOEFF*P ( I , I HI ) + ( 1 . O-CCOEFF) *PBAR ( I ) 
Y2STAR=FN{P2STAR) 

ICOUNT=ICOUNT+l 

IF(Y2STAR .LT. Y(IHI)) GO TO 10 

*****************************************1c*-*.* 

* CONTRACT THE WHOLE SIMPLEX * 

***********************************1f1,1l1c**1c*** 

DO 18 J=1,NN 
DO 17 1=1 ,N 

P(I,J)=(P(I,J)+P(I,ILO))*0.5 

17 XMIN(I)=P(I,J) 

Y(J)=FN(XMIN) 

18 CONTINUE 
ICOUNT=ICOUNT+NN 

IF(ICOUNT .LT. KCOUNT) GO TO 1000 
GO TO 900 

* RETAIN REFLECTiaJ * 
******************************** 

19 CONTINUE 
DO 20 1=1 ,N 

20 P(I,IHI)=PSTAR(I) 

Y(IHI)=YSTAR 
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GO TO 1000 

*****************■*!****************************************************** 

* SELECT THE TW) BEST FUNCTIOI VALUES (YNEWLO AND YSEC) AND THEIR * 

* COORDS. (XMIN AND XSEC) . * 

************************************************************************ 

900 DO 23 J=1,NN 
DO 22 1=1 ,N 

22 XMIN(I)=P(I,J) 

Y(J)=EN (XMIN) 

23 CONTINUE 
YNEWLO=BIGNUM 
DO 24 j=i,m 

IF(Y(J) .GE. YNEWLO) GO TO 24 
YNEWLO=Y(J) 

IBES-^C 

24 CONTINUE 
Y(IBEST)=BIGNUM 
YSEC=BIGNUM 

DO 25 J=1,NN 

IF(Y(J) .GE. YSEC) GO TO 25 
YSEC=Y(J) 

ISEC<r 

25 CONTINUE 
DO 26 1=1 ,N 
XMIN(I)=P(I,IBEST) 

XSEC(I)=P(I,ISEC) 

26 CONTINUE 
RETURN 
END 

FUNCTION FN(XBAR) 

COMMON A^IOL/CAYY,MC,ILINE,IOUT 

COMM(»l /FCOUNT/ICNT,JCNT,KCNT,NDAV,LCNT,IWr,NOIT 

REAL XBAR(15) 

CALL SETUP (XBAR,IOUT, OBJ) 

FN=OBJ 

IF(ILINE.GT.O) WRITE(6,50) JCNT,FN,MC,IWT,NOIT 
50 FORMAT (lOX /'CALL NUMBEF^" ,I5,5X"OBJ="E15.8,10X 
$ ,"NO. VIOLATIONS="I5,5X 
$ /'IWr=",I5,5X/'NOn^"I5) 

RETURN 

END 

SUBROUTINE SETUP (GAINS, IPR, OBJ) 

INTEGER OUTPUT 

REAL GAINS (15) ,GBAR(59) 

COMMOJ /STRAIN/CON (59) 

COMM(»l AARIAB/AMP(15) ,AVE(15) 

COMMON AlOL/CAYY,MC,ILINE,IOUT 
COMMOl /CONSWSU(59) ,SL(59) ,XINEQ(59) 

COMMON /DEVAR/DESIGN(15) ,ITERM(10) ,CST(10) 

COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /PERFAJNAUG,SCF,NVAR,MINEQ 
COMMON /AVOID/FACT(59) ,GNORM(59) 

DO 22 I=1,NVAR 
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DESIGN (I ) =AMP ( I ) *SIN (1 . 5707963*GAINS ( I ) ) +AVE (I ) 

22 CONTINUE 

************************************************************************ 


* TO DEFINE ADDITICWAL DESIGN VARIABLES, INSERT DEFINITIOI CARDS * 

* HERE. SAMPLES ARE GIVEN BELOW: * 

* * 


* 

* 

* 

* 

* 


GX(3)=DESIGN(8) 

W(1)=DESIGN(9) 

W(4)=DESIGN(10) 

W(3)=DESIGN(11) 

GX(5)=DESIGN(12) 


—DESIGN (8) IS MACH NUMBER 
—DESIGN (9) IS WING SWEEP 
—DESIGN (10) IS WING THICKNESS RATIO 
—DESIGN (11) IS WING TAPER RATIO 
—DESIGN (12) IS FUSELAGE DIAMETER 


* 

* 

* 

* 

* 


************************************************************************ 


* ALSO, VARIOUS DESIGN CONSTRAINTS CAN BE ADDED. FOR EXAMPLE, * 

* ADDING THE STATEMENT; * 

* * 


* GX(32)=W(1) * 

* RESTRICTS THE HORIZONTAL TAIL SWEEP ANGLE TO BE * 

* EQUAL TO THE WING SWEEP ANGLE. * 

************************************************************************ 


IF(PX(5) .GT.O) DESIGN(6)=GX(21)*ITERM(4) 

IF(IPR.GT.O) WRITE(6,441) DESIGN 

441 FORMAT (//15X*SET-UP*/20X*DESIGN=*5(T40,5F15. 4/) ) 
OUTPUT=0 

CALL DOCOST(COST,0) 

CALL CNSTRN(O) 

UNAUG=COST*SCF 

IF(IPR.GT.O) WRITE(6,442) UNAUG 

442 FORMAT (20X*CALL TO DOC/CNSTRN COMPLETE ^DOC= $*E15.5) 

********************************************** 


* THIS SECTiai CALCULATES GBAR ARRAY * 

* AND THE PENALTY FOR VIOLATIONS * 

********************************************** 

MC=0 

PENT=0. 

DO 150 I=1,MINEQ 
T=CON(I) 

XINEQ(I)=0. 

GBAR(I)=AMAX1 (T-SU(I) ,SL(I)-T) 

GBAR ( I ) =GBAR ( I ) *FACT ( I ) /GNORM ( I ) 
IF(GBAR(I) .LE.0.0) GO TO 150 
XINEQ(I)=1. 

MC=MC+1 

PENT=PENT+CAYY*GBAR ( I ) *GBAR ( I ) 

150 CONTINUE 

IF(IPR.GT.O) WRITE(6,160) PENT 
160 FORMAT (20X*PENALTY TERM=*E15.5) 

'kis'k'kicirkirk-k-k'kicieic-k'k'kicieicieicicit'k'k'kicicic’k'kie'kic'kie'kicie'k'k 

* AUGMENTED FUNCTION IS CREATED * 

******************************************* 


OBJ=UNAUG+PENT 

RETURN 

END 

SUBROUTINE XOUTPUT(IPRNT) 
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COMMON /CONSTR/SU(59) ,SL(59) ,XINEQ(59) 

COMMON /DEVAR/DESIOa (15) ,ITERM(10) ,CST(10) 

COMMON /DRAG/CDS (6), CDSAP (6) 

COMMON /GEOl/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /GRAVn7/0G(6) 

COMMON /LABELP/ARUN(8) 

COMMON /STAB/DERIVCR(15) ,DERIVAP(15) ,STOR(20) 

COMMON /WrSVE/WTS(20) 

************************************************************************ 

* PRINT OUT PARTS OF THE FUNCTION EVALUATION * 

* * 

* 0=NONE,1=DOC,2=CNSTRN,3=DOC & CNSTRN,4=DOC & CNSTPN & DUMP * 

* 5=DOC(2) & CNSmi(2) ,6=COMMON DUr-IP 

************************************************************************ 

IDUMP=0 

IF(IPRNT.LT.l) GO TO 999 
IF(IPRNT.EQ.l) GO TO 10 
IF(IPRNT.EQ,2) GO TO 20 
IF(IPRNT.EQ.3) GO TO 30 
IF(IPRNT.EQ.4) GO TO 40 
IF(IPRNT.EQ.5) GO TO 50 
IF(IPRNT.EQ.6) GO TO 70 
10 IDOC=l 
ICRN=0 
GO TO 80 
20 ICRN=1 
IDOC=0 
GO TO 80 
30 IDOC=l 
ICRN=1 
GO TO 80 
40 IDOC=l 
ICRN=1 
IDUMP=1 
GO TO 80 
50 IDX=2 
ICRN=2 
IDUMP=0 

80 CALL DOCOST(TERM,IDOC) 

CALL CNSTRN(ICRN) 

IF (IDUMP.lt. 1) GO TO 999 
70 WRITE (6, 1010) 

1010 FORf4AT(*l*30X*COMT4ON DUMP*///) 

WRITE (6, 1012) DESIGN 

1012 F0RTV^T(/15X*DESIGN=*3(T30,5F15.5/) ) 

WRITE (6, 1014) ITERT4 

1014 FORMAT (/15X*ITERM=*3(T30, 5115/)) 

WRITE (6, 1016) PX 

1016 FORMAT(/15X"PX=",3(T30,5F15.4/) ) 

WRITE (6, 1015) CDS 

1015 FORMAT (15X*CDS=*T30,6F15. 9/) 

WRITE (6, 1020) CDSAP 

1020 FORMAT(/15X*CDSAP=*T30,6F15.9/) 

■ ' 2T ' 



WRITE (6, 1025) W 

1025 FORMAT(/15X*W=*5(T30,5F15.9/)) 

\'JRITE(6,1030) HX 

1030 FORMAT(/15X*HX=*5(T30,5F15.9/) ) 

VJRITE(6,1035) GX 

1035 FORMAT (/15X*GX=*8 (T30,5F15.9/) ) 

WRITE (6, 1040) CG 
1040 FORMAT(/15X*OG=*T30,6F15.5) 

WRITE (6, 1045) DERIVCR 
1045 FORMAT (/15X*DERIVCR=*5(T30,5F15. 9/) ) 

WRITE (6,1050) DERIVAP 
1050 FORT4AT (/15X*DERIVAP=*5 (T30 ,5F15.9/) ) 
lfJRITE(6,1055) STOR 

1055 FORMAT (/15X*STOR=*5(T30,5F15. 9/) ) 

VJRITE(6,1060) WTS 

1060 F0Rf^T(/15X*WrS=*5 (T30,5F15.3/) ) 
l'iJRITE(6,1070) CST 

1070 FORMAT (/15X"CST="5 (T30,5F15.3/) ) 

WRITE (4, 1200) ARUN 
1200 FORMAT (8A10) 

************************************************************** 

* WRITE AN OUTPUT TAPE FOR THE PLOTTING PREPROCESSOR * 
************************************************************** 

WRITE (4,*) DESIGN 
WRITE (4,*) ITERM 
WRITE(4,*) W 
WRITE{4,*) HX 
VJRITE(4,*) GX 
999 RETURN 
END 

SUBROUTINE XINPUT 

COrWN /DEVARA)ESIGN(15) , ITERM (10) ,CST(10) 

COM^iON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

C0M40N /GRAVITY/CG(6) 

COMMON AiTSVEA)TS(20) 

********************************** 

* INPUT DESIGN CONSTANTS * 
********************************** 

DATA GX/35*0./ 

READ (5,*) WTS(l) 

READ (5,*) (PX(J) ,J=1,8) 

READ(5,*) (ITERM(I) ,1=1,10) 

READ(5,*) WTS(16) ,OG(5) 

READ(5,*) W(l) ,W(2) ,W(3) ,W(4) ,W(5) 

READ (5,*) W(6) ,W(7) ,W(8) ,W(14) ,W(16) 

READ(5,*) W(17) ,W(18) ,W(19) ,W(20) 

READ(5,*) HX(1) ,HX(2) ,HX(3) ,HX(10) 

READ(5,*) HX(16),HX(17),HX(18),HX(19),HX(20) 

READ(5,*) GX(3) ,GX(4) ,GX(5) ,GX(6) ,GX(7) 

READ(5,*) GX(8),GX(11),GX(12),GX(17),GX(18) 

READ (5,*) GX(19) ,GX(20) ,GX(21) ,GX(22) 

READ(5,*) GX(23) ,GX(24) ,GX(25) ,GX(26) ,GX(27) 

READ(5,*) GX(32) ,GX(34) ,GX(16) 
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************************** 


* WRITE OUT DATA * 

************************** 

WRITE (6, 100) 

100 FORMAT (*1*///30X*/ / / FUNCTION INPUT/ / / *//) 
VJRITE(6,110) WTS(l) 

110 FORMAT(10X*WTS(1)*T40,F15.4) 

WRITE (6, 120) ITERM 
120 FORMAT (10X*ITERM*T40, 1016) 

WRITE (6, 125) PX 

125 FORMAT(10X/'PX=",3(T40,5F15.4/)) 

WRITE(6,130) WTS(16) 

130 FORMAT(10X*WTS(16)*T40,F15.4) 

VJRITE(6,140) W(l) ,W(2) ,W(3) ,W(4) ,W(5) 

140 FORMAT(10X*W. . .1-5=*T40,5F15.4) 

l*JRITE(6,150) W(6) ,W(7) ,W(8) ,W(14) ,W(16) 

150 FORMAT(10X*W...6-8,14,16=*T40,5F15.4) 

WRITE(6,160) W(17) ,W(18) ,W(19) ,W(20) 

160 FORMAT (10X*W. ..17-20=*T40,4F15. 4) 

WRITE(6,170) HX(1) ,HX(2) ,HX(3) ,HX(10) 

170 FORMAT(10X*HX...1-3,10=*T40,4F15.4) 

WRITE(6,175) HX(16) ,H}((17) ,HX(18) ,HX(19) ,HX(20) 
175 FORr4AT(10X*HX. . .16-20-*T40,5F15.4) 

W’RITE(6,180) GX(3) ,GX{4) ,GX{5) ,GX(6) ,GX(7) 

180 FORMAT(10X*GX. . .3-7=*T40,5F15.4) 

VJRITE{6,190) GX(8) ,GX(11) ,GX(12) ,GX(17) ,GX(18) 
190 FORMAT ( 10X*GX ...8,11-12,17, 18=*T40 , 5F15 . 4 ) 
WRITE(6,200) GX(19) ,GX(20) ,GX(21) ,GX(22) 

200 FORMAT (10X*GX. . .19-22=*T40,4F15.4) 

WRITE(6,210) GX(23) ,GX(24) ,GX(25) ,GX(26) ,GX(27) 
210 FORMAT (10X*GX. . .23-27=*T40,5F15.4) 

VJRITE(6,220) GX(32) ,GX (34) 

220 FORMAT (10X*GX. . .32,34=*T40,2F15.4) 

RETURN 

END 

SUBROUTINE DOCOST(UNAUG, OUTPUT) 

INTEGER OUTPUT 

REAL INSUR,PER(12) ,XIOC(12) ,YIOC(12) ,YCOST(9) 
COr^MON /DEVAR/DESIGN(15) ,ITER1(10) ,CST(10) 

COMMON /FCOUNT/ICNT, JCNT,KCNT,NDAV,LCNT, IWT,NOIT 
COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /WTSVE/OTS(20) 
***************************************** 

* INITIALIZE GEOMETRY AND C.G. * 

***************************************** 

CALL GEO 
CALL OGCAL 

********************************************** 

* INITIALIZE AND INCREMENT COUNTERS * 

********************************************** 

DATA JCNT/0/ 

jCNT=JCNT+l 

KCNT=KCNTH 




IF(OUTPUT.GE.l) GO TO 10 

*****<t**<ri!r*************************************** 

* COSTS PER BLOCK HOUR OF DESIGN FLIGHT * 
************************************************* 

20 IOUT=OUTPUT 

IF(OUTPUT.GT.l) IOUT=0 
CALL WEIGHT (lOUT) 

CALL AIRCOST (PRICE, lOUT) 

YRrmT=1.07** (ITERM(6)-1976) 

PRICE=YRMULT*PRICE 

FL=V/TS(18) 

BLKHR=WTS(19) 

*********************************** 

* DIRECT OPERATING COSTS * 

*********************************** 

DEPRE=0.88*PRICE/(14.0*GX(26) ) 

SUPPORT=0.12*PRICE/(14.0*GX{26) ) 

SPARES=0 .06*PRICE/ (14 .0*GX (26) ) 

DELAY=YRMULT*8 . 40 
INSUR=0 . 01*PRICE/ (GX (26) ) 

FCOST=WTS(20) *WTS(16)/6.4 

FCOST=FCOST/BLKHR 

CALL MAINCST (XMCOST , TOUT) 

^0=^8 (5) *0.453592 
FEELAND=YRMULT*1 . 54*WTL/1000 . 0 
FEELAND=FEELAND/BLKHR 

ATT=YRMULT*GX(19) * (0.691*FL+0.00175*FL*FL) 

ATT=ATT/BLKHR 

CRE?//=YRMULT*174*FL+43.5+(0.452*FL+.11299)*(WTS(l)*. 453592/1000.) 

CR0W=CREW/BLKHR 

SERVICE=YRMULT*63 . 0 

CONTROL=YRMULT*82 . 58/BLKHR 

DOC=DEPRE+SUPPORT+SPARES+DELAY+INSUR+FCOST+FEELAND+SERVICE 

$+ATT+CREW+XMCOST+CONTROL 

* PERCENT OF TOTAL * 

**************************** 

PER(l)=DEPRE/DOC 
PER (2) =SUPPORT/DOC 
PER ( 3 ) =SPARES/DOC 
PER(4)=DELAY/DOC 
PER(5)=INSUR/DOC 
PER(6)=FC0ST/D0C 
PER(7)=XT4COST/DOC 
PER (8 ) =FEELAND/DOC 
PER(9)=CREW/DOC 
PER(10)=ATT/DOC 
PER (11) =SERVICE/DOC 
PER (12) =CONTROL/DOC 
DO 30 J=l,12 
30 PER(J)=PER(J)*100. 

TOT=100. 

******************************************************************** 
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* INDIRECT OPERATING COSTS PER BLOCK HOUR OF DESIGN FLIGHT * 
******************************************************************** 

DATA yiOC/lOHMAIN BURDN, 9HFOOD COST,5HMOVIE,8HPASS INS, 

$ 9HMISC PASS,9HADVERTISE,10HCOMMISSION,5HRESER, 

$ 9HPASS HDLG,8HBA3 HDLG,10HCARGO HDLG,9HSERVICING/ 
XIOC(1)=XMCOST*1.05 
IFIRS=. 15*GX (19) *GX (24) 

IECON=GX (19) *GX (24)-IFIRS 
XIOC (2) = (IFIRS*2 . 42+IECON*! .05) 

XI0C(3)=196./BLKHR 
RPM=GX(19) *GX(24) *GX(4)/1000. 

XIOC (4 ) =0 . 52*RPM/BLKHR 
XIOC (5) =GX (19) * . 18/BLKHR 

REVYR=GX (25) *GX (19) *GX (24) *GX (26) *GX (4) /BLKHR 
REVHR=REVYR/GX (26) 

XIOC(6)=.023*REVHR 
XIOC (7 ) =2 . 35*RPM/BLKHR 
PASSPHR=GX (19) *GX (24) /BLKHR 
XIOC (8 ) =4 . 40*PASSPHR 
XIOC(9)=2.87*PASSPHR 
XIOC (10) =1 . 31*PASSPHR 
TONCAR=GX (22) /2000. 

XIOC (11) =131 . 08*TONCAR/BLKHR 

XIOC (12) = (0 . 03*9 . 5+0. 0025) *GX (19) /BLKHR 

TOTIOC=0. 

XIOC ( 1 ) =XIOC ( 1 ) /YRMULT 
XIOC (6) =XIOC (6) /YRMULT 
DO 200 1=1,12 
XIOC ( I ) =XIOC ( I ) *YRMULT 
TOTIOC=TOTIOC+XIOC (I) 

200 CONTINUE 

******:5r********************:*f ********** ******** 

* RETURN ON INVESTMENT CALCULATICNS * 

********************************************** 

XINVEST=0.9*PRICE 
TAXRT=0.48 
COSTHR=DOC+TOTIOC 
PROFIT= (REVHR-COSTHR) *GX (26) 

ROI= (l.-TAXRT) *PROFIT/XINVEST 

FARROI=(.26*PRICE+COSTHR*GX(26) )* (BLKHR/ (GX (19) *GX(24) *GX(26) 

$ *GX(4))) 

************************************* 

* ASSIGN PERFORMANCE INDEX * 

************************************* 

CST(l)=DOC 

CST(2)=DOC*BLKHR 

CST(3)=ROI 

CST(4)=FARROI 

CST(7)=WTS(1) 

INCPH=COSTHR+. 15*XINVEST/ (GX (26) * (1 .-TAXRT) ) 

INCPF=INCPH*BLKHR 

CST(8)=INCPF 

CST(9)=PRICE 
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UNAUG=CST(ITERM(5) ) 

GO TO 40 

*************************** 

* OUTPUT SECTION * ' 

*************************** 

10 WRITE (6, 42) 

l'JRITE(6,44) (DESIGN (JK),JK=1, 12) 

IF(OUTPUT.GT.l) GO TO 20 
WRITE(6,705) 

705 FORMAT (//5X*INPUT CONSTANTS*/) 
l‘JRITE(6,710) (W(J),J=1,8) 

710 FORMAT(10X*WING...SWEEP,INCIDENCE,TAPER RATIO;*T55,3F12.4/ 

$ 10X*THICKNESS,'IWIST,E1,E2,DESIGN CL;*T55,5F12.4) 

1€IITE(6,720) W(14) ,W(19) ,W(17) ,W(18) ,W(16) ,GX(17) ,GX(18) 

720 FORMAT(10X*CM(CR,APP)*T55,2F12.4/ 

$ 10X*DELTA CM(10,25 DEGREES FLAP) :*T55,2F12.4/ 

$ 10X*ANGLE OF ZERO LIFT(0,10,45 DEGREES FLAP) :*T55,3F12.4) 
VJRITE(6,730) W(20) ,GX(11) ,GX(5) ,GX(8) ,GX(12) 

730 FORflAT(10X*DELTA CD (10-45 DEGREES FLAP) :*T55,F12.4/ 

$ 10X*TURBULENCE LENGTH/ROOT 3, FUSE. DIA. :*T55,2F12.4/ 
$10X*CL-MAX (TO) ,CL-MAX(L) :*T55,2F12.4) 

WRITE (6, 740) GX(3) ,GX(4) ,GX(19) ,GX(22) ,GX(23) ,GX(27) 

740 FORf>lAT(10X*MISSION,.. MACH NO., RANGE, NO. PASS;*T55,3F12.4/ 

$ 10X*CARGO WEIGHT:*T55,F12.4/ 

$ 10X*DELTA CG, WTL(MAX)/WrO:*T55,3F12.4) 

WRITE(6,750) GX(6) ,GX(7) ,GX(20) ,GX(21) 

750 FORMAT (10X*ENGINE. . .L,W,WT,TREF:*T55,4F12.4) 

WRITE(6,760) HX(1) ,HX(2) ,HX(3) ,HX(10) ,GX(24) ,GX(25) ,GX(26) 
760 FOPvT4AT(10X*TAIL... TAPER RATIO,THICKNESS,ELE EFF:*T55,3F12.4/ 
$ 10X*ELEVATOR TIME CONSTANT ;*T55,F12. 4/ 

$ 10X*ECONOMICS...LOAD FACT,$/SEAT MI,BLK HR/YR;*T55,3F12.4) 
l*/RITE(6,770) W(9) ,W(10) ,HX(4) ,HX(5) ,HX(6) ,GX(32) 

770 FORMAT (/5X*SOME GEOMETRY CALCULATIONS*/ 

$ 10X*WING. . .SPAN,CMAC:*T55,2F12.4/ 

$ 10X*TAIL. . .SPAN,CMAC,VBAR,SVJEEP*T55,4F12.4) 

1€IITE(6,780) HX(16) ,HX(17) ,HX(18) ,HX(19) ,HX(20) ,GX(35) 

780 FORr«lAT(10X*VERT. TAIL. . .VBAR,TAPER,AR,SWEEP,SR/SV,SV:* 

$ T55,6F12.4) 

GO TO 20 

40 IF (OUTPUT. LT.l) GO TO 100 
IF(OUTPUT.GT.l) GO TO 549 

42 FORMAT (*1*//30X*AIRCRAFT SIZING PROGRAM*//) 

44 FORMAT (5X*DESIGN VARIABLES*/10X*WING AREA (FTXX2)=* 
$T40,F15.4/10X*WING ASPECT RATIO=*T40,F15.4/10X 
$*FUSELAGE LENGTH (FT)=*T40,F15.4/10X 
$*HOR. TAIL AREA (FTXX2)=*T40,F15.4/10X 
$*HOR. TAIL ASPECT RATIO=*T40,F15.4/10X 

$*TOTAL THRUST (LBS)=*T40,F15.4/10X"AFT MOST CG="T40,F15.4/ 

$ 10X"CRUISE MACH NO.="T40,F15.4/10X"SWEEP="T40,F15.4/10X 
$ "WING T/C="T40,F15.4/10X"WING TAPER RATIO="T40,F15.4/ 

$ 10X,"FUSE. DIA=",T40,F15.4) 

WRITE (6, 52) 

52 FORMAT (*1*//30X*DIRECT OPERATING COSTS--DOLLARS/FLT. HOUR*) 
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WRITE(6,54) DEPRE,PER(1) , SUPPORT, PER (2) , SPARES, PER (3) ,DELAY, 

$ PER(4) ,INSUR,PER(5) ,FCOST,PER(6) ,XMCOST,PER(7) ,FEELAND,PER(8) , 

$ CREW,PER(9) ,ATT,PER(10) ,SERVICE,PER(11) ,CONTROL,PER(12) 

54 FORMAT (//10X*DEPREC*T40 , 2F10 . 2/10X*SUPPORT*T40 , 2F10 . 2/10X*SPARES* 
$T40, 

$2F10 . 2/10X*DELAY*T40 ,2F10 . 2/10X*INSURANCE*T40 , 2F10 . 2/10X*FUEL*T40 , 
$2F10 . 2/10X*MAINTENANCE*T40 , 2F10 . 2/10X*LANDING FEE*T40 , 2F10 , 2/lOX 
$*CREW*T40 , 2F10 . 2/10X*ATTENDANTS*T40 , 2F10 . 2/10X*FUEL SERVICE* 

$ T40,2F10.2/ 

$1 OX*CONTROL*T40 , 2F10 . 2) 

WRITE (6, 56) DOC,TOT 

56 FORMAT (/3X*TOTAL DIRECT OPERATING COSTS*T40*$*F9.2,F10.2) 

WRITE (6, 150) 

150 FORMAT (///30X*INDIRECT OPERATING COSTS— DOLLARS/FLT. HOUR*///) 

DO 300 1=1,12 

PER (I ) =100 . 0*XIOC ( I ) AOTIOC 
WRITE(6,152) YIOC(I),XIOC(I),PER(I) 

152 FORT4AT(10X,A10,T40,2F10.2) 

300 CONTINUE 

WRITE (6, 154) TOTIOC,TOT 

154 FORMAT (/5X*TOTAL INDIRECT OPERATING COSTS*T40,2F10.2/) 

549 WRITE (6, 550) REVHR,COSTHR,ROI 

550 FORMAT (*l*///30X*PERF0Rf4ANCE FUNCTION SUMMARY*/// 

$ 10X*REVENUE PER BLOa< HOUR*T50,F12.2/ 

$ 10X*TOTAL COST PER BLOCK HOUR*T50,F12.2/ 

$ 10X*RETURN ON INVESTMENT*T50,F12.4///) 

DATA YCOST/6HDOC/HR,7HDOC/FLT, 3HROI , 4HFARE , lOHSEAT-MI/GA, 

$ 8HL/D (MAX) , 5HMTOGW,4HFARE , 5HPRICE/ 

DO 570 1=1,9 

Vfl^ITE(6,560) I,YCOST(I) ,CST(I) 

560 FORMAT(10X,I5,2X,A10,F12.3) 

570 CONTINUE 
100 RETURN 
END 

SUBROUTINE CGCAL 

COf4MON /DEVAR/DESIGN(15) ,ITERM(10) ,CST(10) 

COr^MON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /STAB/DERIVCR(15) ,DSRIVAP(15) ,STOR(20) 

COmON /GRAVITY/ OG(6) 

******************************* 

* ' ASSIGN CG POSITIONS * 

* (1) AFT-CRUISE * 

* (2) AFT-APPROACH * 

* (3) • EWD-CRUISE * 

* (4) FWD-APPROACH * * 

******************************* 

OG(l)=DESIGN(7) 

CG(2)=CG(1) 

OG(3)=OG(1)-GX(23)/W(10) 

CG(4)=CG(2)-CX(23)/W(10) 

DATA CG(6)/.18/ 

C 

C ALLOW FREE GEAR LOCATION 
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c 

IF(ITERM(7) .GT.O.) OG(5)=CG(l)-CG(6) 

RETURN 

END 

SUBROUTINE GEO 

COMMON /DEVAR/DESIGN(15) ,ITERM(10) ,CST(10) 

COMMON /GEOM/W(20) ,HX(20) ,GX{35) ,PX(15) 
********************************************** 

* WING CONSTANTS OR COMMON VARIABLES * 
********************************************** 

W{9)= (DESIGN (1) *DESIGN (2) ) **0 .5 
W(10)=4.*DESIGN(1)*(1.+W(3)+W(3)*W{3))/ 

$ (3.*W(9) * (l.+W(3) ) **2.0) 

W(11)=TAN(W(1)/GX(2)) 

W(12) =.49640. 45/DESIGN (2) 

W(13)=W(12)-. 9/DESIGN (2) 

W(15)=COS(W(l)/GX(2)) 

*******^************************************************* 

* HORIZONTAL TAIL CONSTANTS OR COMMON VARIABLES * 
********************************************************* 

HX (4)= (DESIGN (4) *DESIGN (5) ) **0.5 

HX (5) =4 . *DESIGN (4) * (1 .+HX (1 ) +HX (1 ) *HX (1 ) ) / 

$ (3 .*HX (4) * (1 .+HX (1) ) **2 .0) 

XTXLT=0.45 

IF(ITERM(4) .EQ.3) XTXLT=0.43 
XLT=XTXLT*DESIGN (3) 

GX(33)=XLT 

HX (6)=DESIGN (4) *XLT/ (DESIGN (1) *W(10) ) 

HX (7) =TAN (GX (32) /GX (2) ) 

T=(1.-HX(1) )/(4.* (l.+HX(l) ) ) 

HX (8 ) =HX (7 ) 4-4 . *T/DESIGN ( 5 ) 

HX(9)=HX(8)-2.*(1.-HX(1))/(DESIGN(5)*(14HX(1))) 

*********************************************** 

* MISC. CONSTANTS OR COMMON VARIABLES * 
*********************************************** 

DATA GX(1) ,GX(2)/3. 14159265, 57. 295779/ 

DATA HX (15) ,GX (13) ,GX (31) /O . ,0 . ,0 ./ 

DATA GX(15)/32.174/ 

GX (35) =HX (16) *W(9) *DESIGN (1) / (XLT*0 .95) 

RETURN 

END 

SUBROUTINE WEIGHT (OUTPUT) 

REAL LT,F1(2) ,F2(2) ,ANS(4) ,DWTA(41) 

COMMON /DEBUG/IDEBIK3,IDEBUG2 

COMMON /FCOUNT/ICNT, JCNT,KCNT,NDAV,LCNT, IWT,NOIT 
COMMON /DEVAR/DESIGN (15) ,ITERM(10) ,CST(10) 

COMMON /DRAG/CDS (6) ,CDSAP (6) 

COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /WTSVE/WrS(20) 

INTEGER OUTPUT 
DATA WFOWTO/.26/ 

DO 20 1=1,41 
20 DWTA(I)=0. 
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************************************************* 

* FUDGE IS A FACTOR FOR WEIGHT OVERRUNS * 
*****************< 1 !******************************* 

FUDGE=1.05 

IF(IDEBUG.GT.O) WRITE(6,337) 

337 FORMAT (*1*/20X*WEIGHT ITERATION LOOPS*) 
************************************** 

* SOME GEOMETRY DEFINITIONS * 
************************************** 

IACT=ITERM(1) 

IGGX=ITERM{2) 

AR=DESIGN(2) 

LT=DESIGN (1) *W(10) *HX (6) /DESIGN (4) 

SHT=DESIGN (4) 

TRV=HX(17) 

SVT=HX (16) *W(9) *DESIGN (1) AT 
TENG=DESIGN (6) /ITERM(4) 

NOIT=0 

DELWTO=20. 

DIV=1. 

WTO=WTS(l) 

WTINIT=WTO 

WTFUEL=WTO*WFOWTO 

FUEL=WTFUEL/6.4 

WTENG=FUDGE* (TENG/GX(21) ) *GX(20) *ITERM(4) 

WTS (7) =WTS (1) / (GX (15) *0 .000889*W(10) ^DESIGN (1) ) 

WTS (9 ) = (WTS (1) -WTFUEL) / (GX (15) *0 . 002378*W(10) *DESIGN (1) ) 
******************************************************************** 

* THIS SECTION CALCULATES l*JTS. INDEPENDENT OF WTO AND FUEL * 
******************************************************************** 

Fl(l)=l. 

Fl(2)=1.15 

WECTL=F1 (IACT+1) *88.46* ( (190. +W (9) ) *4.0E-2) **0.294 
WTSRT=49 . 19* (4 .0*WTENG*1 .OE-3) **0 . 541 
PASS=GX(19)*170. 

WTFURN=39 . 51*GX (19) 

WTFOOD=214.5 

WT02=300.7 

WIWIN=501.55 

WTBGH=144.72 

WTAC=3647. 

WTTR=1500. 

F2(l)=0. 

F2(2)=250. 

ACTCON=F2(IACTH) 

DATA CREW/1700./ 

BAGGAGE=GX (19) *35 .0 
CARGO=GX(22) 

************************************************************ 

* THIS SECTION COMPUTES WEIGHTS DEPENDENT UPON WTO * 

* BEGIN WEIGHT ITERATION THIS SECTION * 

************************************************************ 

50 WIWING=2 . 0*FUDGE*0 . 00428* (DESIGN (1 ) ) **0 . 48*AR* (GX (3) -0 . 05) **0 . 43 
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$*{WrO*1.5*(3.0-PX(8)))**0.04*W(3)**0.14/((100.0*W(4))**0.76 
$* (COS (W(l) /GX (2) ) ) *1 . 54) 

FUDHT=1. 


IF(DESIGN(5) .GT.5.0) FUDHI^l.+(DESIGM(5)-5.0)/10. 

*****************************A*fE*A**********;V***vV*********************** 


* THIS IS A LINEAR CORRRECTION FACTOR FOR TAIL ASPECT RATIOS 5; * 

* INTENDED AS A PENALTY TERT-l TO TRY AND REFLECT THE RANGE OF * 

* TAIL WEIGHT EQUATION VMJDITf * 

************************************************************************ 


TTEMP=HX (5 ) *1 . 5*HX (2 ) 

TET^ (WTO*4 .5)**0.813* (S!IT**n , 584 ) * (HX (4 ) /TTEMP) **0.033 
$*(W(10)AT)**0.28 

VmiTl=2.0* (0.0034*TMP**0c915) *FUDGE 

WTHT2=2.0*. 00563* (lVTO**Oo6) * (SHT**.469) * (DESIGN (5)/.75) **.539 
$ *((1.+HX(1))/HX(2)) **0.692 

V/THT3=1 .0* (4 . 566*1 . OE-4* (DESIGN (4) ** ,48) *DSSIGN (5) * 

$ (l"/TO*4. 5/1.) **0.84) 

************************************************************************ 


* lilTHTl — FROM NICOLAI (REF. 4) * 

* WTHT2 — FROM VDEP (REF. 6) * 

* WTHT3 — ANALOGOUS TO WING WEIGHT EQUATION * 

* * 

* THE THREE HORIZONTAL TAIL I'JEIGHT EQUATIONS ARE THEN AVERAGED * 

************************************************************************ 


WTHT= ( (WrHTl4WTHT2+WTHT3) /3 . 0) *FUDHT 
TEMPV=1 .02* (4 . 5*WTO) ** . 363*SVT**1 ,089* (GX(3) *0.8) ** .601 
$ *LT** (-.726) * (l.+HX(20) ) **.217*HX(18) **.337* 
$(l.+TRV)**.36*(COS(HX(19)/GX(2) ))**(-. 484) 

WTVT= (0 . 19*TEMPV**1 . 014 ) *FUDGE 

TO=10.43* ( .000364* ( (GX (3)+0 .04) *971.15) **2.0) **0.283 
lAfTFUSEl=( (TQ* (DESIGN(3)/GX(5) ) **0.71)* (WTO*1.0E-3) **0.95) *FUDGE 
VJTFUSE2=0 .0796*2 . 1861* (VJTO**0 . 33) * (DESIGN (3) **0.76) 

$ *(GX(5)4GX(5))**1.2 

************************************************************* 

* WTFUSEl — FROM NICOLAI (REF. 4) * 

* WTFUSE2 — FROM VDEP (REF. 6) * 

* * 

* THE TWO FUSELAGE WEIGHT EQUATIONS ARE THEN AVERAGED * 

************************************************************* 


WTFUSE= (WTFUSE1+WTFUSE2) /2 . 

VJTLG=2 . 0* (62 . 21* (WTO*l . OE-3 ) **0 . 84) *FUDGE 
WTCTL= (56.01* (WTO*245 , 6*1 . OE-5) **0 . 576) *FUDGE 

WTINST=2 . 0* ( 1 5 .+ . 0 32*WTO*l . E-3 ) +4 . * (4 , 8+. 006*WTO*1 . OE-3 ) + . 15*l‘JTO*l 
$.0E-3 

wmisc= (0 .771*wr0*l .0E-3*1 . 1) ’-FUDGE 
IF(NOIT.LT.20) GO TO 94 
IF(DELWTO.LT.20.0) GO TO 96 
**************^******************* 

* CALCULATE FUEL WEIGHT * 

94 CALL FUELCAL(WrO,V‘/TFUEL,0) 

FUEL=WTFUEL/6.4 
GO TO 98 
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96 WrFUEL=M’0*WFCWT0 
FUEL=WTFUEL/6.4 

* THIS SECTION CALCULATES llTTS. DEPENDENT UPON FUEL WT. * 

******************************5V********************************* 

98 WTFl=41.6*(FUEL*1.0E-2) **0.818 
lfTF2=7 .91* (FUEL*1 . OE-2) **0 .854 
WTF3=7. 38* (FUEL*1.0E-2) **0.458 
V«’F4=IOGX*28.38* (FUEL*1.0E-2) **0.442 
WTFSYS=WTFl-WTF2+VJTF3-!-li7rF4 

V'JTELEC=1162.66*( (lVTFSYS+931.3)*1.0E-3) **0.506 
************************************?{*************************** 

* THIS SECTION CALCULATES THE EMPTY AND FIXED WEIGHTS * 

***************************************ft****iVA****************** 

Tl=W^WING4WTHT+WTVT+V^’FUSE+l•JTLG+V‘JTCTL+l^iTINST+WTMISC+WECTL4W^SRT 

$-WFURN-H»mD2+WIWIN+l'JTEGH-n^AC+l'JTTR+ACTC^ 

WTEMPTY= (T14WTENG+OTELEC+VJTFSYS) *FUDGE 
l'»ITEMPTY= (1 .0-PX (3) ) *l‘JTEfJlPTY 

WFIXED=WTFOOD+PASS+CRE*J+BAGGAGE+CARGO 

**************************************************** 

* THIS SECTION COMPUTES THE TAKE-OFF VJEIGHT * 
**************************************************** 

WTONEW=WrFIXED4^'JTFUEL-HiJTE.MPTY 

DELWTO=OTONEW-lYTO 

NOn>=NOIT+l 

IlYr=IWT+l 

IF(NOIT.GT.39) GO TO 58 
DIV=1.2 

IF(NOIT.LT.3) DIV=0.85 
IF(NOIT.GT.IO) DIV=1.95 

IF(IDEBUG.GT.O) VJRITE(6,338) NOIT,WrO,OTONEW,DELWTO 
$ ,lVrFUEL,WFOWTO 

338 FORMAT (5X*ITER. NO.=*I5,10X*WTO,l*JTONEJ,DELWTO,WrFUEL,WFOWTO=* 

$ 3F12.2,F12.2,F12.5) 
lVrO=V/rONEWfDELWTO/DIV 
DWTA (NOIT) =DELWTO 

************************************************ 

* ITERATE WEIGHT UNTIL I'JITHIN 0.05 LBS * 

******************************sV;V**************** 

IF(ABS(DELWTO) .GT.0.05) GO TO 50 

57 IF (OUTPUT. EQ.O) GO TO 150 
GO TO 60 

58 WRITE (6, 59) WrO,l'JTONEV‘J,DErj/JTO,VJTINIT,V/FOOTO,OTFUEL,NOIT 

59 FORMAT (/* / / / I'^IIGHT LOOP DID NOT CONVERGE/ / / * 

$ /10X*WrO,WTONEW,DELWTO,VJTINIT,l”JFOVJTO,OTFUEL,NOI'D=*4F12.2, 

$ F12.4,F12.2,I5/) 

WRITE (6, 64) DWTA 

64 FORMAT (10X*DELWTO*9 (T25,5F12.3/) ) 

WRITE (6 ,80) WIWING ,OTKT,l*JTVT,W’TFUSE,OTLG ,lVrCTL,WTINST,WIMISC,WTELEC 
$ ,WrSRT,WrFURN,lYT02 ,Vra”JIM,l'JTEGH,OTAC,WTTR,WTENG ,Ii;ECTL ,WTFSYS , ACTCON 
WRITE (6, 82) WTEMPTY 

WRITE ( 6 , 84 ) PASS , CRE!/'/ , BAGGAGE , WTFOOD , CARGO , VYTFIXED 
WRITE (6, 86) lVrFUEL,WTO,MOXT 



WRITE (6,88) WrHTl,WTHT2,WrHT3,WTFUSEl ,WTFUSE2,WTINIT 
CALL X0UTPUT(6) 

liJRITECerSS?) ICNT,JCNT,KCNT,NDAV,LCNT,IWT,NOIT 
887 FORMAT (/5X"DUMP OF FCOUNT COMMON BLOCK”/10X 
$ "ICNT, JCNT,NDAV,KCNT,LCNT,IWT,NOIT=" ,10X,7I5/) 

WTO=OTONEWfDELWTO/2 . 0 

OTS(l)=WTO 

GO TO 57 

*************************** 

* OUTPUT SECTION * 

*************************** 

60 WRITE(6,80)WIWING,VmiT,WTVT,WrFUSE,WrLG,WrCTL,WriNST,WmiSC,WTELEC 
$ ,WTSRT ,WTFURN , WT02 ,WMIN , WTBGH ,IaITAC jIaTITR, WTENG , WECTL ,WTFSYS , ACTCON 
80 FORrJlAT(*l*//30X,23H***WEIGHT ESTIMATION***//10X*WING*T40,F10.1 
$/10X*HOR. TAIL*T40,F10.1/10X*VERT. TAIL*T40,F10.1/10X*FUSELAGE*T40 
$,F10.1/10X*LANDING GEAR*T40 , FIO . 1/10X*CONTROL SYSTEM*T40,F10.1/10X 
$*INSTRUMENTS*T40 , FIO . 1/10X*MISC INTERIOR*T40 , FIO . 1/10X*ELECTRICAL* 
$T40 , FIO . 1/10X*STARTERS*T40 ,F10 . 1/10X*FURNISHINGS*T40 ,F10 . 1/lOX 
$*OXYGEN*T40 , FIO . 1/10X*V'JINDOWS*T40 , FIO . 1/10X*BAGGAGE HNDLING*T40 , 
$F10.1/10X*AIR CONDITIONING*T40,F10.1/10X*THRUST REVERSER*T40,F10.1 
$/10X*ENGINE*T40,F10.1/10X*ENGINE CONTROLS*T40,F10.1/10X*FUEL SYS* 
$T40,F10.1/10X*ACTIVE CONTROL SYSTEMS*T40,F10.1) 

VJRITE(6,82) WTEMPTY 
82 F0RiW(/5X*EMPTY WEIGHT*T40,F10.1) 

WRITE (6,84) PASS, CREIV, BAGGAGE, OTFOOD, CARGO, WTFIXED 
84 FORMAT (//10X*PASSENGERS*T40 , FIO . 1/10X*CREW*T40 , FIO . 1/lOX 
$*BAGGAGE*T40 , FIO . 1/10X*WTFOOD*T40 , FIO , 1/10X*CARGO*T40 , FIO . 1 
$//5X*FIXED WEIGHT*T40,F10.1) 

WRITE (6,86) WTFUEL,WrO,NOIT 

86 FORMAT (//5X*FUEL*T40, FIO. 1///5X*TAKE-0FF WEIGHT*T40,F10.1///10X 
$*N0. OF ITERATIONS REQUIRED*T40,I10) 

VJRITE(6,88) WrHTl,WTHT2,\fJTHT3,WTFUSEl,WTFUSE2,WTINIT 
88 F0RTJlAT(/5X*WrHT(l,2,3) ,WTFUSE(1,2) ,WINIT=*T45,5F10.2,F12.2) 

CALL FUELCAL(WrO,WTFUEL,l) 

CALL XL0D(GX(9) ,GX(3) ,WTS(11) ,EFF,1, OUTPUT) 
*************************************** 

* ASSIGN DATA BASE VARIABLES * 

*************************************** 

150 l'JTS(l)=WTO 

OTS b) =WTEWPTY 
WTS (4) =.75*WTEMPTY 
WTS(6)=WTFUEL 

IF(NOIT.LT.40) WFCWrO=OTFUEL/WTO 
AENG=0. 

IF(ITERM(4) .EQ.3) AENG=WTENG*0.01*LT*LT/(GX(15) *3.) 

WTS(13)=13.45E6 

OTS(14)=13.33E6 

CALL AT 62 (WrS(ll),ANS) 

WTS (7)=WTS (1) /(ANS(l) *W(10) *DESIGN (1) *GX (15) ) 

WTS (8 ) =WTS (13) *GX (15) /WTS (1) 

CALL AT62 (500., ANS) 

WTS (9) =GX (27) *WTS (1) /(GX (15) *DESIGN (1) *ANS (1) *W(10) ) 

WTS (10) =WTS (14) *GX (15) /(WTS (1) *GX (27) ) 
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IF(IDEBUG2.GT.O) WRITE(6,990) NOIT 
990 FORMAT (T100*NOn^*I10) 

RETURN 

END 

SUBROUTINE CRUALT(WT,CLCR,M,ALT) 

************************************************************************ 

* THIS SUBROUTINE RETURNS ALTITUDE TO SATISFY THE SPECIFIED * 

* IVEIGHT, CL, AND M. ITERATIVE TABLE LOOK UP IS USED. * 

************************************************************************ 

REAL M,HOLD,HNElfJ,ANS(4) 

COMMON /DEVARA>ESIGN(15) ,ITERM(10) ,CST(10) 

CCMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

IC=0 

DATA TO , PO, ALF,R/518 . 14 , 2116 , 229 00356617 , 53 . 3/ 

P=2 . 0*WT/ (1 . 4*DESIGN (1 ) *CLCR*M*M) 

HOLD= (1 . 0- (P/PO) ** (ALF*R) ) *TO/ALF 
50 CALL AT62(HOLD,ANS) 

DP=P-ANS(2) 

DRDH=-1.0*ANS(1) *GX(15) 

HNEW=HOLEH-DP*0 . 9/DRDH 
IC=IC+1 

IF(IC.GT.IOO) GO TO 100 
DALT=HNEW-HOLD 

IF(ABS(DALT).LT.7.5) GO TO 150 

HOLD=HNm 

GO TO 50 

100 WRITE (6, 102) HNEW 

102 FORMAT (*1*////40X,38H***CRUISE ALTITUDE DID NOT CONVERGE***//10X 
$*LAST ALTITUDE=*F12.2) 

150 ALT=HNEW 

152 FORMAT (///10X,*ALTITUDE, NO. OF ITERATIONS=*F12.2,I10) 

RETURN 

END 

SUBROUTINE ENGINE (ALT,M,TC'IM,TSFC) 

C 

C EMPIRICAL ADJUSTMENTS HAVE BEEN MADE TO THIS SUBROUTINE TO 
C MATCH DATA. 

C 

COHMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

REAL M,HN,ANS(4) 

CALL AT62 (ALT,ANS) 

DATA TO,PO/561.2,2116.229/ 

HN=ALT/40000. 

C 

C TCIM FROM REF (7) — AGGARWAL, ET. AL. 

C 

TCm= (1000./41100. )* (30. 06-34. 74*HN+7 .25*M4-12.11*HN*HN 
$ -1 . 12*M*HN-1 3 . 96*M*M-5 . 46*M*HN*HN+1 1 . 82*M*M*HN 
$ -.01*M*M*HN*HN) 

ALTK= ( (ALT-35332 .)* .004/1000 .) +1 . 

IF (ALT.GT. 35332) GO TO 40 
TCm=(1000./41100.)* (34. 71-32. 22*HN-27.4*M 
$ +5.8*HN*HN+36.43*HN*M+2.65*M*M-11 .64*M*HN*HN 
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$ +3.19*M*M*HN-5.72*M*M*HN*HN) 

ALTK=1. 

40 DELTA=ANS(2)/PO 
THETA=ANS(3) AO 
TCm=1.2*TC'IM 

TN=TC'm*41100 ./delta/20000 . 

FN=1 00 0 .*(2.04+3 . 71*TN-4 . 38*I4f 1 . 69*TN*TN+5 . 94*M*TN 
$ +12.99*M*M) 

FNS=FN*ALTK 

F=FNS*DELTA* (THETA) **0.5 
TSFC=F/ (TCIM*4 1100 . ) 

TSFC= (1 .0-PX (2) ) *TSFC 

RETURN 

END 

SUBROUTINE CDZL(ALT,CDO,M,OUTPUT) 

REAL MU(2) ,ANS(4) ,U(2) ,K(2) ,DCD(2) 

INTEGER OUTPUT 

COMiMON /DEVAR/DESIGN(15) ,ITERM(10) ,CST(10) 

COMMON /DRAG/CDS (6) ,CDSAP (6) 

COMMON /GEOiM/l'/(20) ,HX(20) ,GX(35) ,PX(15) 

CALL ATS2(ALT,ANS) 

DATA MU/3.6878E-7,2.9652E-7/ 

IALT=1 

IF (ALT.GT. 10000.) IALT=2 
SHT=DESIGN (4) 

TENG=DESIGN (6) /ITERM (4) 

U(1)=0.6*ANS(4) 

U(2)=M*ANS(4) 

SREF=DESIGN (1) 

********************************** 

* VJING DRAG CALCULATIONS * 

********************************** 

RE^W ( 1 0 ) *U ( I ALT) *ANS ( 1 ) /f^ ( lALT) 

CF=0.455/(ALOG10 (RE//) )**2.58 
CF=CF*(1.0-PX(4)) 

SV'/ET=2.05*SREF 
K(l)=1.2230 
K (2) =1.456 

CDa^J=CF*K ( I ALT) *SWET/SREF 
********************************************** 

* HORIZONTAL TAIL DRAG CALCULATIONS * 

********************************************** 

REH=HX (5) *U (lALT) *ANS (1) /MU (lALT) 
CF=0.074/REH**0.2 
K(l) =1.2560 
K (2) =1.3803 
Sl^/ET=SHT*2.05 
CDOH=CF*K ( lALT) *SWET/SREF 
******************************************* 

* VERTICAL TAIL DRAG CALCULATIONS * 
******************************************* 

TRV=HX(17) 

SV=GX(35) 



BV=(SV*HX(18))**0.5 

CVBAR=4 . *SV* (1 ,+TRV+TRV*TRV) / (3 . *BV* (1 ,+TRV) **2.0) 

REV=CVBAR*U (lALT) *ANS (1) /f4U (lALT) 

CP=.074/REV**0.2 

K(l)=1.1413 

K(2)=1.3504 

S'/JET=2.05*SV 

CDOV=CF*K (lALT) *SWET/SREF 
**************************************** 

* FUSELAGE DRAG CALCULATIONS * 

**************************************** 

REF=DESIGN ( 3 ) *U ( lALT) *ANS ( 1 ) AlU ( lALT) 

CF=0 .0455/ (ALOGIO (REF) ) **2 . 58 
DIA=GX{5) 

FL0D=DESIGN(3)/DIA 

RDIA=DIA/2.0 

ANOSE=2 . *GX (1 ) *RDIA*RDIA 

CONEL= (DESIGN ( 3 ) -RDIA) * . 25 

ABOD=GX (1) *GX (5) * (DESIGN (3) -RDIA-CONEL) 

SLANT= (CONEL*CONEL+RDIA*RDIA) **0 . 5 
ACONE=GX (2) *RDIA*SLANT 
SFUS=ANOSE+ABOD-i-ACONE 

RAT^ (1 .+60 ./ (FLOD*FLOD*FLOD) +. 0025*FLOD) *SFUS/50 . 

CDFFi=CF*RAT 

CDBF=5.0112E-5/(CDFF) **0.5 
CDOF=CDFF*50 . 0/SREF 

******************************************** 

* ENGINE NACELLE DRAG CALCULATIONS * 
******************************************** 

EL=GX (6 ) * (TENG/GX (21) ) **0 . 5 
EDIA=GX (7) * (TENG/GX (21) ) **0 . 5 
SWET=GX (1 ) *EDIA*EL 
REE=EL*U (lALT) *ANS (1) /MU (I ALT) 

CF^0.074/REE**0.2 

TT=2 

IF(ITERM(4) .GT.3) TT=4 
CDOE=l . 0033*CF*SVffiT*TT/SREF 

************************************************************************ 

* SUM UP DRAG COMPONENTS * 

* * 

* CRUD DRAG (IAL'1^2) * 

* LANDING GEAR DOWN (0.0150) + 10 DEGREES FLAPS (0.008) * 

* (IALT=1) * 

************************************************************************ 

DCD(1)=0.023 
DCD (2) =0.0007 

CDO= (CDOW+CDOH+CDOV+CDOV+CDOF+CDBF+CDOE+DCD (lALT) ) *1 .05 
IF(IALT.NE.l) GO TO 76 
********************************** 

* DATA BASE ASSIGNMENTS * 

********************************** 

CDSAP(l)=CDO 

CDSAP(2)=CDOW 



CDSAP(3)=CD0H 
GO TO 78 
76 CDS(l)=CDO 
CDS(2)=CDCF/tf 
CDS(3)=CDOH 
CDS(6)=CDOF 

************************** 


* OUTPUT SECTION * 
************************** 


78 IF(OUTPUT.EQ.O) GO TO 100 

VJRITE(6,80) CDOW,CDOH,CDOV,CDOF,CDBF,CDOE,DCD(IALT) ,CDO 
80 F0PT4AT(*1V/15X,19H***DRK3 ANALYSIS*** ,//10X*WING*T35,F10.4/10X 
$*HORIZONTAL TAIL*T35,F10.4/10X*VERTICAL TAIL*T35,F10.4/10X*FUSELAG 
$E*T35 , FIO . 4/10X*BASE*T35 , FIO . 4/10)C*ENGINE NACELLE*T35 , FIO . 4/lOX 
$ *CRUD/FLAPS*T35 , FIO . 4//5X 

$*AIRCRAFT DRAG*T35,F10.4///10X*INTERFERENCE FACTOR IS 5 PERCENT*) 
l'JRITE(6,82) REW,REH,REV,REF,REE 

82 FORMAT (////15X , 23H***REYNOLD'’ S NUMBERS***//10X*WING*T35 , F12 . 1/lOX , 
$*HORIZONTAL TAIL*T35,F12.1/10X*VERTICAL TAIL*T35,F12.1/10X*FUSELAG 
$E*T35 r F12 . 1/10X*ENGINE*T35 , F12 . 1//) 

IF(IALT.EQ.l) WRITE(6,84) 

84 FORI-1AT (///30X,14H***APPROACH***/) 

100 RETURN 
END 

SUBROUTINE STABCOD(CL,M,ICG,OUTPUT) 


INTEGER OUTPUT 

REAL DCDM(2) ,XACDM(2) ,ANS(4) 

REAL KA,KTR,KH,M,KWB,K 


REAL LF 

CO^^^^ON /DEVAR/DESIGN (15) ,ITERM(10) ,CST(10) 
Ca^iON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COr-MON /GRAVITY/OG(6) 

COrWN /STABA>ERIVCR(15) ,DERIVAP(15) ,STOR(20) 


COMT-CT /WTSVE/WTS(20) 


* THIS SUBROUTINE CALCULATES STABILITY AND CONTROL DERIVATIVES * 

* IN EITHER APPROACH OR CRUISE CONFIGURATIONS * 

* ICR=0 (APPROACH); ICR=1 (CRUISE) 

^^^^^^^^^^^****;Uf***************************************^^********’^***** 


XCG=(OG(1)+OG(3))/2.0 
IF(IOG.LT.7) XCG=OG(ICG) 

ICR=1 

ELASTK=0.875 

IF(M.LT.0.45) ELASTK=0.825 
IF(M.LT.0.46) ICR=0 

****************************************************** 

* THE FOLLOWING SECTION CALCULATES CL (ALPHA) * 
****************************************************** 

KWB=1.0-0.25*GX(5)*GX(5)/(W(9)*W(9))+0.025*GX(5)A^(9) 
BETA= (1 . 0-M*M) **0 . 5 
PIAR=DESIGN (2) *GX (1) 

K=1.0 

XNUM=2.0*PIAR 
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T1=DESIGN (2) *DESIGN (2) *BETA*BETA/ (K*K) 

T2=1.0-W(13) *W(13) /(BETA*BETA) 

CLA!flf=XNUM*ELASTK/ (2 » 0+ (Tl*T2+4 .0) **0 , 5) 

T1=DESIGN (2) *DESIGN (2) / (K*K) 

T2=1.0+W{13)*W{13) 

CLAWM0=XNUM* .725/ (2.0+ (Tl*T2+4 . 0) **0 . 5) 

KA=1 ./DESIGN (2) -1 .0/ (1 .+DESIGN (2) **1.7) 
KTR=(10.-3.*W{3))/7. 

ZTL=GX(34) 

XLT=DESIGN (1) *W(10) *HX (6) /DESIGN (4) 

IOl=(1.0-ZTL/W(9) ) /(2.0*XLTAH9) ) **0.40 
DEDAM0=4.44* (KA*KTR*KH* (W(15) ) **2.0) **1.119 
DEDA=DEDAM0* (CLAW/CLA'/MO) **0.444444 
*1^1= (1.0+HX (9) *HX (9) / (BETA*BETA) ) 

'IM2=BETA*BETA*DESIGN (5) *DESIQ1 (5) 

CLAH=2 . 0*GX (1) *ELASTK*DESIQ1 (5) / (2 . 0+ (TM2*TMl+4 .0) **0 . 5) 

CLAWB=KIVB*CLAW 

ETAH=0.90 

CLA=CLAWB+CLAH*ETAH*DESIGN (4) * (1 .0-DEDA) /DESIGN (1) 
Elfn:NG=W(7) 

* THE FOLLOWING STATET^ENT CALCULATES CD (ALPHA) * 

******************************f?************************** 

CDA=2 . 0*CL*CLA/ (PIAU*E7/JIKG) 

****************************************************** 

* THE FOLLOWING SECTION CALCULATES CM (ALPHA) * 

*******************************5‘t********************** 

)CAaAN.25-.265*M*M/ (DESIGN (2) ) **0.5 
DXACT=XLT/W(10) 

A!fJ=CLAW/GX(2) 

XNRM=DESIGN (1) *W(10) *AW 
TFUS=0. 5836+. 1690* (1.0-DEDA) 

DXACF=-GX (5) *GX (5) *DSSIGN (3) *TFUS/(XNRM*36.5) 

DM'aDF=GX (5) *W(10) / (DESIGN] (1) *AIfJ*290, ) 

DXACE^DXACF+DMWDF 
STOR(ll)=DXACF*AIfJ 
T=DESIGN (6) /ITERM (4) 

VJF=GX (6) * (T/GX (21) ) **0 . 5 
LF=GX(7) * (T/GX (21) ) **0.5 
TT=2. 

IF(ITERM(4) .GT.3) TT=4 
DXACE=-TT*LF*WF*WF/ (15 . 21*XMRM) 

STOR(12)=DXACE*A!fJ 

TA1=CLAH*ETAH*DESIGN (4) * (1 .0-DEDA) /(CLAWB*DESIGN (1) ) 

XACWB=XACWfDXACF+DXACE 

XAC= (XACWB+TA1*DXACT) / (1+TAl) 

DCMCL=XCG-XAC 

CMA=DCMCL*CLA 

A**********************************************’*'*************** 

* THE FOLLOWING SECTION CALCULATES -VELOCITY DERIVATIVES * 

<Hk******************************************A5***************** 

CLU=(li+M)'^M*M*GL/( >-•- 

DATA DCDM,XACDM/0.0202,0.1005,-3.1E-4,5.1E-4/ 
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z=o. 

IF(ICR.GT.O) Z=40000. 

CALL AT62(Z,ANS) 

AT4=ANS(4) 

CDU=DCDM ( ICR+1 ) 2 . *CL*CLU/ (PIAR*EWING*10 . ) 

CKU=-CL*XACDM (ICR+1 ) 

********************************************************* 

* THE FOLLOWING SECTION CALCULATES Q DERIVATIVES * 
********************************************************* 

CLQH=2 . 0*CLAH*ETAH*HX (6 ) 

CLQ=CLQH*1.1 

CDQ=0 

CMQT=-2 .0*CLAH*ETAH*0 ,9*HX (6) *XLTA^ (10) 

*********vV******************************************************** 

* THE FOLLOWING SECTION CALCULATES ALPHA-DOT DERIVATIVES * 

*********************************, V******************************** 


a^O=Cf4QT*l.l 

CLADOT=-2 . 0*CLAH*ETAH*HX (6) *DEDA 
CDADOT=0.0 

a4ADOT=-2 . 0*CLAH*ETAH*HX (6) *XLT*DEDAA^ (10) 
**************************************************************** 

* THE FOLLOWING SECTION CALCULATES CONTROL DERIVATIVES * 
**************************************************************** 

CLDEL=0 . 9*DESIGN (4 ) *CLAH*0 o 46/DESI@I (1 ) 

* THIS IS PER DEGREE CONTROL DEFLECTION, OTHER PER RAD. * 

************************************************************************ 


CDDEL=DESIGN (4) *1 .980E-3/DESIGN (1) 
CKDEL=-CLAH*HX (6 ) *0 . 9*0 , 46 
IF(ICR.LT.l) GO TO 120 
Xr-lU=215. 

IF (OUTPUT.GT.-l ) XMU=WTS (7) 
*********************************** 

* DATA BLOCK ASSIGNMENTS * 

*******;>*************************** 


STOR(l)=CLAW 

STOR(2)=CLAH 

STOR(3)=DEDA 

STOR(4)=XAC 

STOR (5 ) =XAC-CMQ/ (4 . *XMU) 

ST0R(13)=TA1 

DERIVCR(1)=CLA 

DSRIVCR(2)=CDA 

DERIVCR(3)=CMA 

DERIVCR(4)=CLU 

DERIVCR(5)=CDU 

DERIVCR(6)=CMU 

DERIVCR(7)=CLQ 

DERIVCR(8)=CDQ 

DSRIVCR(9)=CMQ 

DERIVCR (10) =CLADOT 

DERIVCR (11) =CDAD0T 

DERIVCR (12) =CMADOT 
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DERIVCR(13)=CLDEL 
DERIVCR(14)=CDDEL 
DERIVCR(15)=CMDEL 
GO TO 150 
120 DERIVAP(1)=CLA 
DERIVAP(2)=CDA 
DERIVAP(3)=CMA 
DERIVAP(4)=CLU 
DERIVAP(5)=CDU 
DERIVAP(6)=CMU 
DERIVAP(7)=CLQ 
DERIVAP(8)=CDQ 
DERIVAP(9)=CMQ 
DSRIVAP(10)=CLADOT 
DERIVAP(ll)=CDADOT 
DERIVAP (12) =CMADOT 
DSRIVAP(13)=CLDEL 
DERIVAP (14) =CDDEL 
DERIVAP (15) =CMDEL 
XMU=70. 

IF(OUTPUT.GT.-l) XMU=l‘JTS(9) 

STOR(6)=CLA 

STOR(7)=CLAH 

STOR(8)=DEDA 

STOR(9)=XAC 

STOR (10) =XAC-CMQ/ (4 . *XMU) 

ST0R(14)=TA1 

150 IF (OUTPUT. LT.l) GO 'TO 200 
************************** 

* OUTPUT SECTION * 

*******;V************ ****** 

IaIRITE(6,162) 

162 FOR''1AT(*1V/20X*STABILITY AND CONTROL DERIVATIVES*//) 
IF(ICR.EQ.O) WRITE(6,163) 

163 FORMAT (25X,14H***APPROACH***//) 

VJRITE(6,164) CL,M,XCG 

164 FORMAT (10X*CL,M,CG. POSITION=*3F15.3//) 

VJRITE(6,166) 

166 FORMAT (/28X*CL*14X*CD*12X*CMV) 

IF(ICR.EQ.l) GO TO 180 
WRITE (6, 174) DERIVAP 

174 FORMAT ( 1 0X*ALPHA*4X , 3F15 . 5/10X*VELOCITY*3F15 . 4/10X*Q*8X , 3F15 . 5/ 
$10X*ALPHA-DOT*3F15 . 5/10JC*ELEVATOR *3F15 . 5//) 

GO TO 190 

180 WRITE (6, 174) DERIVCR 

190 VJRITE(6,192) XAC,DCMCL,CMA,STOR 

192 FORMAT (/10X*NEUTRAL POINT*T35,F12.3/10X*STATIC STABILITY* 

$T35 ,F12 . 3/10X*CM (ALPHA) *T35 rF12.3/10X*STOR=*/ 

$2 (T35f 2F12.3/) ) 

WRITE(6 ,196) DXACE,XACWB,XACW,AW,XNRM,TFUS,DXACF,DMWDF, 

$ CLAWB,CLAH,DEDAM0,XLT,ZTL,CLAWM0,BETA,TA1,CLAW,DEDA 
$ ,XAC,KH,KTR,'IMl,TM2,Kli!3,X0G,KA 
196 FORMAT(5(T30,5F12.4/)) 



200 RETURN 
END 

SUBROUTINE TRIM(CL,M,CDO,CLWING,CLTAIL,CRIT,CRDE,EPS,IPHASE 
$,ICG) 

REAL LTOTOQ,LWOQ,LTLOQ,IW,M 

COMMON /DEVAR/DESIGN(15) ,ITERM(10) ,CST(10) 

COMON /DRAG/CDS (6) ,CDSAP (6) 

COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /GRAVITY/0G(6) 

COITION /STAB/DERIVCR(15) ,DERIVAP(15) ,STOR(20) 
************************************************************************ 

* IPHASE=1 (CR) , =0 (10 DEG FLP,R® CG) , =-l (45 DEG FLP) , * 

* =~2 (45 DEG FLAP, EVJD CG) , =-3 (45 DEG FLAP, EV/D CG, FIND * 

* CLTW) , =-4 (10 DEG FLAP, ElO CG, FIND CLMAX) , =2 (CR, AFT CG) * 

************************************************************************ 

N=0 

************************************************ 

* INITIALIZE VARIABLES FROM DATA BASE * 

************************************************ 

XA= . 25- . 265*M*M/ (DESIGN (2) ) **0 . 5 

XAC= (CG (1 ) +CG ( 3 ) ) /2 .-XA 

IF(ICG.LT.7) XAC=OG(ICG)-XA 

ZAC=0.08 

ZTLC=-.12 

ST=DESIGN(4) 

a4AG^=-.10-.030218*W(ll)-.046875*M 
CMACW=CMACW* (1 .0-PX (1) ) 

CDT=CDS (3) 

AT=STOR(2)/GX(2) 

AV.^STOR(l)/GX(2) 

DEDA=STOR(3) 

AOL=W(16) 

IF(IPHASE.GT.O) GO TO 10 
AT=STOR(7)/GX(2) 

Al*J=STOR(6)/GX(2) 

DEDA=STOR(8) 

CDT=CDSAP(3) 

CMAa*j=CMAa;+w ( 17 ) 

W(19)=CMAOV-W(17) 

AOL=W(16)-K3X(17) 

IF(IPHASE.LT.-3) GO TO 10 
IF(IPHASE.EQ.O) GO TO 10 
Cf4AO/;=CMAO//-W(18) -W(17) 

AOL=W(16)+GX(18) 

CDO=CDO+W(20) 

10 CD=CDO+CL*CL/ (DESIGN (2) *GX(1)*W(7)) 

CRAFUS=STOR (ll)+STOR(12) 

IW=W(2) 

IF(IPHASE.LT.-2) GO TO 200 
LTOTOQ=CL*DESIGN (1 ) 

LTLOQ=-.l*LTOTOQ 

************************************************************************** 
* THIS SECTION C/\LCULATES CL(WIEK3) AND CL (TAIL) GIVEN CL (TOTAL) * 
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non 


************************************************************************** 
20 LW0Q=LT0T0Q-LTL0Q 
CLWING=LWOQ/t)ESIGN ( 1 ) 

IF(CLWING.GT.3.5) CLWING=3.5 
IF(CLWING.LT.O.l) CLWING=0.1 
CLTAIL=LTLOQ/ST 

IF(ABS(CLTAIL) .GT.2.0) CLTAIL=2.0*CLTAIL/ABS(CLTAIL) 

CDTAIL=CDT*DESIGN (1) /ST+CLTML*CLTAIL/ (GX (1) *DESIGN (5) *0.90) 

ALFAW=CLWING/AWfAOL 

IF(ALFAW.GT.16.0) ALFAW=16.0 

CMFUS=CMAFUS* (ALFAi*MW) 

EPS=CLWING*DEDA/AW 

CRDE=0, 

ALFAT=CLTAIL/AT-CRDE 

CTLPWR= (CMACW+CLWING*XAC+CMFUS-AT* (ALFAW-EPS-IW) *HX (6) *.9) / 
$(AT*HX(6)*0.9) 

TAO=0.46 

CRIT=CTLPWR 

CRDE=0. 

IF(CRIT.LT.5.0) GO TO 40 
CRDE= {CTLPWR-5 . 0 ) AAO 
CRIT=5.0 
GO TO 60 

-40 IF(CRIT.GT.-14.0) GO TO 60 
CRDE= (CTLPWR+14 .0) AAO - 
CRIT=-14.0 

60 CNlf;=CLWING*COS ( ( ALFAi-J-IW) /GX (2 ) ) +CD*SIN ( (ALFAW-IW) / 

$GX(2)) 

CCW=CD*COS ( (ALFAW-IW) /GX (2) ) -CLWING*SIN ( (ALFAW-IW) / 

$GX(2)) 

CT^CD 


IF TRIJET, THE UNBALANCE IS ONLY ONE ENGINE 
IF(ITERM(4) .EQ.3) CT=CT/3.0 

CNT= (1 ./ (0 . 9*HX (6 ) ) ) * (CNlrJ*XAC4GaAr*ZACK34FUS-CT*ZTLC4CMACW) 
CLTOEW= (CNT-CDTAIL*SIN ( (ALFAT-CRIT) /GX (2) ) ) /(COS ( (ALFAT 
$-CRIT)/GX(2))) 

DTEST=CLTNEW-CLTAIL 
CLTAIL=CL'INEWfDTEST/ (N+1 ) 

IF (ABS (CLTAIL) .GT.2.0) CLTAIL=2 . 0*CLTAIL/ABS (CLTAIL) 

LTLOQ=CLTAIL*ST 

N=N+1 

IF(ABS(DTEST) .LT.0.003) GO TO 100 
IP(N.GT.25) GO TO GO 
GO TO 20 

80 WRITE(6,82) CLWING,CLTAIL,CRIT,CRDE,CNW,DTEST 
$,ALFAW,EPS,CMFUS 

82 FORMAT,<///20X*TRIM DID NOT C0NVERGE*//3(3F15.4/)//) 

100 IF(CLWING.EQ.3.5) WRITE(6,84) CLWING 
$,CLTAIL,CRIT,CRDE,CNW,DTEST ‘ 

IF(ALFAW.EQ.16.0) l"3RITE(6,84) CLWING 
$, CLTAIL, CRIT,CRDE,CNW,DTEST . 



84 FORMAT (///30X*TRIM HIT ALPHA OR CL LIMIT*//2(10X,3F12.4/) ) 
IF(ALFAW.EQ.16.0) CALL XOUTPUT(6) 

RETURN 

************************************************************************** 
* THIS SECTION CALCULATES CL (TAIL) AND CL (TOTAL) GIVEN CL (WING) * 

************************************************************************** 

200 CLVJING=CL 
CLTAIL=— 1 • 

CDTAIL=CDT*DESIGN (1) /ST4CLTAIL*CLTAIL/(GX(1) *DESIGN (5) *0.90) 

ALFAW=CLWING/Ai^fAOL 

IF(ALFAW.GT.16.0) ALFAW=16.0 

CtJlFUS=CMAFUS* (ALFAW-IW) 

EPS=CLWING*DEDA/AW 

CRDE=0. 

ALFAT=CLTAIL/AT-CRDE 

CTLPWR= (CMACW+CLWING*XAC4CMFUS-AT* (ALFAW-EPS-IW) *HX (6) * .9) / 
$(AT*HX(6)*0.9) 

TAO=0.46 

CRIT=CTLPWR 

CRDS=0. 

IF(CRIT.LT.5.0) GO TO 540 
CRDE= (CTLPWR-5 . 0 ) AAO 
CRIT=5.0 
GO TO 560 

540 IF(CRIT.GT.-14.0) GO TO 560 
CRDE= (CTLPWR+14 ,0) /TAO 
CRIT=-14.0 

560 CN'/J=CLWING*COS ( (ALFAW-IW) /GX (2 ) ) +CD*SIN ( (ALFAW-IW) / 

$GX(2)) 

Caf=CD*COS ( (ALFAW-IW) /GX (2) ) -CLWING*SIN ( (ALFAW-IW) / 

$GX(2)) 

CT=CD 

IF(ITERM(4) .EQ.3) CT=CT/3.0 

CNT= (1 ./ (0 .9*HX (6) ) ) * (CNW*XAC4CCW*ZAC+CMFUS-CT*ZTLC-K:MACW) 

CLTNEVJ= (CNT-CDTAIL*SIN ( ( ALFAT-CRIT) /GX (2) ) ) / (COS ( (ALFAT 
$-CRIT)/GX(2))) 

CLTAIL=CLTNEW 
LTLOQ=CLTAIL*ST 
LWOQ=CLWING*DESIGN (1) 

LTOTOQ=LTLOQfLWOQ 
CL=LTOTOQ/DESIGN ( 1 ) 

RETURN 

END 

SUBROUTINE FUELCAL (WTO , WTFUEL , OUTPUT) 

COf'lMON /DEVAR/DESIGN (15) , ITERM (10) ,CST (10) 

COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /WTSVE/WTS(20) 

REAL ANSI (4) 

INTEGER OUTPUT 
IF (OUTPUT.lt. 1) GO TO 30 
WRITE (6,21) 

21 FORMAT (*1*///30X*/ / / CRUISE ANALYSIS/ / / *//) 

WRITE (6,23) RANGE 


U8 



23 FORMAT (10X*TOTAL MISSION RANGE=*F10.2/10X 
$ *CLIMB DISTANCE= 189.00*/10X 
$ *DESCENT DISTANCE= 113.00*/) 

VmiTE(6,25) 

25 FORMAT (//2X*LBG*2X*MACH NO.*3X*CL*4X*V*7X*TIME*4X*L/D* 

$ 4X*TSFC*2X*T/T (IN) *1X*ALT (BEG) *2X*ALT(END) *2X*WT(BBG) * 

$ 3X*WT(END) *6X*DIST*4X*TC0N*4X*GAT-1A*5X*CLM*/) 

30 l'gTBEG=lVrO*0. 97*0. 965 
GX(31)=.0020 
RANGE=GX(4) 

GX(30)=5. 

TOTTIME=0.78 

ALT1=36000. 

ALT11=ALT1 
VJT1=WTBEG 
RLOVER=RANGE-302 . 

R=RLOVER/10. 

**************************************************** 

* CALCULATE CRUISE PORTION IN 10 SEGMENTS * 
**************************************************** 

DO 40 1=1,10 

CALL CRUFUEL ( WTl , R , WT2 , ALTl , ALT2 , 1 , TIME , OUTPUT , 0 ) 

VJT1=WT2 

ALT1=ALT2 

TOTTIME=TOTTIME+TIME 
IF(I.LT.5) GO TO 40 
IF(I.GT.5) GO TO 40 
VnT4ID=WT2 
OTS(15)=WmiD 
lirrS(ll)=ALT2 
40 CONTINUE 

************************************** 

* RESERVE MISSION ASSUf^PTION * 
************************************** 

VnGRD=WT2 

CALL CRUFUEL (WIGRD, 1000 . , WIRES ,30000 . ,ALTR, 11 ,TM, OUTPUT, 1) 
************************************ 

* INSERT DATA IN DATA BASE * 
************************************ 


FL=TOTTIME 

BLKHR=FL+.327 

SPEED=RANGE/FL 

BU<SPD=RANGE/BLKHR 

WTS(5)=.97*WIGRD 

WTEND=0.97*WTRES 

WTFUEL=WTO-WTEND 

l'JTS(17)=SPEED 

WTS(18)=FL 

WTS(19)=BLKHR 

WTS (20) =VJTBBG-WTS (5) 

WTS (3) =WTS (20) /BLKHR 

WTS (3 ) = (WTBEG-WTS ( 5 ) ) AIME 

FUEL=WTS(20)/6.4 
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XMG=RANGE/FUEL 

XSMG=XMG*GX(19) 

CST(5)=XSMG 

*** A*** A*****^t ************ 

* OUTPUT SECTION * 

************************** 

IF (OUTPUT. EQ.O) GO TO 100 

VaiTE(6,52) WrO,WrBBG,WIMID,WT3RD,WrRES,WrEND,WrFUEL 
52 FORMAT (*1*//40X,26H***FUEL WEIGHT ANALYSIS***///10X*TAKE-OFF*T35, 
$F12 . 2/10X*START-CRUISE*T35 , F12 . 2/10X*MID-CRUISE*T35 , F12 . 2/ 
$10X*END-CRUISE*T35,F12.2/10X*AFTER RESERVE*T35,F12.2/10X 
$*AFTER DESCENT/TAXI*T35,F12.2//5X*NET FUEL WEIGHT (LBS) *T35,F12. 2) 
WRITE (6, 54) ALT11,ALT2,ALTR 

54 FORMAT (////5X*CRUISE ALTITUDES*/! 0X*LEG 1*T35,F12.2/10X*LBG 2* 

$T35 , F12 . 2/10X*RESERVE LEG*T35 , F12 . 2) 
l*JRITE(6,64) FL, SPEED, BLKHR,BLKSPD 
64 FORT>1AT(///10X*FLIGHT LENGTH (HR) *T35,F12.2/ 

$ 10X*AVERAGE SPEED (KTS) *T35,F12.2/ 

$ 10X*BLOCK TIME (HR) *T35,F12.2/ 

$ 10X*BLOCK SPEED (KTS) *T35,F12.2) 
l"JRITE(6,66) V7TS(20) ,FUEL,XMG,XS^'!G 
66 FORT4AT(10X*BLOCK FUEL (LBS) *T35,F12.2/ 

$10X*BLOCK FUEL (GALS) *T35,F12.2/10X 

$*NAUT. MI/GAL*T35,F12.2/10X*NAUT. SEAT MI ./GAL.* 

$T35,F12.2) 

TENG=DESIGN (6) /ITERM (4) 

TREF=GX(21) 

SCP=TENG/TREF 

WRITE(6,68) DESIGN (6) ,ITERM(4) ,TENG,TREF,SCF 
68 FORT-1AT(/10X*INSTALLED THRUST (LBS) *T35,F12.2/ 

$ 10X*NO. OF ENGINES*T35, 112/1 OX 
$ *ENGINE THRUST (LBS) *T35,F12.2/10X 
$ *REFERENCE ENGINE (LBS) *T35,F12.2/10X 
$ *SCALE FACTOR*T35,F12.3) 

100 RETURN 
END 

SUBROUTINE CRUFUEL (WTBBG , RANGE ,WTEND,ALTCR,ALTEND, ICOUNT, TIME 
$ , OUTPUT, IRES) 

************************************************************************ 

* CRUFUEL CALCULATES PERFORI4ANCE DURING CRUISE/CLIMB OR RESERVE * 

* SEGMENT * 

************************************************************************ 

REAL ANS(4) ,LOD,HZ,M 
INTEGER OUTPUT 

COM-'iON /DEVAR/DESIGN(15) ,ITERM(10) ,CST(10) 

COMMON /DRAG/CDS (6) ,CDSAP (6) 

COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COrWN /WTSVE/WS(20) 

M=GX(3) 

GAT4A=0. 

CALL CDZL(ALTCR,CDO,M>0) 

CD0=CD04GX(31) 

XT1=W(8) *W(8) * (1.0/W(6)-1.0/W(7) ) /(GX (1) *DESIGN (2) ) 
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CLM=1 .00* (GX (1) *DESIGN (2) *W(7) * (CDO+XTl) ) **0 .5 

IF(IRES.LT.l) GO TO 10 

CLCR=0.79*CLf’l 

************************************************************************ 

* the following statements calculate the reserve mission which * 

* IS AT 30000 FT AT MAXIMUT4 RANGE. * 

************************************************************************ 

CALL AT62 (30000., ANS) 

VCR= (WIBBG*2.0/ (DESIGN (1) *ANS(1) *CLCR) ) **0.5 
M=VCR/ANS(4) 

CALL XLOD(CLCR,M, 30000., LOD, 1,0) 

HZ=3.0 

ALTEND=30000. 

TIME=RANGE/ (VCR*1 . 467*1 . 1507 ) 

GO TO 88 

************************************************************************ 

* THE FOLLOWING STATEMENTS CALCULATE THE CRUISE/CLIMB AT MAXIMUM * 

* RAtXSE FACTOR. * 

************************************************************************ 

10 CALL AT62(ALTCR,ANS) 

VCR=M*ANS(4) 

WrT=.98*WTBEG 

0=0 . 5*ANS (1) *VCR*VCR 

CLCR=WrBBG/ (Q*DESIGN (1 ) ) 

IF(CLCR.GT.0.79*CLM) GO TO 15 
CLCR=0.79*CLM 

VCR= (WTBEG*2./(ANS (1) *DESIGN (1) *CLCR) ) **0.5 
M=VCR/ANS(4) 

CRCL=.79*CLM 

TIME=RANGE/(VCR/(1. 467*1. 1507) ) 

CALL CRUALT(WTT,CRCL,GX(3) ,ALTD) 

ROC= (ALTD-ALTCR) /(3600.* (TIME-. 02) ) 

GAMA=ASIN (ROCACR) 

GO TO 20 

*************************************************************************** 

* THE FOLLOWING STATEMEISTTS CALCULATE THE CRUISE/CLIMB AT 90% OF * 

* CL FOR (L/D)MAX. * 

*************************************************************************** 

15 IF(CLCR.LT.0.98*CLM) GO TO 18 * 

CRCL=0.9*CLM * 

CALL CRUALT(WTT,CRCL,GX(3) ,ALTD) * 

TIME=RANGE/ (VCR/ (1 . 467*1 , 1507) ) * 

ROC= (ALTD-ALTCR) /(3600.* (TIME-. 02) ) * 

GAMA=ASIN (ROC/VCR) * 

GO TO 20 * 

*************************************************************************** 

18 IF (ICOUNT.lt. 6) GO TO 20 
CRCL=.9*CLM 

CALL CRUALT(WTT,CRCL,GX(3),ALTD) 

ROC= (ALTD-ALTCR) /(3600.* (TIME+0. 02) ) 

GAMA=ASIN (ROC/VCR) 

20 CALL XLOD(CLCR,M,ALTCR,LOD,1,0) 

HZ=ALTCR/10000. 
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88 CALL ENGINE (ALTCR,M,TCIM,TSFC) 

WTEND=WTBE1G/ (EXP (TSFC*RANGE/ (VCR*0 . 59239*LOD) ) ) 

TIME=RANGE/ (VCR/ (1.467*1 . 1507) ) 

T0WAV=DESIGN (6 ) /WTBBG 

T0WRQ= (1 ./TCIM) * (1 . AOD+SIN (GAMA) ) 

TCON==raWAV/TOWRQ 
IF (TCON.LT.GX (30) ) GX (30) =TCON 
IF(ICOUNT.GT.IO) GO TO 92 
CALL CRUALT(WTEND,CLCR,M,ALTEND) 

ALTEND=ALTEND+SIN (GAMA) *VCR*3600 . *TIME 

92 IF(ICOUNT.NE.5) GO TO 99 
*******************:!:************* 

* DATA BASE ASSIGNMENTS * 

idrkie'kitiziciic’Mcfcitfrkic'kitic'k'kirkirkicifit'k’kic'kit 

CDS(4)=LOD 

CDS(5)=CLCRAOD 

GX(9)=CLCR 

CALL CRUALT(WTEND,CLM,GX(3) ,ALTLDM) 

WTS(12)=ALTLDM 

99 IF (OUTPUT. LT.l) GO TO 100 
************************** 

* OUTPUT SECTION * 

************************** 

VJRITE(6,188) ICOUNTrM^CLCR^VCRAIMErLOD^TSFC^TCIM, 

$ ALTCR,ALTEND,WTBEG,OTE!S!D,RANGE,TCON,GAMA,CLr4 
188 FORT4AT(I5,F8.2,F7.3,F8.1,F8.3,F8.2,F7.3,F7.3,5F10.1, 

$ F8.3,F8.4,F8.3) 

100 RETURN 
END 

SUBROUTINE XLOD(CL,XM,ALT,LOD, IPHASE, OUTPUT) 

REAL LOD,ANS(4) ,A1(5) ,A2(4) ,A3(5) ,B1(4) 

INTEGER OUTPUT 

COr^MON /DEVAR/DESIGM (15) ,ITERM(10) ,CST(10) 

COMMON ARAG/CDS(6) ,CDSAP(6) 

COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /GRAVITY/0G(6) 

COMON /STAB/DERIVCR(15) ,DERIVAP(15) ,STOR(20) 

COMI'ION /WTSVE/WTS(20) 

************************************************************************ 

* IPHASE=1 (CR) , =0 (10 DEG FLP,Fb.D CG) , =-l (45 DEG FLP) , * 

* =-2 (45 DEG FLAP, FIaID CG) , =-3 (45 DEG FLAP, EWD OG, FIND * 

* CLT^) , =-4 (10 DEG FLAP, B® CG, FIND CLMAX) , =2 (CR, AFT CG) , * 

************************************************************************ 

DATA Al/. 00102 , . 028817 ,.841,5.714076, 10.706253/ 

DATA A2/1.00131,-. 122063, .030714,-.005556/ 

DATA A3/276 . 559019 ,-1306 . 362579 , 2314 . 351122,-1817 . 552645 , 

$ 534.505085/ 

DATA Bl/-.794,-.296,.812,.lll/ 

IF(IPHi^E.LT.l) XM=.15 
ICG=ITERM(3) 

IF(IPHASE.GT.l) IOG=l 
IF(IPHASE.GT.O) GO TO 25 
CALL AT62(ALT,ANS) 
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F1=GX(27) 

IF(IPHASE.EQ.O.OR.IPHASE.LT.-3) Fl=l. 

V= (2.*F1*IVTS (1) /{ANS (1) *DESIGN (1) *CL) ) **0.5 
XM=V/ANS(4) 

ICG=4 

• IF(IPHASE.EQ.-l.OR.IPHASE.EQ.O) IOG=2 
ICG=2 

IF(IPHASE.LT.-l) ICG=4 
25 CALL CDZL (ALT, CDO,XM, OUTPUT) 

CALL STABCOD (CL, XM, ICG, OUTPUT) 

CALL TRIM ( CL , XM , CDO , CLWING , CLTAIL , CRIT , CRDE , EPS , IPHASE , 
$ICG) 

*********************************** 

* CQMPRESSIBILITY DRAG * 

*********************************** 

DFME=0. 

FKCR=0o 

DWME=0. 

WMCR=0. 

IF(XM.LT.0.65) GO TO 40 
CLWN=CLWING/(W(15) *W(15) ) 

K'MCRM= (W(4)- (B1 (4) *CLWN+B1 (3) ) )/(Bl (2) *CLWN+B1 (1) ) 
l\'MCI^WMCRNAf(15) 

DM=XM-h’MCR 

Z1=A1(1) 

DO 28 1=2,5 
Z1=Z1+A1 (I) *DM** (I-l) 

28 CONTINUE 
DWME=Z1 

BARL=DESIGN (3) / (GX (5) *10 . ) 

Z1=A2(1) 

DO 30 1=2,4 

Z1=Z1+A2 (I) *BARL** (I-l) 

30 CONTINUE 
FMCR=Z1 
Z1=A3(1) 

DM=XM-FMCR 
BARr-l=0.89+DM 
DO 32 1=2,5 

Z1=Z1+A3 (I) *BARM** (I-l) 

32 CONTINUE 

DFME=(Z1-1.)*CDS(6) 

40 E1=W(6) 

E2=W(7) 

IF(IPHASE.LT.l) E2=0.7 
******************************* 

* WING INDUCED DRAG * 

******************************* 

CLO=W(8) 

Tl=CLO*CLO/ (GX (1 ) *DESIGN (2) *E1) 

T2= (CLWING*CLWING-CLO*CLO) / (GX (1) *DESIGN (2) *E2) 
TT1=T1+T2 

************************************** 
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* MUNK'S INTERFERENCE TERM * 

************************************** 

SMU=HX(4)/W(9) 

G= ( . 25*GX (5) -K3X (34) ) /W (9 ) 

G2=G*G 

G3=G2*G 

C0=. 000076+. 006814*G- . 080417*G2+. 247037*G3 

Cl=l . 002161+2 . 242040*G-34 . 140971*G2+73 . 096667*G3 

C2=.000145-24.824801*G+211.181316*G2-442.515185*G3 

C3=-.014537+42.231817*G-375.564896*G2+803.7551111*G3 

C4=. 009817-24 . 947988*G+220 . 010784*G2-472 . 608148*G3 

SIGMA=C0+C1*SMU4C2*SMU*SMU+C3*SMU*SMU*SMU4C4*SMU*SMU*SMU*SMU 

***************************** 

* INTERFERENCE DRAG * 

***************************** 

T3=2.*SIGMA*CLWING*CLTAIL*DESIGN (4) /(SMU*GX (1) *DESIGN (2) 

$ *DESIGN(1)*E1) 

***************************** 

* TAIL INDUCED DRA3 * 

***************************** 

T4=DESIGN (4) *CLTAIL*CLTAIL/ (DESIGN (1) *GX (1) *DESIGN (5) * .8) 
***************************************** 

* INDUCED DRAG FROM HAVING TAIL * 
***************************************** 

DELCL2=CLWING*CLWING-CL*CL 
TDG=DELCL2/(GX (1) *DESIGN (2) *E2)+T3+T4 
CDI=TT1+T3+T4 

CDDEL= (ABS (CRDE) ) *DERIVCR (14) 

IF (IPHASE.lt. 1) CDDEL=(ABS(CRDE)+GX(2)*GX(13) ) *DERIVAP(14) 

CDTOT=CDI+CDO+CDDEL+DWME+DFME 

XTl=CLO*CLO* (1 ./El-1 ./E2) / (GX (1) *DESIGN (2) ) 

CDOP=CDCH-TDG 

********************************* 

* DATA BASE ASSIGNMENTS * 
********************************* 

GX(31)=TDG 

XLDMX=.5* ( (GX (1) *DESIGN (2) *E2) /(CDOP+XTl) ) **0 . 5 

CLLDMX= (GX (1) *DESIGN (2) *W(7) * (CDOP+XTl) ) **0 .5 , 

IF(IPHASE.LT.l) GO TO 48 

HX(14)=CRDE 

CST(6)=XLDMX 

GX(28)=CLWING 

DTAIL=CDS(3)+TDG 

GO TO 50 

48 CDSAP(4)=CL/CDTOT 
DTAIL=CDSAP (3) +TDG 
CDSAP(5)=CDTOT 
HX(11)=CLTAIL 
GX(29)=CLWING 
HX(12)=CRIT 
HX(13)=CRDE 
50 LOD=CL/CDTOT 
************************** 



* OUTPUT SECTION * 

************************** 

IF(OUTPUT.EQ.O) GO TO 100 
WRITE (6, 80) 

80 FORMAT (*1*//20X*ACCURATE LA) ANALYSIS*) 

IF(IPHASE.LT.O) WRITE (6, 81) 

81 FORMAT (//20X*APPROACH WITH 45 DECREES FLAP*) 

IF(IPHASE.LT.“1) WRITE (6, 87) 

87 FORMAT (20X*FORWARD C.G.*/) 

WRITE (6,82) CL,CLWING,CLTAIL,CRDE,CRIT 

82 FORMAT (//10X*CL (REQUESTED) =*T45 , F15 . 3/10X*CL (WING) =*T45 , F15 . 3/lOX 
$*CL (TAIL) =*T45 , F15 . 3/10X*ELEVATOR (DEGREES) =*T45 ,F15 . 2/ 
$10X*STABILIZER (DEGREES) *T45 , F15 . 2/) 

WRITE (6, 83) TT1,T3,T4,TDG,DTAIL,SIGMA 

83 FORMAT (/5X*INDUCED DRAG COMPONENTS*/10X*WING=*T45,F15.4/ 
$10X*INTERFERENCE=*T45 ,F15 . 4/10X*TAIL=*T45 ,F15 .4 

$ /lOX* (TRIM)=*T45,F15.4/10X*(TAIL)=*T45,F15.4/10X*SIGMA*T45, 

$ F15.4/) 

WRITE(6,84) CDI,CDO,CDDEL,DWT-1E,WMCR,DFME,FMCR,CDTOT 

84 FORMAT (5X*DRAG C0EFFICIENTS*/10X*INDUCED=*T45,F15.4/ 

$ 10X*ZER0 LIFT=* 

$T45 , F15 . 4/10X*ELEVATOR=*T45 , F15 .4/lOX 
$ *V'n:NG (MACH) , M(CRIT)=*T45,F15.4,F15.2/10X 
$ *FUSE(MACH), M(CRIT)=*T45,F15.4,F15.2/10X 
$ *TOTAL=*T45,F15.4/) 

WRITE (6, 86) LOD 

86 FORMAT (//15X*LA>=*T45,F15. 3) 

WRITE (6, 92) XLDMX,CLLDMX,OX)P 

92 FORt^T(///15X*MAX. LA>=*T45,F15.3/15X*CL-LA) f4AX=*T45,F15,3/ 

$ 10X*MODIFIED CDO=*T45,F15.4/) 

100 RETURN 
END 

SUBROUTINE AIRCOST (COST, OUTPUT) 

REAL CD (10) ,CP(10) 

REAL T1 (2) ,T3 (2) ,T9 (2) ,T10 (2) ,T32 (2) 

INTEGER OUTPUT 

COMMON A)EVARA)ESIGN(15) ,ITERM(10) ,CST(10) 

CCM-ION An'SVE/WrS(20) 

COMMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

A=14TS(4) 

DATA T1,T3,T9,T10/1., 1.15,1., 2., 1., 1.15,0.,!./ 

DATA T32/1., 1.1/ 

***************************************************** 

* LIST OF PURCHASE PRICE ELEMENTS IN 1974 $ * 

* REF (4) — NICOLAI * 

* ENGINEERING (1) * 

* DEVELOPMENT SUPPORT (2) * 

* FLIGHT TEST (3) * 

* TOOLING (4) * 

* MANUFACTURING LABOR (5) * 

* QUALITY CONTROL (6) * 

* MATERIALS (7) * 

* ENGINE (8) * 
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UCJO UCJU UOCJ 


* AVIONICS (9) * 

* ACTIVE CONTROL SYSTEM (10) * 

***************************************************** 

IACT=0 

IOGX=ITERM(2) 

CD (1) =T1 (IACT+1 ) *10964 . 13*A**0 . 791 
CP (1) =T1 (IACT+1 ) *26567 . 12*A**0 .,791-CD (1) 

CD (2) =1627 . 68*A**0 .873 
CP(2)=0.0 

CD (3) =T3 (IACT+1 ) *T32 (IGGX+1 ) *19 . 16*A**1 . 160 
CP (3) =T3 (IACT+1 ) *T32 (ICGX+1 ) *0 . 0 
CD (4) =15716 . 31*A**0 . 764 
CP (4)=43272 .95*A**0 .764-CD (4) 

CD (5) =13026 . 56*A**0 . 74 
CP (5) =164219,. 18*A**0 .74-CD (5) 

CD(6)=CD(5)*0.13 
CP(6)=CP(5)*0.13 
CD(7)=2733.64*A**0.689 
CP (7) =125960 . 83*A**0 . 689-CD (7 ) 

CD (8 ) = (ITERM (4 ) +1 ) *2 . *1 . 31*169 . 0* (DESIGN (6) /ITERM (4) ) **0 . 8356 

CP (8)=ITERM(4) *250. *1.31*169.0* (DESIGN (6)/ITERT4(4))**0. 8356 
CD (9) =2 . *300000 . *T9 ( IACT+1 ) 

CP (9 ) =250 . *300000 . *T9 ( IACT+1 ) 

ACTIVE CONTROLS PRICE HAS TO BE ESTIMATED 

CD(10)=206250.*2.*ITERM(1) 

CP (10) =206250 . *250 . *ITERM (1) 

TOTD=0. 

TOTP=0. 

CONVERT FROM 1974 $ TO 1976 $ 

DO 50 J=l,10 
CD (J)=CD(J) *1.23077 
CP (J)=CP(J) *1.23077 
TOTD=TOTD4GD(J) 

T0TP=T0TP4GP (J) 

50 CONTINUE 

INCLUDING 10% PERCENT PROFIT 

TOTCOST=TOTD*l . l+TOTP*l . 1 
COST=TOTCOST/250 . 0 
COST^COST* (1 . 0+PX (7 ) *0 . 05) 

************************** 

* OUTPUT SECTION * 

************************** 

IF (OUTPUT. EQ.O) GO TO 100 
WRITE (6,70) 

70 FORMAT (*1*//20X*AIRCRAFT COST ESTIMATES*//43X*DEVELOPMENT* 
$9X*PRODUCTION*/) 

WRITE(6,72) CD(1),CP(1) 
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72 FORMAT (10X*ENGINEERING*T40 , F15 . 2 ,T60 , F15 . 2) 

WRITE(6,74) CD(2),CP(2) 

74 FORMAT (10X*DEVELOPMEMT SUPPORT*T40,F15.2,T60,F15.2) 
Ia/RITE(6,80) CD(3),CP(3) 

80 FORf-lAT(10X*FLIGHT TEST*T40,F15.2,T60,F15.2) 

VJRITE(6,82) CD(4),CP(4) 

82 FORMAT (10X*TOOLING*T40 ,F15 o 2 ,T60 , F15 , 2) 
l'JRITE(6,84) CD(5),C?(5) 

84 FORTJlAT(10X*r'1ANUFAC. LAEOR*T40,F15.2,T60,F15.2) 
WRITE(6,86) CD(6),CP(6) 

86 FORMAT (10X*QUALI1Y C0MTR0L*T40,F15,2,T60,F15o2) 
IVRITE(6,88) CD(7),CP(7) 

88 FORT4AT (10X*MATERIALS*T40 , F15 . 2 ,T60 , F15 , 2) 
liJRITECe^gO) CD(8),CP(8) 

90 FORMAT (10X*ENGINE*T40 ,F15 « 2 ,T60 , F15 . 2) 

WRITE(6,92) CD(9),CP(9) 

92 FORTJIAT (10X*AVIOMICS*T40 ,F15 . 2 ,T60 , F15 . 2) 

VJRITE (6,94) CD (10) ,CP (10) 

94 FORMAT (10X*ACTIVE CONTROLS SYSTEM*T40,F15.2,T60,F15.2) 
VJRITE(6,96) TOTD,TOTP 
96 FORimT (/10X*TOTAL*T40 , F15 , 2 ,T60 , F15 , 2) 

VJRITE(6,78) COST 

78 F0RT-1AT(/5X*T0TAL COST PER AIRCRAFT= $*F12,2) 

100 RETURN 
END 

SUBROUTINE MAINCST(COST,OUTPUT) 

INTEGER OUTPUT 

REAL MCOST(27) ,LCOST(27) ,T9(2) ,XNM(27) 

REAL Tl(2) ,T3(2) ,T3(2) ,T12(2) ,T15(2) 

REAL LCST(27) ,MCST(27) 

COMiMON A>EVARA>ESIGM (15) ,ITERM(10) ,CST(10) 

C0MT40N /WTSVEAJTS(20) 

COMMON /GEQM/l'J(20) ,HX(20) ,GX(35) ,PX(15) 

DATA Tl,T3,T8/lc,lol5,l.,lo2,l,,lo2/ 

DATA T9 ,T12,T15/1.,1. 3, lo, 1.15,1,, 1,15/ 

IACT=ITERM(1) 

ICGX=ITERM(2) 

NENG=ITEra-l(4) 

NPASS=GX(19) 

VJTCWflITS (1) *0.453592 
WTE=IVTS (2) *0.453592 
VJTF=’.irrS (6) *0.456592 

DATA XNM/4HINSP,0RAIR COND,10HAUTO PILOT, 6HCOMMUN, 
$4HELEC,4HFURN,9HFIRE PROT,9HFLT CONTL,4HFUEL, 

$9HHYD POWER,3HICE,5HINSTR,9HLAND GEAR, 

$8HLIGHTING , 5HNAVIG , 6KOXYGEN , 7HPNUEMAT , 
$9HWAT/WASTE,7HAIR APU,9HSTRUCTURE,5HDOORS, 

$8HFUSELAGE , 8HNACELLES , 5HWINGS , 4HSTAB , 

$7HtrJINDOl,'JS, 6HENGINE/ 

* MAINTENANCE COSTS — 1976 DOLLARS/HOUR * 

* REF (20) — AMERICAN AIRLINES * 

* , (1) INSPECTION AND MISC. * 
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* 

(2) 

AIR CONDITIONING 

* 

* 

(3) 

AUTO PILOT 

* 

* 

(4) 

COMMUNICATIONS 

* 

* 

(5) 

ELECTRICAL 

* 

* 

(6) 

EQUIPMENT AND FURNISHINGS 

* 

* 

(7) 

FIRE PROTECTION 

* 

* 

(8) 

FLIGHT CONTROLS 

* 

* 

(9) 

FUEL 

* 

* 

(10) 

HYDRAULIC POWER 

* 

* 

(11) 

ICE AND RAIN 

* 

* 

(12) 

INSTRUMENTS 

* 

* 

(13) 

LANDING GEAR 

* 

* 

(14) 

LIGHTING 

* 

* 

(15) 

NAVIGATION 

* 

* 

(16) 

OXYGEN 

* 

* 

(17) 

PNUEMATICS 

* 

* 

(18) 

WATER/WASTE 

* 

* 

(19) 

AIRBORNE APU 

* 

* 

(20) 

STRUCTURE 

* 

* 

(21) 

DOORS 

* 

* 

(22) 

FUSELAGE 

* 

* 

(23) 

NACELLES/PYLONS 

* 

* 

(24) 

VJINGS 

* 

* 

(25) 

STABILIZERS 

* 

* 

(26) 

WINDOWS 

* 


************************************************ 


LC0ST(1)=T1 (IACr+1) *7.66+0.377*OTE/1000.0 
MC0ST(1)=T1 (IACT+1) *1.21+0. 062*V-«JTE/1000.0 
DATA LCOST(2) ,MCOST(2)/5.1026,4.52/ 

LCOST ( 3 ) =T3 ( IACT+1 ) *11 . 19 

KCOST(3)=T3 (IACT+1) *2.621 

LCOST (4 ) = . 0276*NPASS 

MCOST (4) =0 , 0118*NPi^S 

DATA LCOST(5) ,MCOST(5)/4. 306, 5.748/ 

LCOST (6 ) =9 . 11+0 . 08496*NPASS 
MCOST (6 ) =2 . 38+0 . 05776*NPASS 
LCOST (7 ) = . 213+2 . 29* (2 .+NENG) 

MCOST (7 ) =0 . 365* (2 .+NENG) 

LCOST(8)=T8 (IACT+1) *6.84+0.0035*WTO/1000. 

MCOST (8)=T8 (IACT+1) *3. 876+0. 00655*WTO/1000. 

LCOST (9) =1 . 114+0 . 0262*WTF*T9 (ICGX+1 ) /lOOO . 0 

MCOST (9) =0 . 595+0 . 0123*WTF*T9 (IOGX+1 ) /lOOO . 0 

DATA LCOST(IO) ,MCOST(10) /3. 33, 3,95/ 

LCOST(ll) =. 5089+0 . 0013*WTO/1000 .0 

MCOST (11) = . 0847+. 0037*WTO/1000 . 0 

LCOST (12) =T12 (IACT+1 ) *0 . 509+. 009*WTE/1000 , 0 

MCOST(12)=T12 (IACT+1) *0.235+.0031*WTE/1000.0 

LCOST (13) =4 . 58+. 071*WTO/1000 . 0 

MCOST (13)=4 .961+. 181*WTO/1000. 0 

LCOST (14) =1 . 51+0 . 01152*NPASS 

MCOST (14) =0 . 047+0 . 01392*NPASS 

LCOST(15)=T15 (IACT+1) *10.077 

MCOST (15) =T15 (IACT+1 ) *7 . 166 
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LCOST (16) =. 515+0 . 00265*NPASS 
MCOST (16) =0 . 00752*NPASS 
DATA AC/200./ 

T=DESIGN (6) *4 .448/13ENG 

LCOST (17) =0 . 181+. 0003*AC*T/10000 . 

MCOST (17) =0 . 0019*AC*T/10000 . 

LCOST (18 ) = . 339+0 . 00368*NPASS 

MCOST (18) =0 . 00768*NPASS 

DATA LCOST(19) ,MCOST(19) /.315, .462/ 

LCOST(20)=3.+.0099*WTE/1000. 

DATA MCOST ( 20 )/0./ 

LCOST (21) =1 . 147+0 . 006*NPASS 
MCOST (21) = . 387+0 . 00785*NPASS 
LCOST (22) =1 . 5+0 . 046*WTE/1000 . 

DATA MCOST (22) /O. 5833/ 

NAC=4 

IF(NENG.LT.4) NAC=2 
LCOST (23 ) = . 3366'mC 
MCOST (23) =.1391*NAC 
SW=DESIGN (1 ) / (3 . 201*3 . 281 ) 

DATA LCOST (24) /2. 9475/ 

MCOST (24) =0 . 126+0 . 00506*SW 

DATA LCOST(25) ,MCOST(25)/0. 8321, 0.3737/ 

LCOST (26) =0 .763+0 .00043*NPASS 
MCOST (26) =0 .0362*NPASS 
DO 50 K=l,26 
LCST(K)=LCOST(K)/2.5 
MOST (K) =MCOST (K) /2 . 5 
50 CONTINUE 
TLCOST=0. 

•nviCO3T=0. 

DO 75 K=l,26 
TLCOST=LCST (K ) +TLCOST 
TMCOST=MCST (K) +TMCOST 
75 CONTINUE 

TENG=DESIGN (6) /ITERM (4) 

*************************************************************************** 

* MAIN ENGINE COST * 

*************************************************************************** 

LCOST(27)=(ITERM(4)/(4.0*2.5))*88.5*(TENG/20000.)**0.5 * 

MCOST ( 27 ) = ( ITERM (4 ) / ( 4 . 0*2 , 5) ) *109 . 0* (TENG/20000 . ) **0 . 5 * 

COST^TLCOST+TMCOST+LCOST (27) +MCOST (27) * 

COS'D=COST*(1.0+PX(6)*0.05) * 

*************************************************************************** 
IF (OUTPUT. EQ.O) GO TO 100 
WRITE (6,86) 

86 FORMAT (*1*//30X*MAINTENANCE OPERATING COSTS*//7X*NO. SYSTEM*7X 
$10X*LABOR*7X*MATERIAL* ) 

DO 90 K=l,26 

WRITE(6,88) K,XNM(K),LCOST(K),MOOST(K) 

88 FORMAT(I10,2X,A10,3X,2F15.2) 

90 CONTINUE 

WRITE(6,92)TLCOST,TMCOST,LCOST(27) ,MCOST(27) ,COST 
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92 TORMAT(//5X*LAB0R COST*T35,F15.2/5X*MATERIAL COST*T35,F15.2//5X 
$*ENGINE LABOR COST*T35,F15.2/5X*ENGINE f^IATERIAL COST*T35,F15.2 
$///30X*MAINTOENACE DOC IN 1976 DOLLARS PER HOUR*F15.2) 

100 RETURN 
END 

SUBROUTINE CNSTRN (OUTPUT) 

REAL LODMA,LOD2,ANS(4) ,LFL,NZOACR,NZOAA 
REAL PRAMCR(4) ,PRAMAP(4) ,ACR(5) ,AAP(5) 

REAL TCCR(4) ,TCAP(4) 

CO:4PLEX ROOTCR(4) ,ROOTAP(4) 

COMMON /STRAIN/CON (59) 

INTEGER OUTPUT 

REAL KT1,KT2,KTD0T1,KTD0T2,DERCR(15) , DERAP (15) 

COMMON /CONSTR/SU(59) ,SL(59) ,XINEQ(59) 

C0r4M0N /DEVARA>ESIGN (15) ,ITER?4(10) ,CST(10) 

CO:jIMON /DRAG/CDS (6), CDSAP (6) 

COMiMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMuMON /GRAVITY/GG(6) 

COMMON /STAB/DERIVCR(15) ,DERIVAP(15) ,STOR(20) 

COMMON /WTSVE/WTS(20) 

*************************************************************************** 


* CONSTRAINT IDENTIFICATION 


* 

* NO. DESCRIPTION 

* ******* **************************************************** 

* 1 CRUISE THRUST REQUIREMENT 

* 2 SECOND SEGMENT CLIMB GRADIENT - THRUST REQUIREf^lENT 

* 3 MISSED APPROACH CLIMB GPJ^IENT - THRUST REQUIREMENT 

* 4 LANDING FIELD LENGTH - VJIKG LOADING REQUIRET4ENT 

* 5 .... TAKE-OFF FIELD LENGTH - VUNG LOADING REQUIREMENT 

* 6 LANDING GEAR - AFT CG LIMIT 

* 7,8 STATIC STABILITY - CRUISE, APPROACH 

* 9,10 .MANUEVER MARGIN - CRUISE, APPROACH 

* 11 TAIL LIFT - APPROACH 

* 12 NOSE GEAR UNSTIQ< 

* 13,14 ...DYNAMIC STABILITY - CRUISE, APPROACH 

* 15,16 PHUGOID FREQUENCY - CRUISE, APPROACH 

* 17,18 PHUGOID DAT-1PIKG - CRUISE, APPROACH 

* 19,20 SHORT PERIOD FREQUENCY - CRUISE, APPROACH 

* 21,22 SHORT PERIOD DAMPING - CRUISE, APPROACH 

* 23,24 TIME-TO-DOUBLE (CRUISE, APPROACH) 

* 25,26 TIME-TO-HALF (CRUISE, APPROACH) 

* 27 FLIGHT PATH STABILITY - APPROACH 

* 28,29 VERTICAL GAIN - CRUISE, APPROACH 

* 30,31 TIME SUB THETA 2 - CRUISE, APPROACH 

* 32,33 FREQUENCY**2/VERTICAL GAIN - CRUISE, APPROACH 

* 34 T(l) PARAMETER - APPROACH 

* 35,36 RATIO OF MODE FREQUENCIES - CRUISE, APPROACH 

* 37-39 ELEVATOR VARIANCE - CRUISE 

* 40-42 ELEVATOR VARIANCE - APPROACH 

* 43,44 VARIANCE OF ELEVATOR RATE - CRUISE, APPROACH 

* 45 PASSENGER VOLUME LIMIT 

* 46,47 ELEVATOR DEFLECTIONS - TRIMMED (CRUISE, APPROACH) 


* 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 
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* 48 CRUISE ALTITUDE 

* 49 CRUISE ALTITUDE ((L/D)MAX) 

* 50,51 WING CL - CRUISE, APPROACH 

* 52 TAIL ASPECT RATIO LIMIT - AR (TAIL/AR (WING) 


*************************************************************************** 
IOUT=OUTPUT 
IF (OUTPUT, GT.l) IOUT=0 
CON(1)=GX(30) 

******************************************************************* 

* FIND TAKE-OFF CL-MAX AND SECO!'® SEGMENT CLIMB GRADIENTS * 
******************************************************************* 

CLVJTO=GX(8) 

CALL XLOD(CLWTO,.1,0.,EFF,-4,0) 

NENG=ITERM(4) 

CALL AT62(0.,ANS) 

TOWAV=DESIGN (6) AJTS (1) 

CL2=CLWTO/l,44 

V=(2.*WTS (1) /(DESIGN (1) *ANS (1) *CL2) ) **0.5 
XM=V/ANS(4) 

CALL ENGINE (0 . ,XM,TCn4,TSFC) 

GRAD=0,030 

CALL XLOD(CL2,. 1,0., LOD2, 0,0) 

IF(ITERM(4) .EQ.3) GRAD=0.027 
IF(ITERM(4).EQ.2) GPvAD=0.024 

TOWRQ2= (NENG/ (NENG-1 .))*(! . AOD2+SIN (GRAD) ) * (1 ./TCTM) 

CON (2 ) =TOVAV/TOWRQ2 
***************************** 

* NOSE GEAR UNSTICK * 

***************************** 

ZACT=4.0/W(10) 

XMUF=0,025 

ZACLG=2.*ZACT 

XLT=HX (6) *DESIGN (1) *W(10) /DESIGN (4) 

CALL AT62(0.,ANS) 

VSTALL= (VJTS (1) *2.0/ (DESIGN (1) *ANS(1) *GX(8) ) ) **0.5 
VLO=VSTALL*0.9 

CLW= (STOR (6) * (W(2)-GX (17) -3,0)) /GX (2) 

Q=0 , 5*ANS ( 1 ) *VLO*VLO 
XLW=Q*DESIGN (1) *CLW*1 .2 
X^7/J=1.05*Q*DESIGN (1) * (W(19)+W(17) ) *W(10) 

C1=XLT+ (0 , 25-OG (5) -XMUF*ZACLG) *W(10) 

UNTHRU=DESIGN (6) 

IF(ITERM(4) .EQ.3) in^]THRU=DESIGN(6)/3,0 
X03=CG(3) 

COIT=HX(12) 

C 

C TAIL INCIDENCE SET TO WORSE CASE FOR TAKE-OFF (BOEING) 

C 

IF(GX(16) .NE.-99.0) COIT=GX(16) 

XLTRQ= (XMl*R-ZACT*W(10) * .9*UNTnmJ-K4(10) * (XCG-.25) *XLW 
$ -V/(i0) * (CG (5 ) -CG (3) •iZACLG*}it'rJF) * (VJTS (1) -XLW) ) /Cl 
DET-mX=20, 

XLTAV=STOR (7 ) *HX (6 ) * , 9* ( (1 ,-5T0R (8 ) ) * (W (2) -GX (17) -3 . 0) 



$-W(2)+C0IT-HX (3) *DEMAX) *Q*DESIGN (4) /GX(2) 

C 

C GROUND EFFECTS — DATCOM FUNCTION OF GEOMETRIC ALPHA 
C 

TLGE= . 355*Q*DESIGN (4 ) 

XLTAV2=XLTAV+TLGE 
CLTAV=XLTAV2/ (Q*DESIGN (4 ) ) 

IF (CLTAV.GT. 1 .5) XLTAV2=-1 . 5*Q*DESIGN (4) 

CON (12) =XLTAV2 ALTRQ 

IF(OUTPUT.EQ.l) WRITE(6,548) ZACT,XMUF,ZACLG,ANS,VSTALL,CLW 
$ ,Q,XLW,XMW,Cl,UNTHRU,XOG,XLTRQ,CL2,COIT,XLT 
$ ,XLTAV,CON(12) ,STOR(7) ,STOR(8) ,W(2) ,GX(17) ,HX(6) ,W(10) 

$ ,XLTAV2,TLGE,CLTAV 

548 FORMAT (*1*//20X*DEBUG OF NOSE GEAR UNSTICK*/// 

$10X"ZACT,XMUF, ZACLG="3F15 .4/ 

$ 10X"ANS(1,2,3,4)=" 4F15.6/ 

$ 10X"VSTALL ,CLW,Q,XLW,XMl-f="5F15 . 4/ 

$ 10X"C1,UNTHRU,XCG,XLTRQ,CL2,COIT,XLT="7F12.4/ 

$ 10X"XLTAV,CON(12) ,STOR(7) ,STOR(8) ,W(2)=”5F15.4/ 

$ 10X"GX(17) ,HX(6) ,W(10) ,XLTAV2,TLGE,CLTAV="6F15.4/) 
******************************************************************** 

* FIND APPROACH CL-MAX AND MISSED APPROACH CLIMB GRADIENTS * 

* SET UP APPROACH WITH FORWARD CG AND 45 DEG FLAPS * 

******************************************************************** 

CLWMAX=GX(12) 

CALL XLOD(CLWMAX,.1,0.,EFF,-3,0) 

CLS=CLWMAX 

GRAD=GRAD-.003 

CLA=CLS/1.69 

V= (2.*WTS (1) *GX (27) / (DESIGN (1) *ANS (1) *CLA) ) **0 . 5 
XM=V/ANS(4) 

CALL ENGINE (0.,XM,TCTM,TSFC) 

GX(10)=CLA 

CALL XLOD(CLA,.1,0.,LODMA,-2,0) 

TQWRQ3= (NENG/ (NENG-1 . ) ) * d . AODMA+SIN (GRAD) ) * (1 ,/TCIM) 
TOWAV3=DESIGN (6) / (WPS (1) *GX (27) ) 

CON ( 3 ) ==TOWAV3/TWRQ3 

VA2=GX (27) *WTS (1 ) *498 . 23/ (DESIGN (1) *CLA) 

VA=VA2**0.5 

GX(14)=VA*1.6881 

LFL=0.29875*VA2+25. 

CON(4)=LFL 

TOWAV=DESIGN (6 ) /WTS (1 ) 

TOP= (WTS (1) /DESIGN (1) ) /(CLWTO*TOWAV) 

TOFL= (31. 7*TOP) +910.0 
CON(5)=TOFL 

CON(6)=OG(5)-CG(6)-DESIGN(7) 

VOLPAS=GX(19)*52. 

PASSL=DESIGN (3) -1 . 2*GX (15) - . 5*GX (33) 

VOLAV=GX (1) *GX (5) *GX (5) *PASSL/8 . 

CON (45) =V0LAVA0LPAS 
CLTAIL=HX(11) 

********************************* 
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* DATA BASE ASSIGNMENTS * 


C0N(11)=CLTAIL 

CON(46)=HX(14) 

C0N(47)=HX(13) 

CON(48)=WTS(ll) 

C0N(49)=WTS(12) 

CON(50)=GX(28) 

C0N(51)=GX(29) 

CON (52) =DESIGN (5) /DESIGN (2) 

*************************************************************************** 

* THIS ENDS DESIGN CONSTRAINT SECTION, AND BEGINS HANDLING QUALITY * 

* CONSTRAINT SECTION * 

*************************************************************************** 

CALL XLOD(CLA, .1,500. rEFF,-l,IOUT) 

CALL XLOD(GX(9) ,GX(3) ,WTS(11) ,EFF,2,0) 

CON(7)=CG(l)-STOR(4) 

CON(8)=OG(2)-STOR(9) 

CON(9)=CG(l)-STOR(5) 

CON (10) =OG (2) -STOR (10) 

********************************************** 

* CALCULATE DIMENSIONAL DERIVATIVES * 

********************************************** 


CALL AT62 (WTS (11) ,ANS) 

UOCR=GX(3)*ANS(4) 

CALL DIMDER(UOCR,l,DERCR) 

CALL DIt€)ER(GX(14),0,DERAP) 

IF(OUTPUT.EQ.l) VJRITB(6,20) DERCR, DERAP 
28 FORMAT (*1*//30X*/ / / DIMENSIONAL STABILITY DERIVATIVES/ / / * 
$ //10X*DERCR=*/5(T30,3F12.4/)//10X*DERAP=*/5(T30,3F12.4/) ) 

CON (13) =DERCR (4) *DERCR (3) -DERCR (1) *DERCR (6) 

CON (14) =DERAP (4) *DERAP (3) -DERAP (1) *DERAP (6) 


IFULL=1 

IF(IFULL.GT.O) GO TO 30 

************************************************* 


********************* 


* THIS NEXT SECTION 7\PPROXIMATES AIRCRAFT DYNAMIC PROPERTIES * 
********************************************************************** 

TPHUGCR=DERCR (1) *DERCR (9) -UOCR*DERCR (3) 

OMPH2=GX(15) * (DERCR(3) *DERCR(4)-DERCR(6) *DERCR(1) ) /TPHUGCR 

OMPH= (ABS (OMPH2 ) ) **0 . 5 

CON(15)=OMPH 

TPHUGAP=DERAP (1) *DERAP (9) -GX (14) *DERAP (3) 

OMPHA2=GX (15) * (DERAP (3) *DERAP (4) -DERAP (6) *DERAP (1) ) /TPHUGAP 

OMPHA= (ABS (OMPHA2) ) **0 . 5 

CON(16)=OMPHA 

TTT=-DERCR(5)-(DERCR(6)* (DERCR(2)*UOCR-GX(15) ) )/TPHUGCR 
ZETAPH=TTT/ ( 2 . 0*OMPH) 

CON(17)=ZETAPH 

TTTA=-DERAP (5) - (DERAP (6) * (DERAP (2) *GX (14) -GX (15) ) ) /TPHUGAP 
ZETAPHA=TTTA/ ( 2 . 0*OMPHA) 

CON(18)=ZETAPHA 

OMSP2=DERCR (9) *DERCR (1) -DERCR (3) *UOCR 
OMSP= (ABS (OMSP2 ) ) **0 . 5 
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C0N(19)=0MSP 

0MSPA2=DERAP (9) *DERAP (1) -DERAP (3) *GX (14) 

OMSPA= (ABS (0MSPA2) ) **0 « 5 
CON(20)=OMSPA 

TXT=- (DERCR (1 ) +DERCR (9 ) +UCX:R*DERCR (12) ) 

ZETASP=1XT/ (2 . 0*OMSP) 

C0N(21)=ZETASP 

TXTA=- (DERAP (1) +DERAP (9)-K3X (14) *DERAP (12) ) 

ZETASPA=TXTA/ (2 . 0*OMSPA) 

C0N(22)=ZETASPA 

************************************************************* 

* CALCULATE EXACT DYNAMICS FROM FOURTH OPDER MODEL * 

************************************************************* 

30 CALL LONGRT (DERCR, UOCR,ROOTCRrPRAT^CR,ACR,NXR,TCCR) 

CALL LONGRT (DERAP,VA,ROOTAP,PPJmP,AAP,NOAP,TCAP) 
IF(IFULLoLT.l) GO TO 45 

0MPH=PRAI^1CR(1) 

ZETAPH=PRAMCR(2) 

OM3P=PRAT^CR(3) 

ZETP^P=PRAMCR (4 ) 
o:j!?ha=praT'1Ap (i) 

ZETAPHA=PRAr-mP ( 2 ) 

0MSPA=PRAT'1AP (3) 

ZETASPA=PRAMAP (4 ) 

45 IF (OUTPUT, LTol) GO TO 50 
IF(OUTPUT.GT.l) GO TO 50 
l'JRITE(6,31) 

31 FORMAT (*1*///30X*LONGITUDINAL DYNAMICS*//20X*/ / CRUISE/ / *) 
li!RITE(6,32) ACR 

32 FOPuMAT (/10X*COEFFICIEM'S=*5F15 ,6) 
l'JRITE(6,33) ROOTCR 

33 FOR?-LAT (/10X*ROOTS (REAL , IM/^J3INARY) */4 (2F15 , 4/) ) 
l'JRITE(6,35) PRAMCR 

35 FORMAT (/10X*PHUGOID FREQUENCY*T35,F15,4/10X 
$*PHUGOID D.Ar®ING*T35,F15,4/10X*SHORT PER, FREQ,* 
$T35,F15.4/10X*SH0RT PER, D.AM?IKG*T35,F15,4) 

l«JRITE(6,36) N0CR,TCCR 

36 FORfJlAT(/10X*NO, OF NON-OSCILLATORY ROOTS=*Ill 
$ /10X*TIME CONSTANTS=*T35,4F15,4) 

l'JRITE(6,37) 

37 FORriAT(//20X*/ / APPROACH/ / *) 
l*JRITE(6,32) AAP 
VJRITE(6,33) ROOTAP 
l'JRITE(6,35) PRAMAP 
VJRITE(6,36) N0AP,TCAP 

liJRITE(6,48) GX(8) ,CLWTO,CL2,GX(12) ,CLWMAX,CLA 
48 FORMAT (//10X*CL-MAX TO(W) ,CL-MAX TO (AC) ,CL2*T60,3F10,3/ 

$ 10X*CL-MAX (W), CL-MAX (AC), CLA=*T60,3F10,3) 

50 IF(IFULL.LT.l) GO TO 100 
CON(15)=OMPH 
CON(16)=OMPHA 
C0N(17)=ZETAPH 
CON(18)=ZETAPHA 
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C0N(19)=0MSP 

CON(20)=OMSPA 

C0N(21)=ZETASP 

C0N(22)=ZETASPA 

0MSP2=0MSP*0MSP 

0MSPA2=0MSPA*0MSPA 

IF(OMSP.EQ.O.) OMSP=1.0E-20 

IF(OMSPA.EQ.O.) OMSPA=1.0E-20 

IF(OMPH.EQ.O.) OMPH=1.0E-20 

IF(OMPHA.EQ.O.) OMPHA=1.0E-20 

GREAT=-999. 

100 IF(NOCR.LT.l) GO TO 210 
DO 209 K=1,N0CR 
IF(TCCR(K).LT.0.0) GO TO 209 
IF(TCCR(K) .GT.GREAT) GREAT=TCCR(K) 

209 CONTINUE 

210 TDOU3CR=-0.693/(OMPH*ZETAPH) 

Tl=-0 .693/ (OMSP*ZETASP) 

T2=0. 693/GREAT 

IF (TDOUBCR.GT.Tl .AND.Tl.GT.O .0) TD0UBCR=T1 
IF(T2.LT.0,0) GO TO 219 

IF (TDOUBCR.LT.O.O.OR.TDOUBCR.GT.T2) TDOUBCR=T2 
219 IF(TDOUBCR.LT.O.O) TDOUECF^99. 

CON(23)=TDOUBCR 

GREAT=-999. 

IF(NOAP.LT.l) GO TO 230 
DO 231 K=l,NOAP 
IF(TCAP{K) .LT.0.0) GO TO 231 
IF(TCAP(K) .GT.GREAT) GREAT=TCAP(I<) 

231 CONTINUE 

230 TDOUBAP=-0.693/(OMPHA*ZETAPHA) 

Tl=-0 . 693/ (OMSPA*ZETASPA) 

T2=0. 693/GREAT 

IF (TDOUBAP .GT.Tl , AND.Tl .GT. 0.0) TD0UBAP=T1 
IF(T2.LT.0.0) GO TO 239 

IF (TDOUBAP . LT. 0 , 0 .OR.TDOUBAP .GT.T2) TDOUBAP=T2 
239 IF (TDOUBAP. LE.O.) TDOUEAP=99. 

CON(24)=TDOUBAP 

THALF=-99. 

THALFA=-99. 

IF(OMSP.EQ.1.0E-20) OMSP=-1.0E-20 
IF(ONSPA.EQ.1.0E-20) OMSPA=-1.0E-20 
IF (Or-IPH. EQ . 1 . OE-20) OMPH=-l . OE-20 
IF (OMPHA.EQ. 1 .OE-20) 0:4PHA=-1 .OE-20 
GREAT=-999. 

IF(NOCR.LT.l) GO TO 245 
IF(TD0UBCR.NE.99.) GO TO 250 
DO 242 K=l,NOCR 
IF(TCCR(K).GT.0.0) GO TO 242 
IF (TCCR(K) .GT.GREAT) GREAT=TCCR(K) 

242 CONTINUE 

245 THALF=0.693/(OMPH*ZETAPH) 

T1=0 . 693/ (OMSP*ZETASP) 
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T2=-0. 693/GREAT 

IF(THALF.LT.Tl.AND.Tl.GT.O.O) THALP=T1 
IF(THALF.LT.T2.AND.T1.GT.0.0) THALF=T2 
IF(THALF.LE.O.) THALF=-99. 

250 CON(25)=THALF 
GREAT=-999. 

IF(NOAP.LT.l) GO TO 255 
IF(TDOUBAP.NE.99.) GO TO 260 
DO 252 K=l,NOAP 
IF(TCAP(K) .GT.0.0) GO TO 252 
IF(TCAP(K) .GT.GREAT) GREAT=TCAP (K) 

252 CONTINUE 

255 THALFA=0.693/(OMPHA*ZETAPHA) 

T1=0 . 693/ (OMSPA*ZETASPA) 

T2=-0. 693/GREAT 

IF(THALFA.LT.Tl.AND.Tl.GT.O) THALFA=T1 
IF(THALFA.LT.T2.AND.T2.GT.O) THALPA=T2 
IF(THALFA.LE.O.) THALFA=-99. 

260 CON(26)=THALFA 

TDEL=DERAP (13) /DERAP (15) 

TZT= (DERAP (6 ) *DERAP (4 ) -DERAP (6 ) *DERAP ( 1 ) ) / 

$ (-DERAP(1)+DERAP(3)*TDEL) 

TYT= (DERAP (4 ) -TDEL*DERAP (6) ) /(DERAP (1) -DERAP (3) *TDEL) 

DGDU= (DERAP (5) - (DERAP (2) -GX (15) /GX (14) ) *TYT-DERAP (14) *TZT 
$/DERAP(15))/GX(15) 

CON(27)=DGDU 

ZZT=UOCR* (DERCR(13) *DERCR(3)-DERCR(15) *DERCR(1) ) 

NZOACR=ZZT/ ( (DERCR (15) -DERCR (13) *DERCR (9) /UOCR) *GX (15) ) 
CON(28)=NZOACR 

ZZTA=GX (14) * (DERAP (13) *DERAP (3)-DERAP (15) *DERAP (1) ) 

NZOAA=ZZTA/ ( (DERAP (15) -DERAP (13) *DERAP (9) /GX (14) ) *GX (15) ) 
CON(29)=NZOAA 

TTH2=ZZT/ (UOCR* (DERCR (15) +DERCR (13) *DERCR (12) ) ) 

CON(30)=TTH2 

TTHA2=ZZTA/ (GX (14) * (DERAP (15)+DERAP (13) *DERAP (9) ) ) 

CON(31)=TTHA2 

CON (32) =OMSP2/NZOACR 

CON (33) =OMSPA2/NZOAA 

PIARE=1 ./ (GX (1) *DESIGN (2) *W(7) ) 

TA1=GX (14) / (2 .0*GX (15) * (1 ./CDSAP (4) -2 .*PIARE*GX (10) ) ) 

C0N(34)=TA1 

IF(OMPH*OMPH.GT.O.) GO TO 80 
CON (35) =99. 

GO TO 82 

80 CON(35)=OMSP/CWPH 
82 IF(OMPHA*OMPHA.GT.O.) GO TO 84 
CON (36) =99. 

GO TO 86 

84 CON(36)=OMSPA/OMPHA 

************************************************************************ 

* THE FOLLOWING SECTION CALCULATES THE RESPONSE OF A PITCH * 

* ATTITUDE HOLD/RATE COMMAND AUTOPILOT IN TURBULENCE. * 

************************************************************************ 
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86 WCR2=NZ0ACR 
WCR=WCR2**0.5 
DZETA=.7 

DM=DERCR (15) +DERCR (3) 

DM1=1 ,-DERCR (12) *DERCR (12) *DERCR (15) *WCR2A>M 

DM2=2 .*DZETA*WCR*DERCR (12) *DERCR(15) -DSRCR (15) -DERCR (3) *DERCR (15) 
KTD0T1=2 . *DZETA*WCR* (DERCR (12) *U0CR-1 . ) -DERCR (9)-DERCR (3) *UOCR 
$ +DERCR (12) *WCR2* (1 .-DERCR (12) *UOCR) /DM 
KTD0T1=KTD0T1/ (DM1*DM2) 

KT1=WCR2* (l.-DERCR(12)*U0CR)+KTD0Tl*VTCR2*DERCR(12)*DERCR(15) 
KT1=KT1/DM 

XNUM= ( (DERCR(9) /UOCR) **2,0) * (KT1*KT1-KTDOT1*KTDOT1* (2.0*DZETA 
$*WCRWCR*WCR*GX (11) / (UOCR) ) ) 

TSTO=V:CR2*GX (11) /UOCR 

DEN0M=TST0-(1.+2,*DZETA*TST0A’CR) * (2,*DZETA*WCR+TSTO) 

SIG= (ABS (XNIM/DENOM) ) **0 . 5 
C0M(37)=3,0*SIG 
HX(15)=3.0*SIG 
C 

C ASSUME THAT CRUISE RMS IS 3 FT/SEC 

C 

C0N(38)=KTD0T1 
C0J3(39)=KT1 
WA2=NZ0AA . 

WA=Vm2**0.5 
DM=DERAP (15) +DERAP (3) 

DM1=1 o-DERAP (12) *DERAP (12) *DERAP (15) *WA2/DM 

DM2=2 o*DZETA*WA*DERAP (12) *DERAP (15) -DERAP (15) -DERAP (3) *DERAP (15) 
KTD0T2=2o*DZETA*liJA* (DERAP (12) *GX (14) -1 , ) -DERAP (9) -DERAP (3) *GX (14) 
$ +DERAP (12) *l\'A2* (1 .-DERAP (12) *GX (14) ) /DM 
KTDOT2=KTDOTl/ (DM1*DM2) 

KT2=WA2* (1 .-DERAP (12) *GX (14) )+KTD0Tl*WA2*DERAP (12) *DERAP (15) 
KT2=KT1/DM 

XNUI^1= ( (DERAP (9) /GX (14) ) **2.0) * (KT2*KT2-KTDOT2*KTDOT2* (2.*WA*DZETA 
$+l^fA2*GX (11) /(GX (14) ) ) ) 

TSTO=WA2*GX (11) /GX (14) 

DENOM=TSTO- (1 .+2 .*DZETA*TSTO/li!A) * (2 , 0*DZETA*m+TSTO) 

SIGA= (ABS (XNUM/DEMOM) ) **0 ,5 
COM(40)=7.0*SiaA 
GX(13)=7.0*SIGA 
C 

C ASSUME THAT APPROACH RT-IS IS 7 FT/SEC 

C 

CON(41)=KTDOT2 
COM(42)=KT2 
TA=HX (10) 

TUR=GX(11)/U0CR 

A0=TUR*TA 

A2=TA* (2 . *DZETA*WCR4WCR2*TUR) +1 .+2 . *DZETA*WCR*TUR 
A1=TA* ( 1 . +2 . *DZETA*WCR*TUR+TUR) 

A3=WCR2* (TA+TUR) +2 , *DZETA*WCR 
A4=WCR 

W1=-A1*A4+A2*A3 
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XNUM= (1 . AA) * (DERCR (9) /UOCR) *KTD0T1*KTD0T1*W1 
XNUM=XNUM-A3*TUR* (DERCR (9 ) /UCX:R) ** 2 . 0*KT1*KT1 
DEN=A0*A3*A3+A1* (A1*A4-A2*A3) 

SIGDOT= (ABS (XNUM/DEN) ) **0 . 5 

COM(43)=3.0*SIGDOT 

TUR=GX(11)/GX(14) 

AO=TUR*TA 

A2=TA* (2.*DZETA*WA4WA2*TUR)+1 .+2.*DZETA*TUR*WA 
A1=TA* (1 .+2 . *DZETA*WA*TUR+TUR) 

A3^A2* (TA+TUR) +2 . *DZETA*WA 
A4=WA2 

W1=-A1*A4+A2*A3 

XNUM=XNUM-A3*TUR* (DERAP (9) AOCR) **2.0*KT2*KT2 
DEN=AO*A3*A3+Al* (A1*A4-A2*A3) 

SIGDOTA= (ABS (XNUMA>EN) ) **0 . 5 

CON(44)=7.0*SIGDOTA 

***************************** 

* OUTPUT SECTION * 

IF(OUTPUT.EQ.O) GO TO 999 
VJRITE (6,102) 

102 FORMAT (*1V//20X*AIRCRAFT OPTIMIZATION CONSTRAINTS*//5X 
$*DESIGN CONSTRAINTS*/14X*ID*5X*CONSTRAIOT*T49*VALUE* 
$9X*SL*10X*SU*9X*VIOLATION?* ) 

WRITE(6,112) (I,CON(I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=1,6) 

112 FORMAT (/10X,I5,* CRUISE THRUST*T45,4F12.4/10X,I5,* 2ND SEGMENT C 
$LI^©*T45,4F12.4/10X,I5,* MISSED APPROACH CLIMB*T45,4F12.4/10X,I5 
$* LANDING*T45,4F12.4/10X,I5* TAKE-OFF*T45,4F12.4/10X,I5 

$* LANDING GEAR LIMIT*T45,4F12.4) 

WRITE (6, 116) CON (45) ,SL(45) ,SU(45) ,XINEQ(45) 

116 FORMAT(13X*45 PASSENGER VOLUME*T45,4F12.4) 

WRITE (6, 117) CON (48) ,SL(48) ,SU(48) ,XINEQ(48) 

117 FORMAT (13X*48 CRUISE ALTITUDE*T45,4F12.4) 

WRITE(6,118) CON(49) ,SL(49) ,SU(49) ,XINEQ(49) 

118 FORMAT (13X*49 CRUISE ALTITUDE (L/D (MAX) ) *T45 ,4F12 .4) 

WRITE(6,119) CON(50) ,SL(50) ,SU(50) ,XINEQ(50) , 

$ CON(51),SL(51),SU(51),XINEQ(51) 

119 FORMAT (13X*50 CRUISE VJING CL*T45,4F12.4/ 

$ 13X*51 APPROACH IVING CL*T45,4F12.4) 

WRITE(6,120) CON(52) ,SL(52) ,SU(52) ,XINEQ(52) 

120 FORMAT(13X*52 AR(TAIL)/AR(WING)*T45,4F12.4) 

WRITE (6, 122) 

122 FORMAT (//5X*HANDLING QUALITY CONSTRAINTS*) 

WRITE(6,128) (I,CON(I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=7,10) 

128 FORMAT (/lOX, 15* STATIC STAB. (CR)*T45,4F12.4/10X,I5 
$* STATIC STAB. (AP)* 

$T45,4F12.4/10X,I5* MANEUVER MARGIN (CR)*T45,4F12.4/10X,I5 
$* MANEUVER MARGIN (AP)*T45,4F12.4) 

WRITE (6, 132) (I, CON (I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=11,15) 

132 FORMAT (lOX, 15* TAIL LIFT (AP)*T45,4F12.4/10X,I5* NOSE GEAR UNSTI 
$CK*T45,4F12.4/10X,I5* DYN. STAB. (CR)*T45,4F12.4/10X,I5 
$* DYN. STAB. (AP)*T45,4F12.4/10X,I5* PHUGOID FREQ (CR)*T45, 
$4F12.4) 
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WRITE(6,136) (I,CON(I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=16,20) 

136 FORMAT (10X,I5* PHUGOID FREQ (AP)*T45,4F12.4/10X,I5 

$* PHUGOID DAMPING (CR)*T45,4F12.4/10X,I5* PHUGOID DAMPING (AP)* 
$T45,4F12.4/10X,I5* SHORT PER. FREQ. (CR)*T45,4F12.4/10X,I5 
$* SHORT PER. FREQ. (AP)*T45,4F12.4) 

WRITE(6,141) (I,CON(I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=21,25) 

141 FORT-1AT(10X,I5* SHORT PER. DAIIP (CR)*T45,4F12.4/10X,I5 

$* SHORT PER. DAT'IP (AP)*T45,4F12.4/10X,I5* TIME-TO-DOUBLE (CR)* 
$T45,4F12.4/10X,I5* TIME-TO-DOUBLE (AP)*T45,4F12.4/10X,I5 
$* TIME-TO-HALF (CR)*T45,4F12.4) 

WRITE(6,146) (I,CON(I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=26,30) 

146 FORMAT (10X,I5* TIME-TO-HALF (AP)*T45,4F12.4/10X,I5 

$* FLIGHT PATH STAB. (AP)*T45,4F12.4/10X,I5* VERT. GAIN (CR)* 
$T45,4F12.4/10X,I5* VERT. GAIN (AP)*T45,4F12.4/10X,I5 
$* T(THETA{2)) (CR)*T45,4F12.4) 

WRITE(6,147) (I,CON(I),SL(I),SU(I),XINEQ(I),I=31,40) 

147 FORMAT (lOX, 15* T(THm(2)) (AP) *T45,4F12.4/10X,I5* VvV//NZA (CR)* 
$T45,4F12.4/10X,I5* UW/NZA (AP)*T45,4F12.4/10X,I5* T(l) (AP)* 
$T45,4F12.4/10X,I5* MODE RATIO (CR)*T45,4F12.4/10X,I5 

$* MODE RATIO (AP) *T45,4F12.4/10X,I5* ELE. VAR. (CR)*T45,4F12.4 
$/10X,I5* THETA-DOT GAIN (CR)*T45,4F12.4/10X, 15* THETA GAIN* 
$T45,4F12.4/10X,I5* ELE. VAR. (AP)*T45,4F12.4) 

WRITE(6,152) (I,CON(I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=41,44) 

152 FORMAT (lOX, 15* THETA-DOT GAIN (AP)*T45,4F12.4/10X,I5 

$* THETA GAIN (AP)*T45,4F12.4/10X,I5* ELE-DOT VAR. (CR)*T45, 
$4F12.4/10X,I5* ELE-DOT VAR. (AP)*T45,4F12.4) 

VJRITE(6,157) (I,CON(I) ,SL(I) ,SU(I) ,XINEQ(I) ,1=46,47) 

157 FORMAT (10X,I5* TRIM ELEVATOR (CR)*T45,4F12.4/ 

$10X,I5* TRIM ELEVATOR (AP)*T45,4F12.4) 

WRITE (6,167) 

167 F0RT>1AT(1H1/) 

999 RETURN 
END 

SUBROUTINE DIMDER(U,ICR,C) 

REAL C(15) ,A(15) ,ANS(4) 

C0MT40N /DEVARA>ESIGM(15) ,ITERM(10) ,CST(10) 

Ca-IMON /DRAG/CDS (6) ,CDSAP (6) 

COr^lMON /GEOM/W(20) ,HX(20) ,GX(35) ,PX(15) 

COMMON /STAB/DERIVCR(15) ,DSRIVAP(15) ,STOR(20) 

CC^IMON /WTSVE/liJTS(20) 

************************************************************************ 

* DIMDER CONVERTS FROM NON-DIMENSIONAL TO DIMENSIONAL STABILITY * 

* DERIVATIVES. * 

************************************************************************ 

IF(ICR.EQ.O) GO TO 20 
DO 15 1=1,15 
15 A(I)=DERIVCR(I) 

XMU=OTS(7) 

CALL AT62 (WTS (11) ,ANS) 

DRAG=CDS(5) 

XIY=l*/rS(13) 

CL=GX(9) 

WI^WTS(15) 
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GO TO 30 
20 DO 25 1=1,15 
25 A(I)=DERIVAP(I) 

XMU=WTS(9) 

CALL AT62(500.,ANS) 

DRAG=CDSAP(5) 

XIY=l'7TS(14) 

CL=GX (10) 

V7T=WrS (1) *GX (27) 

30 T1=ANS (1) *DESIGN (1) *0*32.174/1-^ 

C (1) =T1* (-A(l) -DRAG) /2.0 

C(2)=Tl*(CL-A(2))/2.0 

T2=ANS (1) *DESIGN (1) *U/(2.0*XIY) 

C(3)=A(3)*W(10)*T2 
C(4)=T1* (-CL-A(4)) 

C(5)=T1*(-DRAG-A(5)) 

C(6)=W(10)*A(6)*2,0*T2 
C (7 ) =-VJ ( 10 ) *A (7 ) *2 . 0*T2 
C(8)=0. 

C (9)=W(10) *W(10) *A(9) *T2/2.0 
C(10)=-A(10)/(4.0*Xf'5U) 

C(11)=0. 

C (12) =ANS (1) *DESIGN (1) *W(10) *W(10) *A(12) / (4 .0*XIY) 
C(13)=-Tl*U*A(13)/2, 

C(14)=-Tl*U*A(14)/2, 

C(15)=T2*U*W(10)*A(15) 

RETURN 

END 

SUBROUTINE LONGRT (DIM , U , ROOT , PARAM , A , NO , TCONST) 

REAL DIM(15) ,PARAM(4) ,A(5) ,OM(2) ,ZET(2) 

REAL TCONST (4) 

COr^PLEX ROOT(4) ,COM(4) 

************************************************************* 

* LONGRT FINDS ROOTS OF FOURTH ORDER DYNAT4ICS MODEL * 
************************************************************* 

A(l)=l. 

A(2) =-DIM (9) -U*DIM(12)-DIM(1) -DIM(5) 

A(3)=DIM(1) *DIM(9)-DIM(3) *U-DIM(2) *DIM(4)+DIM(5) * (DIM(9)+U*DIM(12) 
$+DIM(l)) 

A(4)=-DIM(5)*(DIM(1)*DIM(9)-U*DIM(3))+DIM(4)*(DIM(2)*DIM(9) 

$+32 , 174*DIM (12) ) -DIM (6) * (U*DIM (2) -32 . 174) 

A(5)=32,174* (DIM(4)*DIM(3)-DIM(6)*DIM(1)) 

IDEGRE=4 

CALL RPOLY(IDEGRE,A,ROOT,IERR) 

IF(IERR.LT.O) WRITE(6,22) IDEGRE 

22 FORMAT (*l*///*ROOT SOLVER BROKE DOWN— ONLY*I5* TERMS FOUND*//) 
OM(1)=0. 

OM(2)=0. 

ZET(l)=lo 

ZET(2)=1. 

DO 25 1=1,4 
TCONST (I) =0. 

PARAM (I) =0. 
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o o o 


25 CONTINUE 
ICOM=0 
ICON=0 


ASSUME THE ROOTS ARE GIVEN IN PAIRS 


DO 35 IROOT=l,4 
D=AIMAG (ROOT(IROOT) ) 

IF(D.EQ.O.O) GO TO 30 

ICOM=ICOMfl 

COM ( ICOM) =ROOT ( IROOT) 

GO TO 35 
30 ICON=ICON+l 

TCONST (ICON)=ROOT (IROOT) 

35 CONTINUE 
JCOM=0 

IF(ICOM.LT,l) GO TO 200 
DO 40 1=1, ICOM, 2 
D=AIMAG(COM(I)) 

E=COM(I) 

JCOM=JCOM+l 

OM ( JCOM) = (D*D+E*E) **0 . 5 
ZET (JCOM) =-E/OM (JCOM) 

40 CONTINUE 
200 NO=ICON 

****************************************** 

* PHUGOID(l,2); SHORT PERIOD (3,4) * 

****************************************** 


IF(OM(l) oGT.OM(2)) GO TO 250 
210 PARAT4(3)=OM(2) 

PARAT4(4)=ZET(2) 

PARAT4(l)=OM(l) 

PARAM(2)=ZET(1) 

GO TO 275 

250 IF(OM(l) .LT.0.8) GO TO 210 
PARAT4(3)=OM(l) 

PARAr-l(4)=ZET(l) 

PARi'J'l(l)=0M(2) 

PARAM(2)=ZET(2) 

275 IF(NO.LT.l) GO TO 999 
sri^L=ioo. 

SMALL2=100. 

DO 300 1=1, NO 

IF(TCONST(I) .LT.SMALL) SMALL=TCONST(I) 

300 CONTINUE 

DO 310 1=1, NO 

IF(TCONST(I) ,LT«SMALL2<,ANDoTC0NST(I) ,NE. SMALL) 
$ SMALL2=TCONST(I) 

310 CONTINUE 

R12=ABS (SMALL*SMALL2) 

ZT=- (SMALL+SMALL2) /(2o0*R12**0 .5) 
IF(R12.EQ.0.0) ZT=0. 

IF(PARAM(3) .EQ.O.) PARAM(3)=R12**0.5 
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IF(PARAM(4) .EQ.1.0) PARM(4)=ZT 
999 RETURN 
END 
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APPENDIX II - SAMPLE INPUT DECK 


MIL 8785B LEVEL III (SIMPACT) 


(Title) 


7 52 

-.2170 .3383 -.1576 


1 

1000. 

4000. 

2 

3. 15 

• 

3 

120. : 

260. 

4 

100. : 

1600. 

5 

2. 15 

• 

6 

10000. 

120000, 

7 

-0.5 : 

L. 

1 

1, 2. 


2 

1. 5. 


3 

1 0 5 O 


4 

1000. 

8000. 

5 

2000. 

10000. 

6 

0. 1.1 

0 

7 

-999. 

999. 

8 

-1.0 

-0,10 

9 

-999. 

999. 

10 

-999. 

999. 

11 

-.8 

.8 

12 

1. 3 

.00 

13 

-999, 

999. 

14 

-999. 

999. 

15 

-999. 

999. 

16 

-999. 

999, 

17 

-999. 

999. 

18 

-999. 

0 999.0 

19 

-999. 

999. 

20 

-999, 

999. 

21 

-999. 

0 999.0 

22 

-999. 

0 999.0 

23 

-999. 

999. 

24 

-999. 

999. 

25 

-999. 

999. 

26 

-999. 

999. 

27 

-3. 

.24 

28 

-999. 

999. 

29 

-999. 

999. 

30 

-999. 

999. 

31 

-999. 

999. 

32 

-999. 

0 999.0 

33 

-999. 

0 999.0 

34 

-999. 

999. 

35 

-999. 

999. 

36 

-999. 

999. 

37 

-999. 

999. 

38 

-999. 

999. 


(Number of design variables. Number of constraints) 
-.1123 -.3620 .0358 .1895 (Initial XBAR) 

(Design variable number, lower bound, upper bound) 


(Constraint number, lov/er bound, upper bound) 
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39 

-999. 

999. 

40 

-999. 

999. 

41 

-999. 

999. 

42 

-999. 

999. 

43 

-999. 

999. 

44 

-999. 

999. 

45 

1. 2.1 

5 

46 

-1. 1 

• 

47 

-1. 1 

.0 

48 

30000. 

46000. 

49 

30000. 

52000. 

50 

.1 .75 

51 

1. 2. 

7 

52 

0. 1. 


2 



1. 

l.OE-6 

I.OEIO 

190000. 


0. 

0. 0. 

• 

O 

• 

o 


(Number of calls to NELMIN) 
(SCF, REQMIN, CAYY, STEP, ILINE) 

( WTS(l) ) 
( PX(1— 8) ) 


1 0 

7 2 1 

1979 

0 

0 0 0 

( ITERM(1— 

10) 

) 

.75 

.65 




( V/TS(16) , CX5(5) 

) 

21.2 

2.0 .38 

.14 

5 

• 

( W(l- 

- 5) 

) 

.98 

.850 0.3 

-.15 


-1. 

( W(6— 8, 14, 

16) 

) 

-.15 

-.25 -. 

12 .025 


( W(17— 

20) 

) 

.4 

.1 .46 . 

2 



( HX(1— 3, 

10) 

) 

.08 

.36 1.8 

35. , 

.3 


( HX(16— 

20) 

) 

.8 

3000. 16. 

667 22. 

58 8.33 

( GX(3- 

- 7) 

) 

2.2 

1443.38 

3.15 


-5.25 -12.725 

( GX(8, 11, 12, 17, 

18) 

) 

200. 

8870. 41100. 

7500. 

( GX(19~ 

22) 

) 

4.0 

.55 .09 

3200. 


.70 

( GX(23~ 

27) 

) 

30. 

0. -99. 




( GX(32, 34, 

16) 

) 


tu 



APPENDIX III - SAMPLE OUTPU^ CORRESPONDING 


HUN HO* ML 87690 LEVEL III tSXHPACT) 


NO* OE VARIABLES* 
HO. OE CuNSIRAlNTS- 
VARIABLES* 


J»XL»XO»AVE»A«P* 

J,XL»XU»AVE*AflP* 

J»XL»XO»AVE»AMP* 

J»Xl»XU»AVE»AHP* 

J.XL»XU»AVE»AnP* 

J*Xl»XU»AVE»ANP* 

J»Xl»XU.AV£>AMP* 

I>SL(I>;SU(I)* 

I,SHI).SU«I»* 

i>SLin*su(n* 

tfSi(n>su(i)* 

i.sLn)»sum* 

i»snn*su(i)* 

I.SUD.SIMD* 
I,SL(t)«Sli(II* 
ifSiiittSod >* 
i»sLtn.su(i>* 
i,SLin»su(n* 
I»SltI».SU»ll* 

i>sL(n»su(ii« 

i.sLiiusum* 

i.SL(n»su(i>* 

i.sicn.sud)* 

i»sidwsud>* 

I»Mdl»SUd)* 

I,Sld)»SUd)* 

I»SLd»»SUd>. 

I,SLd»»SUd)* 

I,Sld)»S(MI)* 

l»Sld>*SUd)* 

(>SLd)>SUd)* 

i,sL(n.sud)* 

I>Sld>>SUdt* 

I»Sld>.SUd)* 

l»SLd»»SUd)* 

I»SLCU»SUdt* 

i>si.d)«sud)* 

l»5l(H»S Jd)* 
I»SLdt.lUd»* 
I.SLdl»SUd)* 

i>sid)«sud)> 


I»SLd)>SUdI. 
I»SLd)»$Ud)* 
I»SLd>i>SUd)* 
I»SLd).SUd)* 
I#$Ld)>SUd )* 
i»SLd}.su(n* 
I>SLdl«SUdl* 
l«SLd)>SLMl)* 
WSLdItSUd)* 
ifSiditSud)* 

i*$im»sud)* 

I*SLd)»SUdl* 
IrSldWSUm* 
l.SUd»SUd)* 
I»SLd)iSUdi* 
I>Sld)>SUd>* 
I»SUl)»SUd>* 
l»Sld»»SUd>* 
HUNBER OF 


92 

-*2170 

.0398 

1000.0000 

3.0000 
120.0000 
100.0000 

2.0000 


6 

lOCOO 

.0000 1 

7 

- 

.5000 

1 

1*C0 

2.00 

2 

1«00 

5.00 

3 

1.00 

5.00 

A 

1000*00 

OOOOeOO 

5 

2000.00 

10009.00 

6 

o.co 

1.00 

7 

-999.00 

999.00 

0 

-1.00 

-.10 

9 

-999.00 

999.00 

10 

-999.00 

999.00 

11 

-.80 

• 00 

12 

1.00 

3.00 

13 

-999.00 

999o00 

1^ 

-999.00 

999.00 

10 

-999.00 

990.00 

It 

-999.00 

999.00 

17 

-999.00 

999.00 

10 

-999.00 

999.00 

19 

-999.00 

999.00 

20 

-999.00 

999.00 

21 

-999.C0 

999.00 

22 

-999.00 

999.00 

23 

-999.00 

999.00 

2^ 

-999.00 

999o00 

25 

-999.00 

999.00 

26 

-999.00 

999.00 

27 

-3.00 

.24 

20 

-999.00 

999,00 

29 

-999.00 

999.00 

30 

-999. CO 

999.00 

31 

-999. CO 

999.00 

32 

-999.00 

999.00 

33 

-999.00 

999.00 

34 

-999.C0 

999.00 

35 

-999.00 

999.00 

36 

-999.00 

999.00 

37 

-999.00 

999.00 

30 

-999.00 

999.00 

39 

-999.00 

999.00 

40 

-999.00 

999.00 

41 

-999.00 

099.00 

42 

-999.00 

999.00 

43 

-999.00 

997.00 

44 

-999.00 

999.00 

45 

1.00 

2.50 

46 

-1.00 

1.00 

47 

-l.CO 

1.00 

40 

30000.00 

46000.00 

49 

30000.00 

52000.00 

50 

.10 

• 75 

51 

l.CO 

2.70 

52 

OoOO 

1.00 


.3363 
.1009 

AOOO.OOOO 
19.0000 
260. COCO 
1600.0000 
19.0000 
00.0000 
1.0000 


-.1976 -.1123 


2900.0000 

9.0000 

190.0000 

690.0000 
e.9000 

69000.0000 

*2900 


REOUESTEO NELNINS* 


TO SAMPLE INPUT 


-.3620 

1900.0000 

6.0000 

70.0000 

790.0000 

6.9000 

99000.0000 

*7900 


T5 



/\/N/\FUNCTXON lNPUT/\/N/\ 


WTSm 

1T£PH 

PX- 


VTSC16) 

W«#*17-20» 

HX# •• 10* 

HX...16-20- 

CX...3-7- 

GX,*.0#U-l2#17,ie- 

GX«.«19-2?- 

CX..e23-27* 

GX*»,32,3A- 

INITAIL XBAR 

INITAII STEPS 

REQHIN»SCF,CAYY 

/\/\/\NELMlN C0HPlE7E/\/\/\ 

XMXN- 

YNEWLO* 

ICnUNT- 

TOTo FUNCTION CAU5» 


TOT, WEIGHT ITERATIONS 
AVE, WT, ITERATIONS PER 

INITAIL X6AR 

IHITAIL STEPS 

REQNIN#SCF/CAYY 

/\/\/\NELMIN CnrtPLETE/N/\/\ 


XttIN* 

YNEWLO* 

ICOUNT- 

TOT« FUNCTION CALLS- 


TOT* WEIGHT ITERATIONS 
AVEo WTo ITERATIONS PEP 


19OOOO*OC00 
10 7 

2 1 

1979 0 

0 0 0 


0*0000 

0*0000 

0.0000 

0*0000 

0*0000 

0*0CD0 

0*0000 

0*0000 

X 

I 

I 

1 

X 

I 

I 

*7500 

21*2000 

2*0000 

• 3800 

*1400 

5*0000 

*9600 

*8500 

*3000 

-*1500 

-1*0000 

-*1500 

•*2500 

•*1200 

• 0250 


*4000 

*1000 

• 4600 

• 2000 


•oeoo 

• 3600 

1*8000 

35*0000 

• 3000 

*8000 

3000*0009 

16.6670 

22*9800 

8*3300 

2*2COO 

1443*3800 

3*1500 

-5*2500 

-12*7250 

200*0000 

8870*0000 

41100*0000 

7500*0000 


4*CCOO 

*5500 

• 0900 

3200*0000 

• 7000 

30*OCCO 

-•2170 *3383 

0*0000 

-*1576 

•*1123 

•*3620 


*0358 *1895 

*2000 *2000 

*2000 

*2CC0 

• 2000 


*2000 *2000 
•lOOOE-06 *1000E401 *1000E+12 

-•15129477 *18686605 

•15769764 

*05231575 

*48770473 

•20364349 *16912692 

*200533E^04 
817 





7 

2973.00 





3*604 





-*1513 

*1869 

-.1577 

*0523 

-•4877 

*2036 

• 1691 




*1000 

*1000 

• 1000 

*1000 

• 1000 

*1000 

*1000 




lOOOE-07 

•lOOOE^Ol 

•1000E413 




-•151^7108 *18662137 •*15771096 *03225537 •*^8798386 

•20313425 *16913646 

*200307E404 
464 
1317 


4421*00 

3*357 


76 



AIRCRAFT SIZING FRQGRAN 


DFSICN VARIABLES 

WING APIA IFTXJf?)- 
WING aspect RATIO- 
FUSELAGE LENGTH (FTI- 
HOR* TAIL AREA (FTXX2I- 
HOR. TAIL ASPECT RATIO- 
TOTAL THRUST UBS)- 
AFT POST CC- 
CRUISE HACH NO.- 
SWEEP- 
WING T/C- 
WING TAPER RATIO- 
FUSE. OIA- 


2146. 4e06 
10.7S38 
172.8356 
911.4627 
3.9914 
62253.2621 
• 4469 
I 
X 
1 
X 
X 


input constants 

WING.. .SWEEPi INC IDFNCE# TAPER RATIO! 
THrCKNESS>TWlSTf El# E2# DESIGN CLi 
CH(CR»APP) 

DELTA Cn(10#25 DEGREES FLAP)! 

ANGLE or ZERO LIFT(0#10#45 DEGREES FLAP)! 
DELTA CO ao-45 DEGREEf FLAP)! 

TURBULENCE LENGTH/ROOT 3#FUSE. OIA.t 
CL-HAXCTO)#CL-HAX(L)l 
nX:.SION...HACM Nn.#RANGE»NO. PASS! 

CARGO UEIGHTl 

DELTA CG» WTL(NAX)/WTOl 

ENGINE...L#W#WT#TREFi 

TAIL...TAPER RATI0#THICKNES$#ELE EFFI 

ELEVATOR TINE CONSTANT! 

ECONOnlCS...LOAO FACT#S/SEAT PX#6LK HR/YR! 
SOME GEOMETRY CALCULATIONS 

VlNG...$PAN»CNACt 

T AIL.. .SPAN#CH AC# VBAR# SWEEP 

VERT. TAIL.. .VBAR# TAPER# AR# SWEEP# SR/SVfSVI 

444WEIGHT ESTIHATI0N44* 


WING 37077.8 
HOP. TAIL 702B.3 
VERT. TAIL 1066.5 
FU9HAGE 47372.1 
LAKDIHG GEAR' 14912.5 
CONlPni SYSTEH 2541.3 
INSTRUMENTS 116.2 
NISC INTERIOR 250.6 
ELECTRICAL 2192.3 
STARTERS 506.9 
FURNISHINGS 7902.0 
OXYGEN 300.7 
WINDOWS 501.6 
BAGGAGE HNOLING 144.7 
AIR CONDITIONING 3647.C 
THRUST REVERSER 1500.0 
ENGINE 18639.1 
ENGINE CONTROLS 219.5 
FUEL SYS 2571.1 
ACTIVE CONTROL SYSTEMS 250.0 

empty weight 156200.2 


PASSENGERS 34000.0 

CREW 1700.0 

BAGGAGE 7000.0 

UTFOOO 214.5 

CARGO 7500.0 

FIXED WEIGHT 50414.5 

/ 


FUEL 74615.1 


TAXE*OFF WEIGHT 281429.7 


21.2000 

2.0000 

.3600 

• 1400 

5.0000 

• 9600 

-.1500 

-.1207 


-.1500 

-.2500 


-1.0000 

-5.2500 

-12.7250 

.0250 

1443.3600 

16.6670 


2.2000 

3.1500 


.6000 

3000.0000 

200.0000 

7500.0000 

4.0000 

•7000 


22.5600 

6.3300 

6670.0000 

• 4C00 

• 1000 

• 4600 

• 2000 

.5500 

.0900 

3200.0000 

151.7661 

15.0926 


60.3166 

16.0370 

2.1663 

• 0600 

• 3600 

1.6000 


• 8500 


41100.0000 


30.0000 

35.0000 


.3000 


.3000 


NO. OF ITERATIONS REQUIREO 


1 


WTNTC1#2#3)#WTFUSE(1#2)#WTINIT- 


6635.45 7629.61 6419.59 57651.36 36692.76 261429.74 




3J2.T617 
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/\/\/\CRUlSC ANAlYSlS/\/\/\ 


TOTAL hISSION RANGE- 3000*00 
CLIHB DISTANCE- 189,00 
DESCENT distance- 113*00 


LEG 

HACH NO* 

CL 

V 

tine 

L/0 

tsfc 

T/T(IN) 

ALTC8EG) 

1 

*79 

*593 

763*4 

*597 

16*76 

*471 

*233 

36CC0*0 

2 

• 80 

*606 

774*8 

*568 

18*58 

• 478 

*217 

37398*1 

3 

* 60 

*606 

774,5 

*568 

16*55 

• 479 

*214 

37714*4 

4 

• 60 

,606 

774,5 

* 56 6 

16*53 

*480 

*211 

36031*6 

5 

*60 

*606 

774*5 

• 566 

18*50 

• 480 

*206 

38349*6 

6 

• 60 

• 606 

774*5 

*568 

16*46 

• 461 

*205 

36668*4 

7 


*674 

774*5 

• 588 

18*31 

*486 

*182 

41197*6 

0 

* 60 

*662 

774*5 

*566 

16*26 

*467 

• 176 

41755*7 

9 

• 80 

• 687 

774*5 

*566 

16*21 

• 486 

• 172 

42272*6 

10 

*60 

• 669 

774*5 

*568 

16*16 

*468 

*169 

42657*4 

11 

*63 

*569 

644*4 

2*620 

19*00 

•425 

• 300 

30000*0 


9*«fuel weight ANALYS1S«^9 


TAKE-OFF 

261429*74 

. start-cruise 

263432,31 

KIO-CRUISE 

244203,64 

END-CRUISE 

225649*37 

AFTER RESERVE 

213005*63 

AFTER OESCENT/TAXl 

206614*66 

NET FUEL weight UBS) 

74815*05 


CRUISE ALTITUDES 


LEG 1 

360C0*OO 

LEG 2 

43026*59 

RESERVE LEG 

30000*00 


FIIGHT length (HR) 

6,67 

AVERAGE SPEED IKTS) 

449*82 

SLOCK TINE (HR) 

7*00 

BLOCK SPECO (KTS) 

426*60 

fcLOCK FUEL (LBS) 

44358*42 

CLOCK FUEL (GALS) 

6931*00 

NAUT* hl/GAL 

*43 

NAUT* SEAT Pl*/GAL* 

86*57 

INSTALLED THRUST (LBS) 

62253*26 

NO* OF ENGINES 

2 

ENGINE THRUST (LBS) 

41126*63 

REFERENCE ENGINE (LBS) 

41100*00 

SCALE FACTOR 

1*001 


ALT(END) 

VT(6EG) 

WT(END) 

OIST 

ICON 

GANA 

CLH 

37398*1 

263432*3 

259519*3 

269*8 

1*346 

• 0007 

•751 

37714*4 

259519.3 

255617*9 

269*6 

1.279 

0*0000 

• 733 

38031*6 

255617*9 

251765*1 

269*6 

1*279 

0*0000 

*734 

36349*6 

251765*1 

247960*6 

269*6 

1*776 

0*0000 

• 735 

36668*4 

247960*6 

244203*6 

269*8 

1*276 

0*0000 

• 735 

41197*6 

244203*8 

2A0A9A.3 

269*6 

1.244 

*0013 

• 736 

41785*7 

240494*3 

236771*9 

269*6 

1*134 

• 00^2 

• 743 

42272*6 

236771*9 

233067*9 

269*8 

1*117 

•0001 

.749 

42657*4 

233067*9 

229545*6 

269*6 

1*106 

*0000 

• 751 

43026*6 

229445*6 

225849*5 

269*6 

1*101 

•ocoo 

• 752 

30000*0 

225649*4 

213005*6 

1000*0 

2*076 

0*0000 

•720 


78 



•♦♦DRAG ANALYSXS««^ 


WING «C07e 
HORIZONTAL TAIL .OOZ9 
VERTICAL TAIL •0011 
FUSELAGE tOOZS 
BASE •OOOZ 
ENGINE NACELLE *0013 
CRUO/FLAPS *0007 

AIRCRAFT DRAG .OIBA 


INTERFERENCE FACTOR IS 9 PERCENT 


#4#REYN0LD*S NUHBERS*«« 


WING 

HORIZONTAL TAIL 
VERTICAL TAIL 
FUSELAGE 
ENGINE 


ZA67AZ27.X 

Z6Z16160*! 

24976l59o6 

2625609B9*? 

369Z6961»6 


STABILITY AND CONTROL DERIVATIVES . 


CL*H#C6« POSITION- «606 •600 *3X4 



Cl 

CD 

CN 

ALPHA 

6^69956 

•2B332 

-X.5X107 

VELOCITY 

• 2X59 

• 00X0 

-•0003 

0 

X6^85623 

0.00000 

-70.17792 

ALPHA-DOT 

-9«190X9 

0.00000 

-26.74639 

ELEVATOR 

•66393 

•00004 

-3.9Z446 


NEUTRAL POINT 

• 940 


STATIC STABILITY 

-•226 


CN(ALPHA) 

-loMX 


STOR» 




5.7X8 

3.890 


• 339 

• 940 


• 604 

4.669 


3«09X 

• 290 


.947 

• 737 


-.028 

-.017 


.172 

• 209 


I 

X 


1 

1 


I 

1 


-.1729 

-.2939 

• 1962 

• 0998 

3232.9997 

• 6994 

-.2769 

• 0040 

9.7163 

3.8903 

.2724 

77.7760 

0.0000 

3.9033 

.6000 

.1720 

9.7X79 

• 3367 

• 9400 

• 9902 

1.2657 

X.6X39 

3.739X 

• 9997 

• 3X44 


• 0758 
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ACCURATE L/D ANALYSIS 


CL(RtOUESTED)- 

• 606 

CLIUINO- 

• 655 

CL(TAIL>- 

••118 

ELFVATOR(OEGREES)- 

0^00 

STAOILIZERIDEGREES) 

-3.09 

INDUCED DRAG CONPOHENTS 

WING- 

• 01A6 

XNTCRFERENCE- 

-•0019 

TAIL* 

• C0C6 

HRIH)* 

• 0009 

ITAID- 

• 0036 

SIGMA 

.3715 

DRAG COEFFICIENTS 

INDUCED- 

• 0133 

ZERO LIFT- 

• 016A 

ELEVATOR- 

0.0000 

WINGinACHI* H(CRIT)- 

• 0011 

FUSE(NACH)* H(CRIT)- 

-•0000 

TOTAL- 

• 0326 

L/0- 

18.A77 


HAX. L/0« 

19.466 

CL-L/0 HAX- 

.736 

lED COO- 

•0X93 


AIRCRAFT COST tSTIHATES 



DEVELOPNFNT 

ENGINEERING 

137888092*30 

PFVELnPNENT SUPPORT 

93303827*90 

PLIGHT TEST 

17877799*66 

TOOLING 

166227160*03 

HANUFAC* LAPOR 

90339216*96 

QUALITY CONTROL 

11766098*21 

HATERIALS 

10693968.66 

ENGINE 

11723236*28 

AVIONICS 

738662*00 

ACTIVE CONTROLS SYSTEM 

907692*63 

TOTAL 

678803912*20 


TOTAL COST RER AIRCRAFT- » X6393A6I.72 


• eo 

• 90 


PRODUCTION 

196227689.86 

0.00 

0*00 

25288A765.99 

1068921091*06 

1363077A1,86 

6712A9929.87 

976936190*39 

92307790*00 

63661978*13 

3237892333*06 
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HAXNTEHANCE OPERATING COSTS 


NO* 

SYSTEM 

1 

XNSI' 

2 

AIR COND 

3 

AUTO PILOT 

4 

COKr^UN 

9 

ELEC 

6 

FURN 

7 

FIRE PROT 

e 

FLT CDNTL 

9 

FUEL 

10 

HYD POWER 

11 

ICE 

12 

INSTR 

13 

LAND GEAR 

14 

LIGHTING 

19 

NAVIG 

16 

OXYGEN 

17 

PNUEMAT 

16 

WAT/WASTE 

19 

AIR APU 

20 

STRUCTURE 

21 

DOORS 

22 

FUSELAGE 

23 

NACELLES 

24 

WINGS 

29 

STAB 

26 

WINDOWS 

LABOR 

COST 

nATEAIU C05T 

ENGINE 

LABOR COST 

ENGINE 

NATERIAl COST 


LABOR 

MATERIAL 

39*92 

9.76 

9*10 

4*92 

13*43 

3*19 

9*92 

2*36 

4*31 

9*79 

26*10 

13*93 

9*37 

1.46 

6*69 

9*49 

2*01 

1.02 

3*33 

3*99 

*67 

• 96 

1*22 

• 49 

13*64 

28*07 

3*61 

2*83 

11*99 

8*24 

1*09 

1*90 

1*28 

6*99 

1*0G 

1*94 

• 32 

*46 

3*70 

0*00 

2*39 

1.96 

4*76 

.90 

*67 

• 20 

2*99 

X.13 

• 83 

• 37 

.09 

7.24 


69*64 


43*64 


29*38 


31*26 


MAINTNENACE DOC JN 1976 DOLLARS PER HOUR 


DIRECT OPERATING COSTS— DOLLARS/FLT# HOUR 


OEPREC 

393*92 

19.63 

SUPPORT 

93*66 

2*68 

SPARES 

26*63 

1*34 

DELAY 

10*29 

• 91 

insurance 

62*61 

3.12 

FUEL 

743*00 

37.06 

NAlNTtNANCE 

166*13 

8*29 

LANDING FEE 

26*6C 

1.34 

CREW 

266*46 

13.29 

attendants 

164*11 

8*16 

FUEL SERVICE 

77ol0 

3.89 

CONTROL 

14*46 

• 72 

DIRECT OPERATING COSTS 

% 2009*07 

100.00 


INDIRECT OPERATING COSTS— OOLLARS/FLT. HOUR 


MAIN BURDN 

174.44 

19*21 

FOOD COST 

160*39 

Ie*93 

MOVIE 

34*32 

3*78 

PASS INS 

30*09 

3*31 

MISC PASS 

6* 30 

*69 

AOVERIISE 

97*64 

10.79 

COMMISSION 

139*79 

14*99 

RESER 

04*79 

9*33 

PASS liOLG 

99*28 

6*09 

BAG HDLG 

29*23 

2*78 

CARGO HDLG 

06.07 

9.40 

SERVICING 

10*07 

1*11 

L INDIRECT OPERATING COSTS 

9Oe;20 

100*00 


166*13 
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P£RFCRH*NCE FUNCTION SUNHARY 


REVENUE PER BLOCK HOUR 

TOTAL COST PER BLOCK HOUR 2913.3^ 

RETURN ON IKVESTHENT #1229 


1 OOC/HR 

2 DOC/FLT 

3 ROI 
h FARE 

5 SEAT-MI/GA 

6 L/DCNAX) 

7 hTOGW 
6 FARE 
9 PRICE 


2CC9.065 

I4027.966 

• 123 

• 096 
86.56S 
I9.A66 

281429*737 

31749.C00 

4336«000 


DEBUG OF NOSE GEAR UNSTICK 


ZACTf XHUFf 2ACLG- .2630 •0230 «530I 

ANS(l#2f 3#4)« *002377 2116.23B4I7 286.130000 1116.449901 

VSTALL#CLW/0*XLWf X«W- 22B.9332 .3642 48*2736 43279.0297 -449306.3408 

CI»UNTHRUfXCGiXLTR0*Cl2>C0IT,XLT- 71#5390 82253,2621 .1819-26773.6962 1.3723 -12.3X03 77*7760 

XLTAV*CONIl2),STOR(7)fSTOR(8)#W«2)* -95242.4919 2.9739 3.0719 .2913 2.0000 

GX(17)#HX(6)^V(10)#XLTAV2*TLGE#CLTAV- -3.2300 2.1883 13.0926 -79622.1333 13620.3364 


-1.8096 
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••♦DRAG ANAIYSIS44* 


WING ,0057 

H0RI70NTAL TAIL *0022 

VERTICAL TAIL ,0000 

rUSElAGE o0022 

BASE ,0002 

engine nacelle ,00X0 

CRUD/FLAPS ,0230 

AIRCRAFT DRAG ,0376 


INTERFERENCE FACTOR IS 9 PERCENT 


• ♦•REYNOLD'S NUnOERS^o.» 


WING 

HORIZONTAL TAIL 
VERTICAL TAIL 
FUSELAGE 
ENGINE 


6^>X039A0o0 

T3'jO9606^io3 

93036692o3 


♦ OOAPPPOACH^'^^ 


STABILITY AND CONTROL DERIVATIVES 


♦ •♦APPP0ACH6«.> 


CL#H#CG* 

POSITION* 

lo7X3 

• 392 


Cl 

CO 

CN 

ALPHA 

A,66A61 

•001A7 

-•A0650 

VELOCITY 

• 0529 

• 0007 

• 0005 

0 

13#219XS 

OeOOOOO 

-61,309A3 

ALPHA-DOT 

-3«A6736 

0,00000 

-17*9712S 

ELEVATOR 

•53636 

eOOOOA 

-2o76A00 


NEUTRAL POINT 

• 5A7 


STATIC STABILITY 

-•100 


CHIALPHA) 

-•^07 


STOB* 

5.7X8 

3.090 


• 339 

• 560 


• 60A 

A.865 


3^051 

• 290 


• 5A7 

• 737 


-•028 

-•017 


• 172 

• 205 


I 

I 


• 4A7 


I X 

X X 


-•2A48 

-•3976 

• 2670 

• 0705 

22C3«?770 

• 7036 

-•3996 

• 0037 

6.0372 

3»0509 

• 2726 

77o7760 

0«0000 

3.5033 

• 9016 

• 2050 

6o0363 

• 2902 

• 5669 

• 9902 

1.2657 

1,2293 

15.3662 

.9997 

•6669 

• 0753 
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ACCUPATE L/0 ANALYSIS 


APPROACH WITH A5 DEGREES FLAP 


ClCREOUtSTEO)- 

1.713 

CL(WtNG)> 

1*761 

CKTAIU* 

-•llA 

£lEVAT0RI0E6R6ESt* 

0*00 

STABlinERtOEGREESt 

-1*35 


INDUCED DRAG COHPONENTS 
WING- 

INTERFERENCE- 

TAIL- 

(TRIH)- 

(TAIL)- 

SIGMi 


• 1303 
-•C0A8 

oOOOS 

e0020 

oOOSO 

• 3715 


DRAG COEFFICIENTS 
INDUCED- 
ZERO LIFT- 
ELEVATOR- 

WING(HACH)> M(CRIT)- 
FUSE(NACH)# M(CRIT)- 
TOTAL- 


• 1260 
• 0626 

• GOOD 
OaOOOO 
OttOOOO 

• 1091 


OoOO 

OoCO 


I/O- 


9^057 


HAX^ L/0- 9^57A 

CL-L/D hAX- 1*350 

MODIFIED COO- •0655 


/\/\/\ClMENSIONAL STAPILITY DERI VATIVES/ V/\/\ 


QFRCK- 

-•A615 

• 0221 

-•OOOA 

-•1126 

-•00A6 

-.0000 

-.0197 

0.0000 

-.3AAA 

• 0030 

0.0000 

-•0002 

-36.3072 

-•0AA6 

-1.5935 


-.5723 

.1201 

-•OOOA 

-.3999 

-.0A30 

• 0000 

-.0208 

0.0000 

-.3632 

• OlOB 

0.0000 

-•OOOA 

-16.7513 

-.0263 

-.596^ 



lONCITUOINAL OYNAnICS 


/V/\CRUlSE/N/\ 


COEFFICIENTS- 1*000000 

*926777 

•443086 

•003445 

*001313 

R00TSCREAL#1HAGINARY) 
>•0008 •OS^A 
■•0008 -•054A 
•A634 .ATTO 
■•4634 -*4779 





PHUGOIO FREOUENCY 

• 0044 




FHUGOID DAMPING 

*0145 




SHORT PER. FPEO* 

*6656 




SHORT PER* DAMPING 

*6961 




NO* OF NON-OSCUIATORY ROOTS- 

0 




TIME CONSTANTS- 

0*0000 

0*0000 

0*0000 

0*0000 

/\/\APPROACH/\/\ 





COEFFICIENTS- 1*000000 

*977315 

*264314 

*035203 

*004930 

ROOTSIREAL* IMAGINARY) 
•*0321 *1277 
■•0321 -,1277 
•*4867 *2134 
•*4007 -*2134 





PHUGOID FREOUENCY 

• 1317 




PHUGOlO DAMPING 

• 2439 




SHORT PER* FREQ* 

*5332 




SHORT PER* DAMPING 

*9164 




HO* OF HON-OSCILLATDRY ROOTS- 

0 




TIME CONSTANTS- 

0*0000 

0*0000 

0*0000 

0*0000 

CL-MAX TO(U)#Cl-MAX TO (AC)#CL2 


2*200 

1*976 

1*372 

CL-MAX (V)# CL-MAX (AC)# CLA- 


3*150 

2*095 

1*713 
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AIRCRAFT OPTIMIZATION CONSTRAINTS 


DESICN constraints 

10 CONSTRAINT 


VALUE SL SU VIOLATION? 


1 

2 

3 

A 

5 

6 
A9 
A8 
A9 

50 

51 

52 


CRIIISE THAUST 

1.1009 

1.0000 

2.0000 

O.OCOO 

2ND SEGKENT CLIMB 

1.0000 

1.0000 

5.0000 

O.COOO 

HISSED APPROACH CLIMB 

1.3067 

l.OOCO 

5.0000 

O.OCOO 

LANDING 

7999.9572 

1000.0000 

eooo.cooo 

0.0000 

TAKE-OFF 

8106.5719 

2000.0000 

10000.0000 

O.OCOO 

LANDING GEAR LIMIT 

.0231 

0.0000 

1.0000 

0.0000 

PASSENGER VOLUME 

1.0000 

1.0000 

2.5000 

0*0000 

CPUISE ALTITUDE 

38666.6002 

30000.0000 

66000.0000 

0.0000 

CRUISE ALTITUDECL/D(HAXn 

62707.2395 

30000.0000 

52000.0000 

0.0000 

CRUISE WING CL 

.6552 

•icoo 

.7500 

0.0000 

APPROACH UING CL 

1.6799 

1.0000 

2.7000 

0.0000 

AR(TAIL)/AP(W1N6) 

.3719 

0.0000 

1.0000 

0.0000 


HANDLING QUALITY CONSTRAINTS 


7 

6 

9 

10 

n 

12 

13 

lA 

15 

16 
17 
16 

19 

20 
21 
22 
.’3 
2A 

25 

26 
27 
26 

29 

30 

31 

32 

33 
3A 

35 

36 

37 
30 
39 
AO 
A1 
A2 
A3 
AA 
A6 
A7 


STATIC STAB* CCR) 

STATIC STAB, (AP) 
MANEUVER MARGIN (CP) 
MANEUVER MARGIN (AP) 
TAIL LIFT (AP) 

NUSE GEAR UNSTICK 
DYN# STAB. (CR) 

OYN. STAB. (AP) 

PHUGOIO FRfO (CR) 
PHUGOTO FREO (AP) 
PHUGOTD DAMPING (CR) 
PHUGOIO DAMPING (AP) 
SHORT PER. FREO. (CP) 
SHORT pro. FPEO. (AP) 
SHORT PER. DAMP (CR) 
SHORT PER. DAMP (AP) 
TIME-TO-DOUBLE (CP) 
TIME-TO-DOUBLE (AP) 
TIME-TO-HALF (CR) 
TIMf-TO-HALF (AP) 
FLIGHT PATH STAB. CAP) 
VERT. GAIN (CR) 

VERT# GAIN (AP) 
T(THETA(2)) (CR ) 
T(THETA(2)) (AP) 

UU/NZA (CR) 

WU/NZA (AP) 

T(l) (AP) 
mode RATIO (CR) 
mode PATIO (AP) 

ELE« VAP. (CR) 
7HETA-OOT GAIN (CR) 
THETA GAIN 
ELE« VAP. (AP) 
THFTA-DOT GAIN (AP) 
THETA GAIN (AP) 
ELE-OOT VAR. (CR) 
ELE-OOT VAP. (AP) 

TRIM ELEVATOR (CR) 
TRIM ELEVATOR (AP) 


-.0930 -999.0000 999.0000 0.0000 

-.1000 -1.0000 -.1000 0.0000 ‘ 

-.1383 -999.0000 999.0000 0.0000 

-.2900 -999.0C0O 999.0000 0.0000 

-.3897 -.8000 .8000 0.0000 

2.9739 1.0000 3.0000 0.0000 

.0000 -999.0000 999.0000 O.OCOO 

.0002 -999.0000 999.0C0O O.OCOO 

.0566 -999.0000 999.0000 0.0000 

.1317 -999.0000 999.0C00 0.0000 

,0165 -999.0000 999.0COO 0.0000 

.2639 -999.0000 999.C000 0.0000 

,6656 -999.0000 999.0000 O.OCOO 

•5332 -999.0000 999.0000 0.0000 

.6961 -999.0000 999,0000 O.OCOO 

.9166 -999.0000 999,0000 0.0000 

99,0000 -999.0000 999,0000 0.0000 

99,0000 -999.0CP0 999,0000 0.0000 

879.9306 -999.0000 999.0000 0.0000 

21.5825 -999.0000 999.0000 0.0000 

-.0016 -3.0000 .2600 0.0000 

10.7991 -9<59.0000 999.0000 0.0000 

6,6631 -999.0000 999.0000 0.0000 

•6568 -999.0000 999.0000 0.0000 

-.0613 -999.0000 999.0000 0.0000 

.0610 -999.00C0 999.0000 0*0000 

,0612 -999.0000 999,0000 0.0000 

-670.5609 -999,0000 999.0000 0.0000 

12.2270 -999.0000 999.0000 0*0000 

6.0696 -999.0000 990.0000 0.0000 

.0011 -999.0000 999,0000 0.0000 

-2.8326 -99O.C000 999,0000 0,0000 

-7.5609 -999,0000 999.OC00 0,0000 

,0096 -999.0000 999.0000 0.0000 

-6.7296 -999.0000 999.0000 0.0000 

12.6395 -999.0000 999,0000 0.0000 

,3106 -999.0000 999.0000 0.0000 

.6266 -999.0000 999.0000 0.0000 

0.0000 -1.0000 1.0000 0.0000 

O.cooo -1.0000 1.0000 0.0000 
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COHMOH OUtiP 





OESXGN- 

21A6.46060 

10.73378 

172.03261 

911.49267 

3.99136 


82?53«26214 

•^A692 

1 

X 

X 


I 

Z 

I 

X 

X 

1TER«* 

1 

0 

7 

2 

1 


1979 

0 

0 

0 

' 0 


PX» 

0.0000 

0.0000 

o.cooo 

0.0000 

c.oooo 



0.0000 

0.0000 

o.cooo 

X 

X 



I 

I 

I 

I 

I 


CDS* 

.018A39289 

#007789039 

•002919233 

10.502198924 

#032756961 

#002502042 

COSAP- 

.037644705 

•005683200 

•002194607 

9.056559926 

#189143434 

X 

w« 

21.200000000 

2.000000000 

•360000000 

#140000000 

5*000000000 



.960000000 

•850000000 

•300000000 

151.788142905 

15.092617897 



•367874441 

•537923713 

#454076207 

-.190000000 

#932323800 



•1.000000000 

-.150000000 

-#290000000 

-#120716547 

•025000000 


MX- 

•400000000 

•100000000 

•460000000 

60.316604382 

16.037019979 



2.188322980 

. 577390275 

•604725140 

.469975411 

.2COCOCOGO 



-.113905980 

-1.349290330 

OoOOOOOOOOO 

0.000000000 

#001076020 

- 


.060000000 

•360000000 ' 

. loGOCOOOOOO 

35.000000000 

o’occoecoo 



GX- 

3.141592650 

22#580000000 

1443.380COOOOO 

-99.0UOOOOOOC 

411OO.COOOOO0OO 

32CO.OCOOOOOOO 

•000464779 

57.295779000 

e#330000000 

3.150000000 

-5.250000000 

7500.0C000C000 

oTCOOGOOOO 

30.000000000 

.800000000 

2o20000C000 

.009396652 

-12.725000000 

4,000000000 

#636912543 

77.776023746 

3000.000000000 

.606073010 

275.000973653 

200.000000000 

•5300000CO 

1.761201049 

0.000000000 

16.667000000 

lo7129nOG41 

32.174CC0C00 

0870.000000000 

•090C0C000 

1.100904500 

352.761689100 


CG- 

•44692 

•44692 

#10109 

•18169 

.65000 

•IC003 

OERIVCR- 

6.599575550 

-.000309099 

•283323942 

16.056229399 

-•623274330 

o.ooooocooo 

•213493621 

-78.177918252 

#001015132 

-5.190192106 



0.000000000 

-26.746354637 

•603934363 

•000840005 

-3.924404329 


DERIVAP- 

4.664806036 

•000531027 

0.000000000 

•501471571 

13.219153412 

-17.971247917 

-#486502225 

0.000000000 

#536361547 

• 0529'30270 
-61.309434652 
•000840003 

•000650740 

-3.487362363 

-2o764004C04 


STOR- 

5.717865050 

4.864e06036 

-.027830200 

1 

3.890306526 

3.050093306 

-.017255523 

I 

#336700306 

.290192456 

#172003120 

X 

#939955222 

•546927754 

#205004590 

I 

#505263796 

#736076092 

1 

H- 


VT5- 

281429.737 

74815.053 

38668.400 

.750 

156200.104 

431.364 

42707.240 

449.029 

75429.489 

1537.650 

12450000.000 

6.669 

117150.138 

80.692 

13330000.000 

6.996 

2X9073.G04 

2177o044 

24420SoC35 

44350.421 


CST- 

2005.065 

19.690 

14027.960 

201429#737 

• 123 
31749.000 

• 096 
4538.000 

86.560 

X 


• 

-.15147188 .18662137 -.15771096 .052i 

.20313425 .16913646 

29537 -.48790586 
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APPENDIX IV - Procedure File used to Execute OPDOT on the 
Langley Research Center C^omputer System 


OPDOT,T7770 ,CM70000 . 

USER(820235N) 

CHARGE, 101264, LRC. 

GET,0PD0T1. 

GET,INPUT=SM10. (SMIO is 

FTN ( I=0PD0T1 , OPT=2 , R=0 ) 

ATTACH (FTNMLIB/UN=LIBRARY,NA) 
LDSET(PRESET=ZERO,LIB=FTNMLIB) (All 

LGO. 

REWIND (TAPE4) 

REIVIND(LGO) 

GET,PPB. 

PPB. 

REWIND (TAPE?) 

ATTACH (LR030SF/UN=LIBRARY) 
GET,ABS2290/UN=181500N. 

ABS2290,TAPE7. 

PLOT.CALPOST,ll 

CONT.//BLANK PAPER, LEROY .3 PEN, 

CONT. BLACK INK, MULTIPLE PLOT MODE// 
EXIT. 


RM 1174 ARBUCKLE/SLIWA 


static margin 10% case — APPENDIX II) 

variables set to zero since program is 

already operational) 

(PPB is the binary code for the 
OPDOT plotting preprocessor) 


(Reference 28) 
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APPENDIX V - KEY PROGRAM VARIABLES 


SUBROUTINE 

VARIABLE 

DESCRIPTION [UNITS, IF APPLICABLE] 

SII4PACT 

AMP 

Amplitude of sinusoid transformation Z to X domain 

(Main 

AVE 

Ave. of sinusoid transformation Z to X domain 

Program) 

CAYY 

Penalizing weight for constraint violation 


FACT 

= 1 if constraint is to he considered, = 0 othen-mse 


GNORM 

Constraint normalization — ave. of boundaries 


MC 

Nvtmher of constraint violations 


MINEQ 

Number of constraint functions 


NONEL 

Number of NELMIN optimizations 


NVAR 

Number of independent design variables 


REQMIN 

Convergence criteria 


SCF 

Scale factor to multiply performance index 


SL 

Lower constraint boundary 


STEP 

Initial optimizer step size 


SU 

Upper constraint boundary 


XBAR 

Independent design variable in Z domain 


XBARO 

Initial value for independent design 


XINEQ 

= 1 if constraint is violated, = 0 othen^rise 


XL 

Lowest allowable value of independent design variable 
in X domain 


XMIN 

Vector of optimum independent design variable in 
Z domain 


XU 

Upper value of independent design variable in 
X domain 


YNEWLO 

Optimimi perform_ance index from NELMIN 

NELMN 

IHO 

Vertex with highest performance index 


ILO 

Vertex with lowest performance index 


P 

Coordinates of simplex 


PBAR 

Centroid of simplex 


START 

Initial independent design variables 

FN 

ILINE 

= 1 output; = 0 no output 

SETUP 

COST 

Performance index 


IPR 

= 1 output; = 0 no output 


OBJ 

Augmented performance index 


PENT 

Penalty contribution to OBJ 


UNAUG 

Unaugmented performance index 

DOCOST 

ATT 

Attendant’s cost [$/hr] 


BLKHR 

Block hours of design mission [hr] 


CONTROL 

Cost of logistics control [$/hr] 


COSTHR 

Total cost per hour [$/hr] 


CREW 

Crew cost [$/hr] 


DELAY 

Delay cost [$/hr] 


DEPRE 

Depreciation cost [$/hr] 
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APPENDEX V - cont. 


SUBROUTINE 

VARIABLE 

DESCRIPTION [UNITS, IF APPLICABLE] 


FARROI 

Fare for 15^ ROI [$/pass .-mile] 


FCOST 

Fuel cost [$/hr] 


FEELAND 

Landing fees [$] 


FL 

Flight length (lift-off to touchdown) 


lECON 

Number of economy seats 


IFIRS 

Number of first class seats 


INSUR 

Insurance cost [$/hr] 


PASSPHR 

Passengers per hour 


PER 

Percent of total cost [^] 


PRICE 

Purchase price of airplane [$] 


PROFIT 

Profit of operations [$] 


REVHR 

Revenue hours 


REVYR 

Revenue years _ 


ROI 

Return on investment [year 


RPM 

Revenue passenger miles 


SERVICE 

Servicing cost [$] 


SPARES 

Spares cost [$] 


SUPPORT 

Support cost [ $ ] 


TAXRT 

Tax rate 


TONCAR 

Tons of cargo [tons] 


TOT 

Total operating cost [$/hr] 


XINVEST 

Investment cost [$] 


XIOC 

Indirect operating cost [$/hr] 


YRICFLT 

Year inflation factor 

WEIGHT 

ACTCON 

Weight of active control system [lbs] 


BAGGAGE 

Weight of passenger's baggage [lbs] 


CARGO 

Cargo weight [lbs] 


CREW 

Crew weight [lbs] 


DELWTO 

Difference between last gross weight and new gross 
weight [lbs] 


FUDGE 

Weight overrun fudge factor 


FUEL 

Mission fuel [gallons] 


LT 

Length of tail [ft] 


WOIT 

Number of weight iterations 


PASS 

Passenger weight [lbs] 


SHT 

Horizontal tail voliime coefficient 


SVT 

Vertical tail volume coefficient 


TRV 

Vertical tail taper ratio 


WECTL 

Weight of electrical system [lbs] 


WTAC 

Air conditioning weight [lbs] 


WTENG 

Weight of engine [lbs] 


WTFOOD 

Weight of food [lbs] 


WTFSYS 

Fuel system weight [lbs] 


WTFUEL 

Fuel weight [lbs] 


WTFURN 

Weight of furnishings [lbs] 
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appendix V - cont. 


SUBROUTINE 

VARIABLE 

DESCRIPTION [UNITS, IF APPLICABLE] 


WTFUSE 

Fuselage weight [ihs] 


WTHT 

Horizontal tail weight [ihs] 


WTINST 

Instrument weight [lbs] 


WTLG 

Landing gear weight [lbs] 


WTMISC 

Miscellaneous weight [lbs] 


WTOZ 

Weight of oxygen system [lbs] 


WTSRT 

Weight of engine starters [lbs] 


WTVT 

Weight of vertical tail [lbs] 


WTWIN 

Weight of windows [lbs] 


WTVfING 

Weight of wing [lbs] 

CRUALT 

ALF 

Standard atmospheric property a, reference 


ALT 

Cruise altitude [ft] 


ANS 

Vector of atmosphere as a function of altitude 


CLCR 

Cruise lift coefficient 


DALT 

Difference between old altitude guess and new [ft] 


DP 

A pressure [ibs/ft^] 


DRDH 

density gradient [slug/ft^] 


KNEW 

New altitude guess [ft] 


HOLD 

Old altitude guess [ft] 


IC 

Number of iterations in table look up 


M 

Mach number 


P 

Pressure [ibs/ft^] 


PO 

Sea level pressure [ibs/ft^] 


R 

Gas constant 


TO 

Sea level temperature [®R] 


WT 

Cruise weight [lbs] 

ENGIflE 

ALT 

Operating altitude [ft] 


DELTA 

Pressure ratio, P/Pq 


F 

Fuel flow rate [lbs /hr] 


FN 

Normalized fuel flow 


HN 

Normalized altitude 


TCTM 

Cruise thrust over installed thrust 


THETA 

Temperature ratio, T/Tq 


TN 

Normalized thrust 


TSFC 

Thrust specific fuel consumption 

CDZL 

ABOD 

Surface area of fuselage body [ft^] 


ACONE 

Surface area of tail cone [ft^] 


ANOSE 

Surface a.rea of nose cone [ft^] 


BV 

Span of vertical tail [ft] 


CDBF 

Fuselage bluff body drag coefficient 


CDFF 

Fuselage friction drag coefficient 


CDOE 

Engine drag coefficient 


CDOF 

Fuselage drag coefficient 
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APPENDIX V - cont. 


SUBROUTINE 

VARIABLE 

DESCRIPTION [UNITS, IF APPLICABLE] 


CDOH 

Horizontal tail drag coefficient 


CDOV 

Vertical tail drag coefficient 


CDOW 

Wing drag coefficient 


CF 

Friction coefficient 


CONEL 

Length of tail cone [ft] 


CVBAR 

Vertical tail mean aerodynamic chord [ft] 


DIA 

Fuselage diameter [ft] 


EDIA 

Engine diameter [ft] 


EL 

Engine length [ft] 


FLOD 

Ratio of fuselage length to diameter 


MU 

Viscosity [sec~2] 


SREF 

Reference area [ft 2] 

STABCOD 

AW 

Wing lift curve slope [deg~^] 


BETA 

Mach number correction factor 


CLAWB 

of wing body 


DCDM 

9Cd/9M 


DCMDCL 

Cm/9CL, static stability 


DEDA 

Downwash gradient, de/da 9 ^ 



DEDAMO 

Incompressible downwash gradient, ^ 

M=0 


DXACE 

Shift in aerodynamic center due to engines 


ELASTK 

Elasticity correction factor 


ETAH 

Tail efficiency, 


XAC 

Aerodynamic center 


XACDM 

Shift in aerodynamic center due to compressibility 


XACW 

Wing aerodynamic center 

TRIM 

ALFAT 

Angle-of-attack of tail [deg] 


ALFAW 

Angle-of-attack of wing [deg] 


AOL 

Angle-of- zero -lift [deg] 


AT 

Tail lift curve slope [ deg ^ ] 


AW 

Wing lift curve slope [deg 


CLTAIL 

Tail lift coefficient 


CLWING 

Wing lift coefficient 


CMACW 

Wing pitching moment coefficient 


CMFUS 

Fuselage pitching moment coefficient 


CNT 

Normal tail force coefficient 


CNW 

Normal wing force coefficient 


CRDE 

Elevator deflection [deg] 


GRIT 

Stabilizer deflection [deg] 


CTLPWR 

Control power coefficient 


DEDA 

Downwash gradient 


EPS 

Downwash angle [deg] 


LTLOQ 

Tail lift normalized by dynamic pressure 


LTOTOQ 

Total lift normalized by dynamic pressure 


LWOQ 

Wing lift normalized by dynamic pressure 
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APPENDIX V - cont. 


SUBROUTINE 

VARIABLE 

DESCRIPTION [UNITS, IF APPLICABLE] 


TAO 

f 

Elevator control effectiveness parameter 

FUELCAL 

BLKSPD 

Block speed for mission [knots] 


RANGE 

Range of mission [naut-mi] 


WTMID 

Mid-point weight [lbs] 


WIRES 

Weight after reserve mission [lbs] 


XMG 

Miles per gallon of fuel 


XSMG 

Seat miles per gallon of fuel 

CRUFUEL 

ALTCR 

Cruise altitude [ft] 


ALTD 

Desired altitude [ft] 


ALTEND 

Altitude at end of leg [ft] 


ALTLDM 

Alt5.tude for maximiim glide ratio '[ft] 


CLCR 

Cruise lift coefficient 


CLM 

Lift coefficient for maximum glide ratio 


GAMA 

Climb gradient 


IRES 

=1, reserve leg, = 0 otherwise 


LOD 

Glide ratio 


ROC 

Rate-of-climb , [ft/hr] 


TOWAV 

T/W available 


TOWRQ 

T/W required 

XLOD 

CLLDMX 

Lift coefficient for maximum glide rs.tio 


CLO 

Design lift coefficient of wing section 


CLWI 

Wing lift coefficient normal to leading edge 


EPS 

Wing doimwash angle at tail [deg] 


SIGI-IA 

Wing-ta.il interference factor 

i 

TDG 

Tail drag contribution 


WMCR 

Critical Mach number of \/ing 


XLDI4X 

Maairaum glide ratio 

MAINCST 

lcost(i) 

Labor cost of Ith component [$/hr] 


MCOST(I) 

Material cost of Ith component [$/hr] 

CN3TRN 

AAP 

Coefficients of longitudinal characteristic 



polynomial in approach 


ACR 

Coefficients of longitudinal characteristic 



polynomial in cruise 


CLA 

Approach lift coefficient 


CLS 

Stall lift coefficient with fiai flaps 


CLTAIL 

Required tail lift coefficient 


CLTAV 

Available tail lift coefficient for nose gea,r imstir 


CL2 

Second segment climb lift coefficient 


DERAP 

Dim.ensional stability derivatives in approach 


DERCR 

Dimensional stability derivatives in cruise 


DGDU 
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APPENDIX V - cont. 


SUBROUTINE 

VARIABLE 

DESCRIPTION [UNITS, IF APPLICABLE] 


DZETA 

ROOTAP 

ROOTCR 

UNTHRU 

Desired short period damping ratio 
Longitudinal roots in approach 
Longitudinal roots in cruise 

Unbalanced thrust component during take-off roll 


APPENDIX VI - MAP OF COMMON BLOCKS USED WITHIN 
DESIGN SECTION OF PROGRAM 


COI^G^ON 

BLOCK 

ARRAY 

NO. 

DESCRIPTION [UNITS, IF APPLICABLE] 

DEVAR 

CST 

1 

Direct operating cost per block hour [$/hr] 



2 

Direct operating cost per flight [$/flight]. 
Return on investiment [year ^] 



3 



1* 

Fare for a ROI of .15 [$/seat-naut. mi] 



5 

Fuel efficiency [seat-naut. mi/gal] 


1 

6 

Maximum glide ratio in cruise 



T 

Maocimum take-off gross wt. [lbs] 


DESIGN 

1 

* Wing area [ft^] 



2 

Wing 3,spect ratio 



3 

Fuselage length [ft] 



h 

[ Horizontal tail area [ft^] 



5 

Horizontal tail aspect ratio 



6 

Installed thrust [lbs] 



7 

Aft most center-of-gravity in cruise [% tiAC] 


ITERM 

1 

= 1 if active controls, = 0 otherwise 



2 

= 1 if center-of-gravity control, = 0 othervrise 



3 

Cruise C.G. Selector: = 7 if midway betvreen FliTD 

& AFT limits 



4 

Number of engines 



5 

Element of CST used in optimizer 



6 

Year of evaluation 



7 

= 1 if landing gear is movable, = 0 otherwise 

DRAG 

CDS 

1 

Cruise Total parasite drag coefficient 



2 

Cruise Wing drag coefficient 



3 

Cruise Horizontal tail drag coefficient 



4 

Cruise Glide ratio 



5 

Cruise Total drag coefficient 



6 

Cruise Fuselage drag coefficient 
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APPENDIX VI - cont. 


COMMON 

BLOCK 

ARRAY 

NO. 

DESCRIPTION [UNITS, IF APPLICABLE] 


CDSAP 

1 

Approach Total parasite drag coefficient with 




take-off flaps 



2 

Approach Wing drag coefficient 



3 

Approach Horizontal tail drag coefficient 



h 

Approach Glide ratio with landing flaps 



5 

Approach Total drag coefficient 



6 

Approach Fuselage drag coefficient 

GEOM 

GX 

1 

PI (tt) = 3 .IU 15927 



2 

Radian to degree conversion factor = 57.295779 



3 

Cruise Ma-ch numher 



U 

Design range [naut. miles] 



5 

Fuselage diameter [ft] 



6 

Reference engine length [ft] 



T 

Reference engine diameter [ft] 



0 

Maximum wing lift coefficient with flaps in take- 




off position 



9 

Trimmed lift coefficient in cruise 



10 

Trammed lift coefficient 



11 

Characteristic turbulence length divided by /s 




[ft] 



12 

Maximum wing lift coefficient with full flaps 



13 

Variance of elevator deflection in approa.ch [ra.d] 



iH 

Velocity in approach [ft/ sec ] 



15 

Acceleration constant, g [ft/sec^] 



l6 

Take-off stabilizer position, [deg], if -99 then 




use trim position for climbout 



IT 

Angle of zero lift in take-off [deg] 



18 

Angle of zero lift with landing flaps [d.eg] 



19 

Number of pa.ssenger seats 



20 

Reference engine weight [lbs] 



21 

Thrust of reference engine [lbs] 



22 

Weight of cargo [lbs] 



23 

Allowable C.G. range [% MAC] 



2h 

Pa.ssenger load factor 



25 

Fare [ $/ seat-mile ] 



26 

Airplane utilization [hrs/yr] 



27 

Ratio of maximum landing weight to gross weight 



28 

Wing lift coefficient in cruise 



29 

VJing lift coefficient in approach 



30 

Minimimi ratio of thrust available to thrust 




required in cruise 



31 

Trim drag correction for estimating L/Dj^g^ 



32 

Sweep of horizontal tail [deg] 



33 

Distance between horizontal tail and wing [ft] 
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APPENDIX VI - cont. 




COMMON 

BLOCK 

ARRAY 

NO. 

DESCRIPTION [UNITS, IF APPLICABLE] 



34 

Height of horizontal tail above fuselage 




centerline [ft] 



35 

Vertical tail area [ft^] 


HX 

1 

Horizontal ta.il taper ratio 



2 

Horizontal tail thickness ratio 



3 

Elevator control effectiveness factor 



1 + 

Horizontal tail span [ft] 



5 

Horizontal tail mean aerodynamic chord [ft] 



6 

Horisonta,! tail volume coefficient 



7 

Tan (tail qua.rter chord sweep angle) 



8 

Tan (taJ.l 3.eading edge sweep angle) 


1 

9 

Tan (tail half chord sweep angle) 



10 

Elevator servo time constant [sec] 



11 

Approach tail lift coefficient 



12 

Tail incidence in approach [deg] 



13 

Trimmed elevator position in approach [deg] 



lU 

Trimmed elevator position in cruise [deg] 



15 

Va.ria.nce of elevator in cruise [rad] 



16 

Vertical tail volume coefficient 



17 

Vertical tail taper ratio 



18 

Vertical tail aspect ratio 



19 

Vertical oall sweep at quarter chord [deg] 



20 

Ratio of rudder area to vertical tail area 


w 

1 

Wing sweep at quarter chord [deg] 



2 

Incidence of wing [deg] 



3 

Wing taper ratio 



k 

Wing thickness ratio 



5 

Wing geometric twist [deg] 



6 

Oswald’s efficiency factor in drag bucket 



7 

Osward’s efficiency factor out of drag bucket 



8 

Design lift coefficient for wing section 



9 

Wing span [ft] 



10 

Wing mean aerodynamic chord, MAC [ft] 



11 

Tan (wing quarter chord sweep angle) 



12 

Tan (wing leading edge sweep angle) 



13 

Tan (wing half chord sweep angle) 



1I+ 

Wing pitching moment coefficient at cruise 



15 

Cos (wing quarter chord sweep angle) 



16 

Wing angle-of-zero-lift with no flaps [deg] 



17 

Increase in wing pitching moment coefficient at 




take-off 



18 

Increase in wing pitching moment coefficient at 




landing 
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APPENDIX VI - cont. 


COMMON 

BLOCK 

ARRAY 

NO. 

? DESCRIPTION [UNITS, IF APPLICABLE] 



19 

Wing pitching moment coefficient at a^pproach 



20 

Increase in drag due to flap deflection 


PX 

1 

Pitching moment reduction coefficient 



2 

TSFC reduction coefficient 



3 

Weight reduction coefficient 



h 

Drag reduction coefficient 



5 

Disconnect rubber engine (l.O = yes; 0.0 = no) 



6 

Maintenance cost boost {% Boost = PX(6) ** ^%) 



7 

Purchase price boost {% Boost = PX(T) 5%) 



8 

Maneuver loaxl alleviation (Reduction in design 




3-imit load in g’s) 

GRAVITY 

CGS 

1 

Aft cruice C.G. position [% RAC] 



2 

Aft approach C.G- position [% MAC] 



3 

Fon-rard cruise C.G. position [% MAC] 



1+ 

ForTrard approach C.G. powition [%¥AC] 



5 

Landing gear position [% MAC] 



6 

Required distance between C.G. and landing gear 




MAC] 

STAB 

DERIVAP 

1 

Approach 

JLi 

a 

S 

a 



2 



3 

U 

a 



h 

"l 

U 



5 

S 

u 



6 

S/r 




u 



7 

^L 



8 

s 

q 



■9 

q 



10 

"l. 

a 



11 

a 



12 

a 
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APPENDIX VI - cont 


COMMON 

BLOCK 

ARRAY 

NO. 

DESCRIPTION [UNITS, IF APPLICABLE] 



13 




lU 




15 

X 


DERIVCR 

1 

Cruise C^^ 




a 



2 

S 

a 



3 

M 

a 



k 

^L 

u 



5 

s 

u 



6 

"m 

u 



7 

"l 



8 

S 

<1 



9 

c,, 

M 

q 



10 

"l. 

a 



11 

S. 

a 



12 

Cm 

M. 

a 



13 

"^5 



lU 




15 



STOR 

1 

Cruise C^^ 

a wing 



2 

c. 

a horizontal tail 



3 

9e/8a 



h 

Stick fixed neutral point \.% MAC] 



5 

Stick fixed maneuver point [% MAC] 
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APPENDIX VI - cont . 


COMMON 

BLOCK 

ARRAY 

NO. 

DESCRIPTION [UNITS, IF APPLICABLE] 



6 

Approach Cj^ 




a wing 



7 

Cl 




a horizontal tail 



8 

9e/9a 



9 

Stick fixed neutral point [% MAC] 



10 

stick fixed maneuver point \% M<^C] 



11 

C fuselage 




mac 



12 

C engine 




mac 

OTSVE 

WTS 

1 

Maximum gross weight at take-off [ihs] 



2 

Empty weight [lbs] 



3 

Fuel f].ow rate [ibs/block hour] 



h 

Manufacturers airframe weight [lbs] 



5 

Landing weight after mission [lbs] 



6 

Fuel weight including reserves [lbs] ^ 



7 

Aircraft specific density (cruise), 



8 

Radius of gyration squared (cruise), ly/m [ft] 


1 

9 ^ 

Aircraft specific density (approach), 



10 

Radius of gyration squares (approach), ly/m [ft] 



11 

AJ-titude at mid cruise [ft] 



12 

L/D at cruise altitude 




max 



13 

Pi.tch moment of inertia (cruise), [slug- ft] 



2.h 

Pitch moment of inertia (approach), [s3.ug-ft] 



15 

Weight at mid cruise [lbs] 



l6 

Fuel cost [$/gal] 



17 

Cruise velocity [knots] 



18 

Flight tim.e for mission [hrs] 



19 

Block time [hrs] 



20 

Weight of fuel to fly economic mission [3.bs] 


99 






^ • . 


APPENDIX VII - ASSUMPTIONS USED IN CALCULATING 
TRANSPORT DESIGN FACTORS 


MISSION; 


Cruise Mach Number .80 

Divergence Mach Number .8L 

Design Range 65 OO km 

Number of Seats 200 

Cargo 33^00 N 

Maximiam Lift Coefficient 3.15 

Landing Field Requirement 2 UU 0 m 

Take-Off Field Requirement 3050 m 


GEOIffiTRY: 


Wing Sweep Angle 26.4 deg 

Wing Thickness Ratio .12 

Wing Taper Ratio .38 

Wing Incidence Angle 2 deg 

Wing Geometric Twist 5 deg 

Tail Thickness Ratio .10 

Tail Sweep Angle 30 deg 

Tail Taper Ratio .4 

Vertical Tail Sweep 35 deg 

Ratio of Rudder Area to Vertical Tail Area .30 

Ratio of Elevator Chord to Horizontal Tail Chord .25 
Ratio of Flap Span to Wing Span .6 

Maximum Flap Deflection 45 deg 

Fuselage Diameter 5.08 m 

Height of Aerodynamic Center Above c.g. .08 MAC 

Height of Thrust Vector Above c.g. -.12 liAC 

Height of Horizontal Tail Above c.g. 0 

Number of Engines 2 


ECONOMICS: 

Fuel Cost 
Load Factor 
Fare 

Utilization Rate 
Depreciation Period 
Residual Value 
Tax Rate 
Year of Study 

Assumed Annual Inflation Rate 
Number of Prototype Aircraft 


20(^/liter 

.55 

9<^/seat-naut . mi 
3200 hr/yr 
l4 yr 

12 percent 

.48 

1979 

.07 

2 
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APPENDIX VII - cont. 


Aircraft Fleet Size 

250 


Initial Production Rate 

. 5 /month 


Full Production Rate 

5 /month 

$/hr 

Engineering Rate 

19.55 

Tooling Rate 

lit. 00 '7lt 

$/hr 

Labor Rate 

10.90 '7l+ 

$/hr 

Engines for Test Aircraft 

Ratio of Manufacturer's Airframe Weight to 

3 


Take-Off Wt. 

.75 



MISCELLANEOUS : 

Maximum Dynamic Pressure 

5.13 N/m^ 

Pressurized Volume 

178.2 m3 

Number of Pilots 

3 

Number of Attendants 

8 

Air Conditioning Flow Rate 

200 kg/min 

Autopilot Channels (w/MLIX) 

5 

General Capacity 

750 kilovolt-amperes 

Maintenance Complexity Factor 

1.6 

Hydraulics Volume Flow Rate 

300 liters /min 

Number of Inertial Platform Systems 

1 

Ratio of APU-on Time to Engine on Time 

.1 

Curved Windshield 

Ratio of First Class to Economy Seating 

.15 

Maximum Speed 

U83 knots 

Supercritical Airfoil Technology 
Airfoil Design Lift Coefficient 

.5 

Some Nonlinear Aerodynamics Terms 
Baseline Engine 

CF-6 

Elevator Servo Time Consto-nt 

.1 sec 
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OPTIMAL DESIGN METHODOLOGY 



ECONOMICS 
FIXED GEOMETRIES 
TECHNOLOGY LEVEL 
NON-LINEAR AERO 
BASELINE ENGINE 
ETC. 


Figure 1.- Generalized flow diagram for OPDOT 





Figure 2.- Schematic showing primary calling sequence of subroutines used to 
evaluate the performance index and constraint functions. 



















MISSION PROFILE 



Figure 3.- Mission profile used in OPDOT 





MIL 8785B LEVEL III (SIMPflCT) 
X Z -45.0 10.0-20.0 


10. OORT 


A 


Figure U.- Concluded. 
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